The linking of plate tectonics and evolutionary divergence

It is exciting to be living at a time when the big questions in biology can be investigated using modern genetics and computing [1]. Bauzà-Ribot et al.[2] take on one of the fundamental drivers of biodiversity, the effect of continental drift in the formation of the world’s biota 3 and 4, employing next-generation sequencing of whole mitochondrial genomes and modern Bayesian relaxed molecular clock analysis. Bauzà-Ribot et al.[2] conclude that vicariance via plate tectonics best explains the genetic divergence between subterranean metacrangonyctid amphipods currently found on islands separated by the Atlantic Ocean. This finding is a big deal in biogeography, and science generally [3], as many other presumed biotic tectonic divergences have been explained as probably due to more recent transoceanic dispersal events [4]. However, molecular clocks can be problematic 5 and 6 and we have identified three issues with the analyses of Bauzà-Ribot et al.[2] that cast serious doubt on their results and conclusions. When we reanalyzed their mitochondrial data and attempted to account for problems with calibration 5 and 6, modeling rates across branches 5 and 7 and substitution saturation [5], we inferred a much younger date for their key node. This implies either a later trans-Atlantic dispersal of these crustaceans, or more likely a series of later invasions of freshwaters from a common marine ancestor, but either way probably not ancient tectonic plate movements.

Time and space in biogeography : response to Parenti & Ebach (2013)

A recent Guest Editorial by Parenti & Ebach (2013, Journal of Biogeography, 40, 813–820) disagrees with the methods or interpretations in two of our recent papers. In addition, the authors open a debate on biogeographical concepts, and present an alternative philosophy for biogeographical research in the context of their recently described biogeographical subregion called ‘Pandora’. We disagree with their approach and conclusions, and comment on several issues related to our differing conceptual approaches for biogeographical research; namely, our use of molecular phylogenetic analyses, including time estimates; and Parenti & Ebach's reliance on taxon/general area cladograms. Finally, we re-examine their ‘tests’ supporting the existence of ‘Pandora’.

A phylogenetic analysis and taxonomic revision of the oribatid mite family Malaconothridae (Acari: Oribatida), with new species of Tyrphonothrus and Malaconothrus from Australia

Hitherto, the Malaconothridae contained Malaconothrus Berlese, 1904 and Trimalaconothrus Berlese, 1916, defined by the possession of one pre-tarsal claw (monodactyly) or by three claws (tridactyly) respectively. However, monodactyly is a convergent apomorphy within the Oribatida and an unreliable character for a classification. Therefore we undertook a phylogenetic analysis of 102 species as the basis for a taxonomic review of the Malaconothridae. We identified two major clades, equivalent to the genera Tyrphonothrus Knülle, 1957 and Malaconothrus. These genera are redefined. Trimala-conothrus becomes the junior subjective synonym of Malaconothrus. Some 42 species of Trimalaconothrus are recom-bined to Malaconothrus and 15 species to Tyrphonothrus. Homonyms created by the recombinations are rectified. The replacement name M. hammerae nom. nov. is proposed for M. angulatus Hammer, 1958, the junior homonym of M. an-gulatus (Willmann, 1931) and the replacement name M. luxtoni nom. nov. is proposed for M. scutatus Luxton, 1987, the junior homonym of M. scutatus Mihelč ič, 1959. Trimalaconothrus iteratus Subías, 2004 is an unnecessary replacement name and is a junior objective synonym of Malaconothrus longirostrum (Hammer 1966). Malaconothrus praeoccupatus Subías, 2004 is a junior objective synonym of M. machadoi Balogh & Mahunka, 1969. Malaconothrus obsessus (Subías, 2004), an unnecessary replacement name for Trimalaconothrus albulus Hammer 1966 sensu Tseng 1982, becomes an available name for what is in fact a previously-undescribed species of Malaconothrus. We describe four new species of Tyrphonothrus: T. gnammaensis sp. nov. from Western Australia, T. gringai sp. nov. and T. maritimus sp. nov. from New South Wales, and T. taylori sp. nov. from Queensland. We describe six new species of Malaconothrus: M. beecroftensis sp. nov., M. darwini sp. nov. M. gundungurra sp. nov. and M. knuellei sp. nov. from New South Wales, M. jowettae sp. nov. from Norfolk Island, and M. talaitae sp. nov. from Victoria.

