Genetic variation in cytokine-related genes and migraine susceptibility

Migraine is classified by the World Health Organization (WHO) as being one of the top 20 most debilitating diseases. According to the neurovascular hypothesis, neuroinflammation may promote the activation and sensitisation of meningeal nociceptors, inducing the persistent throbbing headache characterized in migraine. The tumor necrosis factor (TNF) gene cluster, made up of TNFα, lymphotoxin α (LTA), and lymphotoxin β (LTB), has been implicated to influence the intensity and duration of local inflammation. It is thought that sterile inflammation mediated by LTA, LTB, and TNFα contributes to threshold brain excitability, propagation of neuronal hyperexcitability and thus initiation and maintenance of a migraine attack. Previous studies have investigated variants within the TNF gene cluster region in relation to migraine susceptibility, with largely conflicting results. The aim of this study was to expand on previous research and utilize a large case-control cohort and range of variants within the TNF gene cluster to investigate the role of the TNF gene cluster in migraine. Nine single nucleotide polymorphisms (SNPs) were selected for investigation as follows: rs1800683, rs2229094, rs2009658, rs2071590, rs2239704, rs909253, rs1800630, rs1800629, and rs3093664. No significant association with migraine susceptibility was found for any of the SNPs tested, with further testing according to migraine subtype and gender also showing no association for disease risk. Haplotype analysis showed that none of the tested haplotypes were significantly associated with migraine.

Impacts of climatic factors on evolution of molecular diversity and the natural distribution of wild stocks of the giant freshwater prawn (Macrobrachium rosenbergii)

Natural distributions of most freshwater taxa are restricted geographically, a pattern that reflects dispersal limitation. Macrobrachium rosenbergii is unusual because it occurs naturally in rivers from near Pakistan in the west, across India and Bangladesh to the Malay Peninsula, and across the Sunda Shelf and Indonesian archipelago to western Java. Individuals cannot tolerate full marine conditions, so dispersal between river drainage basins must occur at limited geographical scales when ecological or climatic factors are favorable. We examined molecular diversity in wild populations of M. rosenbergii across its complete natural range to document patterns of diversity and to relate them to factors that have driven evolution of diversity in this species. We found 3 clades in the mitochondrial deoxyribonucleic acid (mtDNA) data set that corresponded geographically with eastern, central, and western sets of haplotypes that last shared a common ancestor 1 × 106 y ago. The eastern clade was closest to the common ancestor of all 3 clades and to the common ancestor with its congener, Macrobrachium spinipes, distributed east of Huxley's Line. Macrobrachium rosenbergii could have evolved in the western Indonesian archipelago and spread westward during the early to mid-Pleistocene to India and Sri Lanka. Additional groups identified in the nuclear DNA data set in the central and western clades probably indicate secondary contact via dispersal between regions and modern introductions that have mixed nuclear and mtDNA genes. Pleistocene sea-level fluctuations can explain dispersal across the Indonesian archipelago and parts of mainland southeastern Asia via changing river drainage connections in shallow seas on wide continental shelves. At the western end of the modern distribution where continental shelves are smaller, intermittent freshwater plumes from large rivers probably permitted larval dispersal across inshore areas of lowered salinity.

Feral pig populations are structured at fine spatial scales in tropical Queensland, Australia

Feral pigs occur throughout tropical far north Queensland, Australia and are a significant threat to biodiversity and World Heritage values, agriculture and are a vector of infectious diseases. One of the constraints on long-lasting, local eradication of feral pigs is the process of reinvasion into recently controlled areas. This study examined the population genetic structure of feral pigs in far north Queensland to identify the extent of movement and the scale at which demographically independent management units exist. Genetic analysis of 328 feral pigs from the Innisfail to Tully region of tropical Queensland was undertaken. Seven microsatellite loci were screened and Bayesian clustering methods used to infer population clusters. Sequence variation at the mitochondrial DNA control region was examined to identify pig breed. Significant population structure was identified in the study area at a scale of 25 to 35 km, corresponding to three demographically independent management units (MUs). Distinct natural or anthropogenic barriers were not found, but environmental features such as topography and land use appear to influence patterns of gene flow. Despite the strong, overall pattern of structure, some feral pigs clearly exhibited ancestry from a MU outside of that from which they were sampled indicating isolated long distance dispersal or translocation events. Furthermore, our results suggest that gene flow is restricted among pigs of domestic Asian and European origin and non-random mating influences management unit boundaries. We conclude that the three MUs identified in this study should be considered as operational units for feral pig control in far north Queensland. Within a MU, coordinated and simultaneous control is required across farms, rainforest areas and National Park Estates to prevent recolonisation from adjacent localities.

