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.

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.

Identification of the semaphorin receptor PLXNA2 as a candidate for susceptibility to schizophrenia

The discovery of genetic factors that contribute to schizophrenia susceptibility is a key challenge in understanding the etiology of this disease. Here, we report the identification of a novel schizophrenia candidate gene on chromosome 1q32, plexin A2 (PLXNA2), in a genome-wide association study using 320 patients with schizophrenia of European descent and 325 matched controls. Over 25,000 single-nucleotide polymorphisms (SNPs) located within approximately 14,000 genes were tested. Out of 62 markers found to be associated with disease status, the most consistent finding was observed for a candidate locus on chromosome 1q32. The marker SNP rs752016 showed suggestive association with schizophrenia (odds ratio (OR) = 1.49, P = 0.006). This result was confirmed in an independent case-control sample of European Americans (combined OR = 1.38, P = 0.035) and similar genetic effects were observed in smaller subsets of Latin Americans (OR = 1.26) and Asian Americans (OR = 1.37). Supporting evidence was also obtained from two family-based collections, one of which reached statistical significance (OR = 2.2, P = 0.02). High-density SNP mapping showed that the region of association spans approximately 60 kb of the PLXNA2 gene. Eight out of 14 SNPs genotyped showed statistically significant differences between cases and controls. These results are in accordance with previous genetic findings that identified chromosome 1q32 as a candidate region for schizophrenia. PLXNA2 is a member of the transmembrane semaphorin receptor family that is involved in axonal guidance during development and may modulate neuronal plasticity and regeneration. The PLXNA2 ligand semaphorin 3A has been shown to be upregulated in the cerebellum of individuals with schizophrenia. These observations, together with the genetic results, make PLXNA2 a likely candidate for the 1q32 schizophrenia susceptibility locus.

Association study of the dystrobrevin-binding gene with schizophrenia in Australian and Indian samples

Numerous studies have reported association between variants in the dystrobrevin binding protein 1 (dysbindin) gene (DTNBP1) and schizophrenia. However, the pattern of results is complex and to date, no specific risk marker or haplotype has been consistently identified. The number of single nucleotide polymorphisms (SNPs) tested in these studies has ranged from 5 to 20. We attempted to replicate previous findings by testing 16 SNPs in samples of 41 Australian pedigrees, 194 Australian cases and 180 controls, and 197 Indian pedigrees. No globally significant evidence for association was observed in any sample, despite power calculations indicating sufficient power to replicate several previous findings. Possible explanations for our results include sample differences in background linkage disequilibrium and/or risk allele effect size, the presence of multiple risk alleles upon different haplotypes, or the presence of a single risk allele upon multiple haplotypes. Some previous associations may also represent false positives. Examination of Caucasian HapMap phase II genotype data spanning the DTNBP1 region indicates upwards of 40 SNPs are required to satisfactorily assess all nonredundant variation within DTNBP1 and its potential regulatory regions for association with schizophrenia. More comprehensive studies in multiple samples will be required to determine whether specific DTNBP1 variants function as risk factors for schizophrenia.

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.

Separate and interacting effects within the catechol-O-methyltransferase (COMT) are associated with schizophrenia