Comparing genotyping algorithms for Illumina's Infinium whole-genome SNP BeadChips

Background Illumina's Infinium SNP BeadChips are extensively used in both small and large-scale genetic studies. A fundamental step in any analysis is the processing of raw allele A and allele B intensities from each SNP into genotype calls (AA, AB, BB). Various algorithms which make use of different statistical models are available for this task. We compare four methods (GenCall, Illuminus, GenoSNP and CRLMM) on data where the true genotypes are known in advance and data from a recently published genome-wide association study. Results In general, differences in accuracy are relatively small between the methods evaluated, although CRLMM and GenoSNP were found to consistently outperform GenCall. The performance of Illuminus is heavily dependent on sample size, with lower no call rates and improved accuracy as the number of samples available increases. For X chromosome SNPs, methods with sex-dependent models (Illuminus, CRLMM) perform better than methods which ignore gender information (GenCall, GenoSNP). We observe that CRLMM and GenoSNP are more accurate at calling SNPs with low minor allele frequency than GenCall or Illuminus. The sample quality metrics from each of the four methods were found to have a high level of agreement at flagging samples with unusual signal characteristics. Conclusions CRLMM, GenoSNP and GenCall can be applied with confidence in studies of any size, as their performance was shown to be invariant to the number of samples available. Illuminus on the other hand requires a larger number of samples to achieve comparable levels of accuracy and its use in smaller studies (50 or fewer individuals) is not recommended.

Familial recurrence risks for multiple sclerosis in Australia

BACKGROUND: Genetic susceptibility to multiple sclerosis (MS) has been recognised for many years. Considerable data exist from the northern hemisphere regarding the familial recurrence risks for MS, but there are few data for the southern hemisphere and regions at lower latitude such as Australia. To investigate the interaction between environmental and genetic causative factors in MS, the authors undertook a familial recurrence risk study in three latitudinally distinct regions of Australia. METHODS: Immediate and extended family pedigrees have been collected for three cohorts of people with MS in Queensland, Victoria and Tasmania spanning 15° of latitude. Age of onset data from Queensland were utilised to estimate age-adjusted recurrence rates. RESULTS: Recurrence risks in Australia were significantly lower than in studies from northern hemisphere populations. The age-adjusted risk for siblings across Australia was 2.13% compared with 3.5% for the northern hemisphere. A similar pattern was seen for other relatives. The risks to relatives were proportional to the population risks for each site, and hence the sibling recurrence-risk ratio (λ(s)) was similar across all sites. DISCUSSION: The familial recurrence risk of MS in Australia is lower than in previously reported studies. This is directly related to the lower population prevalence of MS. The overall genetic susceptibility in Australia as measured by the λ(s) is similar to the northern hemisphere, suggesting that the difference in population risk is explained largely by environmental factors rather than by genetic admixture.

Genome-wide meta-analysis identifies novel multiple sclerosis susceptibility loci

Objective: To perform a 1-stage meta-analysis of genome-wide association studies (GWAS) of multiple sclerosis (MS) susceptibility and to explore functional consequences of new susceptibility loci. Methods: We synthesized 7 MS GWAS. Each data set was imputed using HapMap phase II, and a per single nucleotide polymorphism (SNP) meta-analysis was performed across the 7 data sets. We explored RNA expression data using a quantitative trait analysis in peripheral blood mononuclear cells (PBMCs) of 228 subjects with demyelinating disease. Results: We meta-analyzed 2,529,394 unique SNPs in 5,545 cases and 12,153 controls. We identified 3 novel susceptibility alleles: rs170934T at 3p24.1 (odds ratio [OR], 1.17; p ¼ 1.6 � 10�8) near EOMES, rs2150702G in the second intron of MLANA on chromosome 9p24.1 (OR, 1.16; p ¼ 3.3 � 10�8), and rs6718520A in an intergenic region on chromosome 2p21, with THADA as the nearest flanking gene (OR, 1.17; p ¼ 3.4 � 10�8). The 3 new loci do not have a strong cis effect on RNA expression in PBMCs. Ten other susceptibility loci had a suggestive p < 1 � 10�6, some of these loci have evidence of association in other inflammatory diseases (ie, IL12B, TAGAP, PLEK, and ZMIZ1). Interpretation: We have performed a meta-analysis of GWAS in MS that more than doubles the size of previous gene discovery efforts and highlights 3 novel MS susceptibility loci. These and additional loci with suggestive evidence of association are excellent candidates for further investigations to refine and validate their role in the genetic architecture of MS.