Genome-wide SNP detection, validation, and development of an 8K SNP array for apple

As high-throughput genetic marker screening systems are essential for a range of genetics studies and plant breeding applications, the International RosBREED SNP Consortium (IRSC) has utilized the Illumina Infinium® II system to develop a medium- to high-throughput SNP screening tool for genome-wide evaluation of allelic variation in apple (Malus×domestica) breeding germplasm. For genome-wide SNP discovery, 27 apple cultivars were chosen to represent worldwide breeding germplasm and re-sequenced at low coverage with the Illumina Genome Analyzer II. Following alignment of these sequences to the whole genome sequence of 'Golden Delicious', SNPs were identified using SoapSNP. A total of 2,113,120 SNPs were detected, corresponding to one SNP to every 288 bp of the genome. The Illumina GoldenGate® assay was then used to validate a subset of 144 SNPs with a range of characteristics, using a set of 160 apple accessions. This validation assay enabled fine-tuning of the final subset of SNPs for the Illumina Infinium® II system. The set of stringent filtering criteria developed allowed choice of a set of SNPs that not only exhibited an even distribution across the apple genome and a range of minor allele frequencies to ensure utility across germplasm, but also were located in putative exonic regions to maximize genotyping success rate. A total of 7867 apple SNPs was established for the IRSC apple 8K SNP array v1, of which 5554 were polymorphic after evaluation in segregating families and a germplasm collection. This publicly available genomics resource will provide an unprecedented resolution of SNP haplotypes, which will enable marker-locus-trait association discovery, description of the genetic architecture of quantitative traits, investigation of genetic variation (neutral and functional), and genomic selection in apple.

Genetic bottlenecks in time and space: reconstructing invasions from contemporary and historical collections

Herbarium accession data offer a useful historical botanical perspective and have been used to track the spread of plant invasions through time and space. Nevertheless, few studies have utilised this resource for genetic analysis to reconstruct a more complete picture of historical invasion dynamics, including the occurrence of separate introduction events. In this study, we combined nuclear and chloroplast microsatellite analyses of contemporary and historical collections of Senecio madagascariensis, a globally invasive weed first introduced to Australia c. 1918 from its native South Africa. Analysis of nuclear microsatellites, together with temporal spread data and simulations of herbarium voucher sampling, revealed distinct introductions to south-eastern Australia and mid-eastern Australia. Genetic diversity of the south-eastern invasive population was lower than in the native range, but higher than in the mid-eastern invasion. In the invasive range, despite its low resolution, our chloroplast microsatellite data revealed the occurrence of new haplotypes over time, probably as the result of subsequent introduction(s) to Australia from the native range during the latter half of the 20th century. Our work demonstrates how molecular studies of contemporary and historical field collections can be combined to reconstruct a more complete picture of the invasion history of introduced taxa. Further, our study indicates that a survey of contemporary samples only (as undertaken for the majority of invasive species studies) would be insufficient to identify potential source populations and occurrence of multiple introductions.

Evidence of associations between cytokine gene polymorphisms and quality of life in patients with cancer and their family caregivers