Several lines of evidence have implicated the catechol-O-methyltransferase (COMT) gene as a candidate for schizophrenia (SZ) susceptibility, not only because it encodes a key dopamine catabolic enzyme but also because it maps to the velocardiofacial syndrome region of chromosome 22q11 which has long been associated with SZ predisposition. The interest in COMT as a candidate SZ risk factor has led to numerous case-control and family-based studies, with the majority placing emphasis on examining a functional Val/Met polymorphism within this enzyme. Unfortunately, these studies have continually produced conflicting results. To assess the genetic contribution of other COMT variants to SZ susceptibility, we investigated three single-nucleotide polymorphisms (SNPs) (rs737865, rs4633, rs165599) in addition to the Val/Met variant (rs4680) in a highly selected sample of Australian Caucasian families containing 107 patients with SZ. The Val/Met and rs4633 variants showed nominally significant associations with SZ (P<0.05), although neither of the individual SNPs remained significant after adjusting for multiple testing (most significant P=0.1174). However, haplotype analyses showed strong evidence of an association; the most significant being the three-marker haplotype rs737865-rs4680-rs165599 (global P=0.0022), which spans more than 26 kb. Importantly, conditional analyses indicated the presence of two separate and interacting effects within this haplotype, irrespective of gender. In addition, our results indicate the Val/Met polymorphism is not disease-causing and is simply in strong linkage disequilibrium with a causative effect, which interacts with another as yet unidentified variant approximately 20 kb away. These results may help explain the inconsistent results reported on the Val/Met polymorphism and have important implications for future investigations into the role of COMT in SZ susceptibility.

The use of genetic correlations to evaluate associations between SNP markers and quantitative traits

Open-pollinated progeny of Corymbia citriodora established in replicated field trials were assessed for stem diameter, wood density, and pulp yield prior to genotyping single nucleotide polymorphisms (SNP) and testing the significance of associations between markers and assessment traits. Multiple individuals within each family were genotyped and phenotyped, which facilitated a comparison of standard association testing methods and an alternative method developed to relate markers to additive genetic effects. Narrow-sense heritability estimates indicated there was significant additive genetic variance within this population for assessment traits ( h ˆ 2 =0.28to0.44 ) and genetic correlations between the three traits were negligible to moderate (r G = 0.08 to 0.50). The significance of association tests (p values) were compared for four different analyses based on two different approaches: (1) two software packages were used to fit standard univariate mixed models that include SNP-fixed effects, (2) bivariate and multivariate mixed models including each SNP as an additional selection trait were used. Within either the univariate or multivariate approach, correlations between the tests of significance approached +1; however, correspondence between the two approaches was less strong, although between-approach correlations remained significantly positive. Similar SNP markers would be selected using multivariate analyses and standard marker-trait association methods, where the former facilitates integration into the existing genetic analysis systems of applied breeding programs and may be used with either single markers or indices of markers created with genomic selection processes.

New SNPs for population genetic analysis reveal possible cryptic speciation of eastern Australian sea mullet (Mugil cephalus)

Sustainable management of sea mullet (Mugil cephalus) fisheries needs to account for recent observations of regional-scale differentiation. Population genetic analysis is sought to assess the situation of this ecologically and economically important fish species in eastern Australian waters. Here, we report (i) new population genetic markers [single nucleotide polymorphisms (SNPs) and potential microsatellites], (ii) first estimates of spatial genetic differentiation and (iii) prospective power tests for designing more comprehensive studies. Six DNA samples from three sampling regions (North Queensland, South Queensland and central New South Wales) on the eastern coast of Australia were used to prepare restriction site associated DNA (RAD) tag libraries from genomic DNA digested with EcoRI and MseI. A pooled sample of regional RAD tag libraries was sequenced using the Roche GS-FLX Titanium platform. A total of 172837 raw reads (17.4Mbp) were retrieved, 95500 of which were used to discover 1267 SNPs and 1417 microsatellites. A subset of 161 SNPs was validated based on 63 additional DNA samples genotyped using the Sequenom MassArray (iPLEX Gold chemistry). Altogether 92 SNPs (57%) were confirmed, with 40% of these marking fixed variants between northern and southern sampling regions. Our preliminary findings indicate a multispecies fishery stock of M. cephalus in eastern Australian waters, but suggest that strong genetic differentiation occurs north of major fishing grounds. Low potential differentiation within major fishing grounds (e.g. FST=0.0025) can be resolved with a likely power 67% by using standard sample sizes of 50 and validated subsets of available markers.

Smoking behaviour modifies IL23r-associated disease risk in patients with Crohn's disease