Towards a translational medical research ecosystem

In this paper we introduce a novel design for a translational medical research ecosystem. Translational medical research is an emerging field of work, which aims to bridge the gap between basic medical science research and clinical research/patient care. We analyze the key challenges of digital ecosystems for translational research, based on real world scenarios posed by the Lab for Translational Research at the Harvard Medical School and the Genomics Research Centre of the Griffith University, and show how traditional IT approaches fail to fulfill these challenges. We then introduce our design for a translational research ecosystem. Several key contributions are made: A novel approach to managing ad-hoc research ecosystems is introduced; a new security approach for translational research is proposed which allows each participating site to retain control over its data and define its own policies to ensure legal and ethical compliance; and a design for a novel interactive access control framework which allows users to easily share data, while adhering to their organization's policies is presented.

A polymorphism in the HLA-DPB1 gene is associated with susceptibility to multiple sclerosis

We conducted an association study across the human leukocyte antigen (HLA) complex to identify loci associated with multiple sclerosis (MS). Comparing 1927 SNPs in 1618 MS cases and 3413 controls of European ancestry, we identified seven SNPs that were independently associated with MS conditional on the others (each ). All associations were significant in an independent replication cohort of 2212 cases and 2251 controls () and were highly significant in the combined dataset (). The associated SNPs included proxies for HLA-DRB1*15:01 and HLA-DRB1*03:01, and SNPs in moderate linkage disequilibrium (LD) with HLA-A*02:01, HLA-DRB1*04:01 and HLA-DRB1*13:03. We also found a strong association with rs9277535 in the class II gene HLA-DPB1 (discovery set , replication set , combined ). HLA-DPB1 is located centromeric of the more commonly typed class II genes HLA-DRB1, -DQA1 and -DQB1. It is separated from these genes by a recombination hotspot, and the association is not affected by conditioning on genotypes at DRB1, DQA1 and DQB1. Hence rs9277535 represents an independent MS-susceptibility locus of genome-wide significance. It is correlated with the HLA-DPB1*03:01 allele, which has been implicated previously in MS in smaller studies. Further genotyping in large datasets is required to confirm and resolve this association.

Multiple sclerosis susceptibility-associated SNPs do not influence disease severity measures in a cohort of Australian MS patients

Recent association studies in multiple sclerosis (MS) have identified and replicated several single nucleotide polymorphism (SNP) susceptibility loci including CLEC16A, IL2RA, IL7R, RPL5, CD58, CD40 and chromosome 12q13–14 in addition to the well established allele HLA-DR15. There is potential that these genetic susceptibility factors could also modulate MS disease severity, as demonstrated previously for the MS risk allele HLA-DR15. We investigated this hypothesis in a cohort of 1006 well characterised MS patients from South-Eastern Australia. We tested the MS-associated SNPs for association with five measures of disease severity incorporating disability, age of onset, cognition and brain atrophy. We observed trends towards association between the RPL5 risk SNP and time between first demyelinating event and relapse, and between the CD40 risk SNP and symbol digit test score. No associations were significant after correction for multiple testing. We found no evidence for the hypothesis that these new MS disease risk-associated SNPs influence disease severity.

Genome-wide association study identifies new multiple sclerosis susceptibility loci on chromosomes 12 and 20