PURPOSE/OBJECTIVES: To identify latent classes of individuals with distinct quality-of-life (QOL) trajectories, to evaluate for differences in demographic characteristics between the latent classes, and to evaluate for variations in pro- and anti-inflammatory cytokine genes between the latent classes. DESIGN: Descriptive, longitudinal study. SETTING: Two radiation therapy departments located in a comprehensive cancer center and a community-based oncology program in northern California. SAMPLE: 168 outpatients with prostate, breast, brain, or lung cancer and 85 of their family caregivers (FCs). METHODS: Growth mixture modeling (GMM) was employed to identify latent classes of individuals based on QOL scores measured prior to, during, and for four months following completion of radiation therapy. Single nucleotide polymorphisms (SNPs) and haplotypes in 16 candidate cytokine genes were tested between the latent classes. Logistic regression was used to evaluate the relationships among genotypic and phenotypic characteristics and QOL GMM group membership. MAIN RESEARCH VARIABLES: QOL latent class membership and variations in cytokine genes. FINDINGS: Two latent QOL classes were found: higher and lower. Patients and FCs who were younger, identified with an ethnic minority group, had poorer functional status, or had children living at home were more likely to belong to the lower QOL class. After controlling for significant covariates, between-group differences were found in SNPs in interleukin 1 receptor 2 (IL1R2) and nuclear factor kappa beta 2 (NFKB2). For IL1R2, carrying one or two doses of the rare C allele was associated with decreased odds of belonging to the lower QOL class. For NFKB2, carriers with two doses of the rare G allele were more likely to belong to the lower QOL class. CONCLUSIONS: Unique genetic markers in cytokine genes may partially explain interindividual variability in QOL. IMPLICATIONS FOR NURSING: Determination of high-risk characteristics and unique genetic markers would allow for earlier identification of patients with cancer and FCs at higher risk for poorer QOL. Knowledge of these risk factors could assist in the development of more targeted clinical or supportive care interventions for those identified.

One and the same: Integrative taxonomic evidence that Bactrocera invadens (Diptera: Tephritidae) is the same species as the Oriental fruit fly Bactrocera dorsalis

The invasive fruit fly Bactrocera invadens Drew, Tsuruta & White, and the Oriental fruit fly Bactrocera dorsalis (Hendel) are highly destructive horticultural pests of global significance. Bactrocera invadens originates from the Indian subcontinent and has recently invaded all of sub-Saharan Africa, while B. dorsalis principally occurs from the Indian subcontinent towards southern China and South-east Asia. High morphological and genetic similarity has cast doubt over whether B. invadens is a distinct species from B. dorsalis. Addressing this issue within an integrative taxonomic framework, we sampled from across the geographic distribution of both taxa and: (i) analysed morphological variation, including those characters considered diagnostic (scutum colour, length of aedeagus, width of postsutural lateral vittae, wing size, and wing shape); (ii) sequenced four loci (ITS1, ITS2, cox1 and nad4) for phylogenetic inference, and; (iii) generated a cox1 haplotype network to examine population structure. Molecular analyses included the closely related species, Bactrocera kandiensis Drew & Hancock. Scutum colour varies from red-brown to fully black for individuals from Africa and the Indian subcontinent. All individuals east of the Indian subcontinent are black except for a few red-brown individuals from China. The postsutural lateral vittae width of B. invadens is narrower than B. dorsalis from eastern Asia, but the variation is clinal, with subcontinent B. dorsalis populations intermediate in size. Aedeagus length, wing shape and wing size cannot discriminate between the two taxa. Phylogenetic analyses failed to resolve B. invadens from B. dorsalis, but did resolve B. kandiensis. Bactrocera dorsalis and B. invadens shared cox1 haplotypes, yet the haplotype network pattern does not reflect current taxonomy or patterns in thoracic colour. Some individuals of B. dorsalis/B. invadens possessed haplotypes more closely related to B. kandiensis than to conspecifics, suggestive of mitochondrial introgression between these species. The combined evidence fails to support the delimitation of B. dorsalis and B. invadens as separate biological species. Consequently, existing biological data for B. dorsalis may be applied to the invasive population in Africa. Our recommendation, in line with other recent publications, is that B. invadens be synonymized with B. dorsalis.

Genetic divergence and echolocation call frequency in cryptic species of Hipposideros larvatus s.l. (Chiroptera : Hipposideridae) from the Indo-Malayan region

The intermediate leaf-nosed bat (Hipposideros larvatus) is a medium-sized bat distributed throughout the Indo-Malay region. In north-east India, bats identified as H. larvatus captured at a single cave emitted echolocation calls with a bimodal distribution of peak frequencies, around either 85 kHz or 98 kHz. Individuals echolocating at 85 kHz had larger ears and longer forearms than those echolocating at 98 kHz, although no differences were detected in either wing morphology or diet, suggesting limited resource partitioning. A comparison of mitochondrial control region haplotypes of the two phonic types with individuals sampled from across the Indo-Malay range supports the hypothesis that, in India, two cryptic species are present. The Indian 98-kHz phonic bats formed a monophyletic clade with bats from all other regional populations sampled, to the exclusion of the Indian 85-kHz bats. In India, the two forms showed 12–13% sequence divergence and we propose that the name Hipposideros khasiana for bats of the 85-kHz phonic type. Bats of the 98-kHz phonic type formed a monophyletic group with bats from Myanmar, and corresponded to Hipposideros grandis, which is suggested to be a species distinct from Hipposideros larvatus. Differences in echolocation call frequency among populations did not reflect phylogenetic relationships, indicating that call frequency is a poor indicator of evolutionary history. Instead, divergence in call frequency probably occurs in allopatry, possibly augmented by character displacement on secondary contact to facilitate intraspecific communication.