Background and Aim The etiology of Crohn's disease (CD) implicates both genetic and environmental factors. Smoking behavior is one environmental risk factor to play a role in the development of CD. The study aimed to assess the contribution of the interleukin 23 receptor (IL23R) in determining disease susceptibility in two independent cohorts of CD, and to investigate the interactions between IL23R variants, smoking behavior, and CD-associated genes, NOD2 and ATG16L1. Methods Ten IL23R single-nucleotide polymorphisms (SNPs) were genotyped in 675 CD cases, and 1255 controls from Brisbane, Australia (dataset 1). Six of these SNPs were genotyped in 318 CD cases and 533 controls from Canterbury, New Zealand (dataset 2). Case–control analysis of genotype and allele frequencies, and haplotype analysis for all SNPs was conducted. Results We demonstrate a strong increased CD risk for smokers in both datasets (odds ratio 3.77, 95% confidence interval 2.88–4.94), and an additive interaction between IL23R SNPs and cigarette smoking. Ileal involvement was a consistent marker of strong SNP–CD association (P ≤ 0.001), while the lowest minor allele frequencies for location were found in those with colonic CD (L2). Three haplotype blocks were identified across the 10 IL23R SNPs conferring different risk of CD. Haplotypes conferred no further risk of CD when compared with single SNP analyses. Conclusion IL23R gene variants determine CD susceptibility in the Australian and New Zealand population, particularly ileal CD. A strong additive interaction exists between IL23R SNPs and smoking behavior resulting in a dramatic increase in disease risk depending upon specific genetic background.

Sequence diversity and differential expression of major phenylpropanoid-flavonoid biosynthetic genes among three mango varieties

Background: Mango fruits contain a broad spectrum of phenolic compounds which impart potential health benefits; their biosynthesis is catalysed by enzymes in the phenylpropanoid-flavonoid (PF) pathway. The aim of this study was to reveal the variability in genes involved in the PF pathway in three different mango varieties Mangifera indica L., a member of the family Anacardiaceae: Kensington Pride (KP), Irwin (IW) and Nam Doc Mai (NDM) and to determine associations with gene expression and mango flavonoid profiles. Results: A close evolutionary relationship between mango genes and those from the woody species poplar of the Salicaceae family (Populus trichocarpa) and grape of the Vitaceae family (Vitis vinifera), was revealed through phylogenetic analysis of PF pathway genes. We discovered 145 SNPs in total within coding sequences with an average frequency of one SNP every 316bp. Variety IW had the highest SNP frequency (one SNP every 258bp) while KP and NDM had similar frequencies (one SNP every 369bp and 360bp, respectively). The position in the PF pathway appeared to influence the extent of genetic diversity of the encoded enzymes. The entry point enzymes phenylalanine lyase (PAL), cinnamate 4-mono-oxygenase (C4H) and chalcone synthase (CHS) had low levels of SNP diversity in their coding sequences, whereas anthocyanidin reductase (ANR) showed the highest SNP frequency followed by flavonoid 3'-hydroxylase (F3'H). Quantitative PCR revealed characteristic patterns of gene expression that differed between mango peel and flesh, and between varieties. Conclusions: The combination of mango expressed sequence tags and availability of well-established reference PF biosynthetic genes from other plant species allowed the identification of coding sequences of genes that may lead to the formation of important flavonoid compounds in mango fruits and facilitated characterisation of single nucleotide polymorphisms between varieties. We discovered an association between the extent of sequence variation and position in the pathway for up-stream genes. The high expression of PAL, C4H and CHS genes in mango peel compared to flesh is associated with high amounts of total phenolic contents in peels, which suggest that these genes have an influence on total flavonoid levels in mango fruit peel and flesh. In addition, the particularly high expression levels of ANR in KP and NDM peels compared to IW peel and the significant accumulation of its product epicatechin gallate (ECG) in those extracts reflects the rate-limiting role of ANR on ECG biosynthesis in mango. © 2015 Hoang et al.