To identify multiple sclerosis (MS) susceptibility loci, we conducted a genome-wide association study (GWAS) in 1,618 cases and used shared data for 3,413 controls. We performed replication in an independent set of 2,256 cases and 2,310 controls, for a total of 3,874 cases and 5,723 controls. We identified risk-associated SNPs on chromosome 12q13-14 (rs703842, P = 5.4 x 10(-11); rs10876994, P = 2.7 x 10(-10); rs12368653, P = 1.0 x 10(-7)) and upstream of CD40 on chromosome 20q13 (rs6074022, P = 1.3 x 10(-7); rs1569723, P = 2.9 x 10(-7)). Both loci are also associated with other autoimmune diseases. We also replicated several known MS associations (HLA-DR15, P = 7.0 x 10(-184); CD58, P = 9.6 x 10(-8); EVI5-RPL5, P = 2.5 x 10(-6); IL2RA, P = 7.4 x 10(-6); CLEC16A, P = 1.1 x 10(-4); IL7R, P = 1.3 x 10(-3); TYK2, P = 3.5 x 10(-3)) and observed a statistical interaction between SNPs in EVI5-RPL5 and HLA-DR15 (P = 0.001).

A pharmacogenomic evaluation of migraine therapy

Migraine is a common idiopathic primary headache disorder with significant mental, physical and social health implications. Accompanying an intense unilateral pulsating head pain other characteristic migraine symptoms include nausea, emesis, phonophobia, photophobia and in approximately 20-30% of migraine cases, neurologic disturbances associated with the aura phase. Although selective serotonin (5-HT) receptor agonists (i.e., 5-HT(1B/1D)) are successful in alleviating migrainous symptoms in < or = 70% of known sufferers, for the remaining 30%, additional migraine abortive medications remain unsuccessful, not tested or yet to be identified. Genetic characterization of the migrainous disorder is making steady progress with an increasing number of genomic susceptibility loci now identified on chromosomes 1q, 4q, 5q, 6p, 11q, 14q, 15q, 17p, 18q, 19p and Xq. The 4q, 5q, 17p and 18q loci involve endophenotypic susceptibility regions for various migrainous symptoms. In an effort to develop individualized pharmacotherapeutics, the identification of these migraine endophenotypic loci may well be the catalyst needed to aid in this goal. In this review the authors discuss the present treatment of migraine, known genomic susceptibility regions and results from migraine (genetic) association studies. The authors also discuss pharmacogenomic considerations for more individualized migraine prophylactic treatments.

A genetic analysis of serotonergic biosynthetic and metabolic enzymes in migraine using a DNA pooling approach

Migraine is a common debilitating primary headache disorder with significant mental, physical and social health implications. The brain neurotransmitter 5-hydroxytryptamine (5-HT; serotonin) is involved in nociceptive pathways and has been implicated in the pathophysiology of migraine. With few genetic studies investigating biosynthetic and metabolic enzymes governing the rate of 5-HT activity and their relationship to migraine, it was the objective of this study to assess genetic variants within the human tryptophan hydroxylase (TPH), amino acid decarboxylase (AADC) and monoamine oxidase A (MAOA) genes in migraine susceptibility. This objective was undertaken using a high-throughput DNA pooling experimental design, which proved to be a very accurate, sensitive and specific method of estimating allele frequencies for single nucleotide polymorphism, insertion deletion and variable number tandem repeat loci. Application of DNA pooling to a wide array of genetic loci provides greater scope in the assessment of population-based genetic association study designs. Despite the application of this high-throughput genotyping method, negative results from the two-stage DNA pooling design used to screen loci within the TPH, AADC and MAOA genes did not support their role in migraine susceptibility.

Association of the NuMA region on chromosome 11q13 with breast cancer susceptibility