Atlantic-Mediterranean and within-Mediterranean molecular variation in Coris julis (L. 1758) (Teleostei, Labridae)

Sequence variation in the mitochondrial control region was studied in the Mediterranean rainbow wrasse (Coris julis), a species with pronounced pelagic larval phase inhabiting the Mediterranean Sea and the adjacent coastal eastern Atlantic Ocean. A total of 309 specimens from 19 sampling sites were analysed with the aim of elucidating patterns of molecular variation between the Atlantic and the Mediterranean as well as within the Mediterranean Sea. Phylogeographic analyses revealed a pronounced structuring into a Mediterranean and an Atlantic group. Samples from a site at the Moroccan Mediterranean coast in the Alboran Sea showed intermediate frequencies of “Mediterranean” and “Atlantic” haplotypes. We recognised a departure from molecular neutrality and a star-like genealogy for samples from the Mediterranean Sea, which we propose to have happened due to a recent demographic expansion. The results are discussed in the light of previous studies on molecular variation in fish species between the Atlantic and the Mediterranean and within the Mediterranean.

Identity-by-Descent Mapping to Detect Rare Variants Conferring Susceptibility to Multiple Sclerosis

Genome-wide association studies (GWAS) have identified around 60 common variants associated with multiple sclerosis (MS), but these loci only explain a fraction of the heritability of MS. Some missing heritability may be caused by rare variants that have been suggested to play an important role in the aetiology of complex diseases such as MS. However current genetic and statistical methods for detecting rare variants are expensive and time consuming. 'Population-based linkage analysis' (PBLA) or so called identity-by-descent (IBD) mapping is a novel way to detect rare variants in extant GWAS datasets. We employed BEAGLE fastIBD to search for rare MS variants utilising IBD mapping in a large GWAS dataset of 3,543 cases and 5,898 controls. We identified a genome-wide significant linkage signal on chromosome 19 (LOD = 4.65; p = 1.9×10-6). Network analysis of cases and controls sharing haplotypes on chromosome 19 further strengthened the association as there are more large networks of cases sharing haplotypes than controls. This linkage region includes a cluster of zinc finger genes of unknown function. Analysis of genome wide transcriptome data suggests that genes in this zinc finger cluster may be involved in very early developmental regulation of the CNS. Our study also indicates that BEAGLE fastIBD allowed identification of rare variants in large unrelated population with moderate computational intensity. Even with the development of whole-genome sequencing, IBD mapping still may be a promising way to narrow down the region of interest for sequencing priority. © 2013 Lin et al.

The interleukin 1 gene cluster contains a major susceptibility locus for ankylosing spondylitis

Ankylosing spondylitis (AS) is a common and highly heritable inflammatory arthropathy. Although the gene HLA-B27 is almost essential for the inheritance of the condition, it alone is not sufficient to explain the pattern of familial recurrence of the disease. We have previously demonstrated suggestive linkage of AS to chromosome 2q13, a region containing the interleukin 1 (IL-1) family gene cluster, which includes several strong candidates for involvement in the disease. In the current study, we describe strong association and transmission of IL-1 family gene cluster single-nucleotide polymorphisms and haplotypes with AS.

Matters Arising: Role of HLA genes in familial spondyloarthropathy

Said-Nahal and colleagues report an intriguing finding of an association with HLA-DR4 independent of B27 in families with ankylosing spondylitis (AS),1 a finding highlighted by an accompanying editorial.2 The approach of studying B27 positive and B27 negative haplotypes may prove powerful in identifying further cis or trans encoded genes involved in AS. However, the reported association of DR4 with AS is quite a surprising finding given that no difference was noted in B27-DR4 haplotype frequencies in patients and ethnically matched healthy controls. Many previous studies have not reported any such association,3–9 including a similar preliminary study by the same authors.10 Although these studies were mainly case-control studies, population stratification is highly unlikely to cause a false negative finding if the effect size of the reported association with DR4 is as high as Said-Nahal and colleagues describe...