Genomic deletions and mutations resulting in the loss of eight genes reduce the in vivo replication capacity of Meleagrid herpesvirus 1

Meleagrid herpesvirus 1 (MeHV-1 or turkey herpesvirus) has been widely used as a vaccine in commercial poultry. Initially, these vaccine applications were for the prevention of Marek’s disease resulting from Gallid herpesvirus 2 infections, while more recently MeHV-1 has been used as recombinant vector for other poultry infections. The construction of herpesvirus infectious clones that permit propagation and manipulation of the viral genome in bacterial hosts has advanced the studies of herpesviral genetics. The current study reports the construction of five MeHV-1 infectious clones. The in vitro properties of viruses recovered from these clones were indistinguishable from the parental MeHV-1. In contrast, the rescued MeHV-1 viruses were significantly attenuated when used in vivo. Complete sequencing of the infectious clones identified the absence of two regions of the MeHV-1 genome compared to the MeHV-1 reference sequence. These analyses determined the rescued viruses have seven genes, UL43, UL44, UL45, UL56, HVT071, sorf3 and US2 either partially or completely deleted. In addition, single nucleotide polymorphisms were identified in all clones compared with the MeHV-1 reference sequence. As a consequence of one of the polymorphisms identified in the UL13 gene, four of the rescued viruses were predicted to encode a serine/threonine protein kinase lacking two of three domains required for activity. Thus four of the recovered viruses have a total of eight missing or defective genes. The implications of these findings in the context of herpesvirus biology and infectious clone construction are discussed.

A transcriptomic analysis of the kidney tissue of Tra catfish (pangasianodon hypophthalmus) reared in saline condition: De novo assembly, annotation, SNP discovery

Pangasianodon hypophthalmus is a commercially important freshwater fish used in inland aquaculture in the Mekong Delta, Vietnam. The current study using Ion Torrent technology generated EST resources from the kidney for Tra catfish reared at a salinity level of 9 ppt. We obtained 2,623,929 reads after trimming and processing with an average length of 104 bp. De novo assemblies were generated using CLC Genomic Workbench, Trinity and Velvet/Oases with the best overall contig performance resulting from the CLC assembly. De novo assembly using CLC yielded 29,940 contigs, and allowing identification of 5,710 putative genes when comppared with NCBI non-redundant database. A large number of single nucleotide polymorphisms (SNPs) were also detected. The sequence collection generated in our study represents the most comprehensive transcriptomic resource for P. hypophthalmus available to date.

Genomic deletions and mutations resulting in the loss of eight genes reduce the in vivo replication capacity of Meleagrid herpesvirus 1

Meleagrid herpesvirus 1 (MeHV-1 or turkey herpesvirus) has been widely used as a vaccine in commercial poultry. Initially, these vaccine applications were for the prevention of Marek’s disease resulting from Gallid herpesvirus 2 infections, while more recently MeHV-1 has been used as recombinant vector for other poultry infections. The construction of herpesvirus infectious clones that permit propagation and manipulation of the viral genome in bacterial hosts has advanced the studies of herpesviral genetics. The current study reports the construction of five MeHV-1 infectious clones. The in vitro properties of viruses recovered from these clones were indistinguishable from the parental MeHV-1. In contrast, the rescued MeHV-1 viruses were significantly attenuated when used in vivo. Complete sequencing of the infectious clones identified the absence of two regions of the MeHV-1 genome compared to the MeHV-1 reference sequence. These analyses determined the rescued viruses have seven genes, UL43, UL44, UL45, UL56, HVT071, sorf3 and US2 either partially or completely deleted. In addition, single nucleotide polymorphisms were identified in all clones compared with the MeHV-1 reference sequence. As a consequence of one of the polymorphisms identified in the UL13 gene, four of the rescued viruses were predicted to encode a serine/threonine protein kinase lacking two of three domains required for activity. Thus four of the recovered viruses have a total of eight missing or defective genes. The implications of these findings in the context of herpesvirus biology and infectious clone construction are discussed.