The development of breast cancer is a complex process that involves multiple genes at many stages, from initial cell cycle dysregulation to disease progression. To identify genetic variations that influence this process, we conducted a large-scale association study using a collection of German cases and controls and >25,000 SNPs located within 16,000 genes. One of the loci identified was located on chromosome 11q13 [odds ratio (OR)=1.85, P=0.017]. The initial association was subsequently tested in two independent breast cancer collections. In both sample sets, the frequency of the susceptibility allele was increased in the cases (OR=1.6, P=0.01). The susceptibility allele was also associated with an increase in cancer family history (P=0.1). Fine mapping showed that the region of association extends approximately 300 kb and spans several genes, including the gene encoding the nuclear mitotic apparatus protein (NuMA). A nonsynonymous SNP (A794G) in NuMA was identified that showed a stronger association with breast cancer risk than the initial marker SNP (OR=2.8, P=0.005 initial sample; OR=2.1, P=0.002 combined). NuMA is a cell cycle-related protein essential for normal mitosis that is degraded in early apoptosis. NuMA-retinoic acid receptor alpha fusion proteins have been described in acute promyelocytic leukemia. Although the potential functional relevance of the A794G variation requires further biological validation, we conclude that variations in the NuMA gene are likely responsible for the observed increased breast cancer risk.

Genetic polymorphisms in DPF3 associated with risk of breast cancer and lymph node metastases

Background Several studies have identified rare genetic variations responsible for many cases of familial breast cancer but their contribution to total breast cancer incidence is relatively small. More common genetic variations with low penetrance have been postulated to account for a higher proportion of the population risk of breast cancer. Methods and Results In an effort to identify genes that influence non-familial breast cancer risk, we tested over 25,000 single nucleotide polymorphisms (SNPs) located within approximately 14,000 genes in a large-scale case-control study in 254 German women with breast cancer and 268 age-matched women without malignant disease. We identified a marker on chromosome 14q24.3-q31.1 that was marginally associated with breast cancer status (OR = 1.5, P = 0.07). Genotypes for this SNP were also significantly associated with indicators of breast cancer severity, including presence of lymph node metastases ( P = 0.006) and earlier age of onset ( P = 0.01). The association with breast cancer status was replicated in two independent samples (OR = 1.35, P = 0.05). High-density association fine mapping showed that the association spanned about 80 kb of the zinc-finger gene DPF3 (also known as CERD4 ). One SNP in intron 1 was found to be more strongly associated with breast cancer status in all three sample collections (OR = 1.6, P = 0.003) as well as with increased lymph node metastases ( P = 0.01) and tumor size ( P = 0.01). Conclusion Polymorphisms in the 5' region of DPF3 were associated with increased risk of breast cancer development, lymph node metastases, age of onset, and tumor size in women of European ancestry. This large-scale association study suggests that genetic variation in DPF3 contributes to breast cancer susceptibility and severity.

Variation in the vitamin D receptor gene is associated with multiple sclerosis in an Australian population

Multiple Sclerosis (MS) is a chronic inflammatory demyelinating disease of the central nervous system (CNS) resulting in accumulating neurological disability. The disorder is more prevalent at higher latitudes. To investigate VDR gene variation using three intragenic restriction fragment length polymorphisms (Apa I, Taq I and Fok I) in an Australian MS case-control population. One hundred and four Australian MS patients were studied with patients classified clinically as Relapsing Remitting MS (RR-MS), Secondary Progressive MS (SP-MS) or Primary Progressive MS (PP-MS). Also, 104 age-, sex-, and ethnicity-matched controls were investigated as a comparative group. Our results show a significant difference of genotype distribution frequency between the case and control groups for the functional exon 9 VDR marker Taq I (p(Gen) = 0.016) and interestingly, a stronger difference for the allelic frequency (p(All) = 0.0072). The Apa I alleles were also found to be associated with MS (p(All) = 0.04) but genotype frequencies were not significantly different from controls (p(Gen) = 0.1). The Taq and Apa variants are in very strong and significant linkage disequilibrium (D' = 0.96, P < 0.0001). The genotypic associations are strongest for the progressive forms of MS (SP-MS and PP-MS). Our results support a role for the VDR gene increasing the risk of developing multiple sclerosis, particularly the progressive clinical subtypes of MS.