Replication of association of IL1 gene complex members with ankylosing spondylitis in Taiwanese Chinese

Objective: To test the association of interleukin 1 (IL1) gene family members with ankylosing spondylitis (AS), previously reported in Europid subjects, in an ethnically remote population. Methods: 200 Taiwanese Chinese AS patients and 200 ethnically matched healthy controls were genotyped for five single nucleotide polymorphisms (SNPs) and the IL1RN.VNTR, markers previously associated with AS. Allele, genotype, and haplotype frequencies were compared between cases and controls. Results: Association of alleles and genotypes of the markers IL1F10.3, IL1RN.4, and IL1RN.VNTR was observed with AS (p<0.05). Haplotypes of pairs of these markers and of the markers IL1RN.6/1 and IL1RN.6/2 were also significantly associated with AS. The strongest associations observed were with the marker IL1RN.4, and with the two-marker haplotype IL1RN.4-IL1RN.VNTR (both p = 0.004). Strong linkage disequilibrium was observed between all marker pairs except those involving IL1B-511 (D′ 0.4 to 0.9, p<0.01). Conclusions: The IL1 gene cluster is associated with AS in Taiwanese Chinese. This finding provides strong statistical support that the previously observed association of this gene cluster with AS is a true positive finding.

Association of variants in MMEL1 and CTLA4 with rheumatoid arthritis in the Han Chinese population

Background: The genome-wide association study era has made great progress in identifying susceptibility genes and genetic loci for rheumatoid arthritis (RA) in populations of White European ancestry. However, few studies have tried to dissect disease aetiopathogenesis in other ethnic populations. Objective: To investigate these associations in the Han Chinese population. Methods: Haplotypes from the HapMap database Chinese population were used to select tag-single-nucleotide polymorphisms (SNPs) (r2 =0.8) across 19 distinct RA genomic regions. A two phase case-control association study was performed, with 169 SNPs genotyped in phase I (n=571 cases, n=880 controls), and 64 SNPs achieving p<0.2 in the first phase being genotyped in phase II (n=464 cases, n=822 controls). Association statistics were calculated using permutation tests both unadjusted and adjusted for the number of markers studied. Results: Robust association was detected for MMEL1 and CTLA4 , and modest association was identified for another six loci: PADI4 , STAT4 , PRDM1 , CDK6 , TRAF1-C5 and KIF5A-PIP4K2C. All three markers genotyped in MMEL1 demonstrated association, with peak signal for rs3890745 (p=2.6×10 -5unadjusted, p=0.003 adjusted, OR=0.79). For CTLA4 , significance was detected for three of five variants showing association, with peak association for marker rs12992492 (p=4.3×10-5 unadjusted, p=0.0021 adjusted, OR=0.77). Lack of association of common variants in PTPN22 with RA in Han Chinese was confirmed. Conclusion: This study identifies MMEL1 and CTLA4 as RA susceptibility genes, provides suggestive evidence of association for a further six loci in the Han Chinese population and confirms lack of PTPN22 association in Asian populations. It also confirms the value of multiethnic population studies to help dissect disease aetiopathogenesis.

Analysis of killer immunoglobulin-like receptor genes in ankylosing spondylitis

Objectives: To assess the possible association of killer immunoglobulin-like receptor (KIR) genes, specifically KIR3DL1, KIR3DS1 and KIR3DL2, with ankylosing spondylitis (AS). Methods: 14 KIR genes were genotyped in 200 UK patients with AS and 405 healthy controls using multiplex polymerase chain reaction. Sequence-specific oligonucleotide probes were used to subtype 368 cases with AS and 366 controls for 12 KIR3DL2 alleles. Differences in KIR genotypes and KIR3DL2 allele frequencies were assessed using the χp2p test. Results: KIR3DL1 and KIR3DS1 gene frequencies were very similar in cases with AS and controls (odds ratio = 1.5, 95% confidence interval 0.8 to 3.0, and odds ratio = 1.02, 95% confidence interval 0.2 to 5.3, respectively). KIR3DL2 allele frequencies were not significantly different between cases with AS and controls. Conclusions: Neither the KIR gene content of particular KIR haplotypes nor KIR3DL2 polymorphisms contribute to AS.