Computational Methods for Detecting Large-Scale Chromosome Rearrangements in SNP Data

Large-scale chromosome rearrangements such as copy number variants (CNVs) and inversions encompass a considerable proportion of the genetic variation between human individuals. In a number of cases, they have been closely linked with various inheritable diseases. Single-nucleotide polymorphisms (SNPs) are another large part of the genetic variance between individuals. They are also typically abundant and their measuring is straightforward and cheap. This thesis presents computational means of using SNPs to detect the presence of inversions and deletions, a particular variety of CNVs. Technically, the inversion-detection algorithm detects the suppressed recombination rate between inverted and non-inverted haplotype populations whereas the deletion-detection algorithm uses the EM-algorithm to estimate the haplotype frequencies of a window with and without a deletion haplotype. As a contribution to population biology, a coalescent simulator for simulating inversion polymorphisms has been developed. Coalescent simulation is a backward-in-time method of modelling population ancestry. Technically, the simulator also models multiple crossovers by using the Counting model as the chiasma interference model. Finally, this thesis includes an experimental section. The aforementioned methods were tested on synthetic data to evaluate their power and specificity. They were also applied to the HapMap Phase II and Phase III data sets, yielding a number of candidates for previously unknown inversions, deletions and also correctly detecting known such rearrangements.

GPRA and the asthma locus on chromosome 7p14-p15

The basis of this work was the identification of a genomic region on chromosome 7p14-p15 that strongly associated with asthma and high serum total immunoglobulin E in a Finnish founder population from Kainuu. Using a hierarchical genotyping approach the linkage region was narrowed down until an evolutionary collectively inherited 133-kb haplotype block was discovered. The results were confirmed in two independent data sets: Asthma families from Quebec and allergy families from North-Karelia. In all the three cohorts studied, single nucleotide polymorphisms tagging seven common gene variants (haplotypes) were identified. Over half of the asthma patients carried three evolutionary closely related susceptibility haplotypes as opposed to approximately one third of the healthy controls. The risk effects of the gene variants varied from 1.4 to 2.5. In the disease-associated region, there was one protein-coding gene named GPRA (G Protein-coupled Receptor for Asthma susceptibility also known as NPSR1) which displayed extensive alternative splicing. Only the two isoforms with distinct intracellular tail sequences, GPRA-A and -B, encoded a full-length G protein-coupled receptor with seven transmembrane regions. Using various techniques, we showed that GPRA is expressed in multiple mucosal surfaces including epithelial cells throughout the respiratory tract. GPRA-A has additional expression in respiratory smooth muscle cells. However, in bronchial biopsies with unknown haplotypes, GPRA-B was upregulated in airways of all patient samples in contrast to the lack of expression in controls. Further support for GPRA as a common mediator of inflammation was obtained from a mouse model of ovalbumin-induced inflammation, where metacholine-induced airway hyperresponsiveness correlated with elevated GPRA mRNA levels in the lung and increased GPRA immunostaining in pulmonary macrophages. A novel GPRA agonist, Neuropeptide S (NPS), stimulated phagocytosis of Esterichia coli bacteria in a mouse macrophage cell line indicating a role for GPRA in the removal of inhaled allergens. The suggested GPRA functions prompted us to study, whether GPRA haplotypes associate with respiratory distress syndrome (RDS) and bronchopulmonary dysplasia (BPD) in infants sharing clinical symptoms with asthma. According to the results, near-term RDS and asthma may also share the same susceptibility and protective GPRA haplotypes. As in asthma, GPRA-B isoform expression was induced in bronchial smooth muscle cells in RDS and BPD suggesting a role for GPRA in bronchial hyperresponsiveness. In conclusion, the results of the present study suggest that the dysregulation of the GPRA/NPS pathway may not only be limited to the individuals carrying the risk variants of the gene but is also involved in the regulation of immune functions of asthma.