Combined mitochondrial and nuclear DNA sequences resolve the interrelations of the major Australasian marsupial radiations

Australasian marsupials include three major radiations, the insectivorous/carnivorous Dasyuromorphia, the omnivorous bandicoots (Peramelemorphia), and the largely herbivorous diprotodontians. Morphologists have generally considered the bandicoots and diprotodontians to be closely related, most prominently because they are both syndactylous (with the 2nd and 3rd pedal digits being fused). Molecular studies have been unable to confirm or reject this Syndactyla hypothesis. Here we present new mitochondrial (mt) genomes from a spiny bandicoot (Echymipera rufescens) and two dasyurids, a fat-tailed dunnart (Sminthopsis crassicaudata) and a northern quoll (Dasyurus hallucatus). By comparing trees derived from pairwise base-frequency differences between taxa with standard (absolute, uncorrected) distance trees, we infer that composition bias among mt protein-coding and RNA sequences is sufficient to mislead tree reconstruction. This can explain incongruence between trees obtained from mt and nuclear data sets. However, after excluding major sources of compositional heterogeneity, both the “reduced-bias” mt and nuclear data sets clearly favor a bandicoot plus dasyuromorphian association, as well as a grouping of kangaroos and possums (Phalangeriformes) among diprotodontians. Notably, alternatives to these groupings could only be confidently rejected by combining the mt and nuclear data. Elsewhere on the tree, Dromiciops appears to be sister to the monophyletic Australasian marsupials, whereas the placement of the marsupial mole (Notoryctes) remains problematic. More generally, we contend that it is desirable to combine mt genome and nuclear sequences for inferring vertebrate phylogeny, but as separately modeled process partitions. This strategy depends on detecting and excluding (or accounting for) major sources of nonhistorical signal, such as from compositional nonstationarity.

Family-level relationships among the Australasian marsupial “herbivores” (Diprotodontia: Koala, wombats, kangaroos and possums)

The marsupial order Diprotodontia includes 10 extant families, which occupy all terrestrial habitats across Australia and New Guinea and have evolved remarkable dietary and locomotory diversity. Despite considerable attention, the interrelations of these families have for the most part remained elusive. In this study, we separately model mitochondrial RNA and protein-coding sequences in addition to nuclear protein-coding sequences to provide near-complete resolution of diprotodontian family-level phylogeny. We show that alternative topologies inferred in some previous studies are likely to be artifactual, resulting from branch-length and compositional biases. Subordinal groupings resolved herein include Vombatiformes (wombats and koala) and Phalangerida, which in turn comprises Petauroidea (petaurid gliders and striped, feathertail, ringtail and honey possums) and a clade whose plesiomorphic members possess blade-like premolars (phalangerid possums, kangaroos and their allies and most likely, pygmy possums). The topology resolved reveals ecological niche structuring among diprotodontians that has likely been maintained for more than 40 million years.

Tinamous and moa flock together : mitochondrial genome sequence analysis reveals independent losses of flight among ratites

Ratites are large, flightless birds and include the ostrich, rheas, kiwi, emu, and cassowaries, along with extinct members, such as moa and elephant birds. Previous phylogenetic analyses of complete mitochondrial genome sequences have reinforced the traditional belief that ratites are monophyletic and tinamous are their sister group. However, in these studies ratite monophyly was enforced in the analyses that modeled rate heterogeneity among variable sites. Relaxing this topological constraint results in strong support for the tinamous (which fly) nesting within ratites. Furthermore, upon reducing base compositional bias and partitioning models of sequence evolution among protein codon positions and RNA structures, the tinamou–moa clade grouped with kiwi, emu, and cassowaries to the exclusion of the successively more divergent rheas and ostrich. These relationships are consistent with recent results from a large nuclear data set, whereas our strongly supported finding of a tinamou–moa grouping further resolves palaeognath phylogeny. We infer flight to have been lost among ratites multiple times in temporally close association with the Cretaceous–Tertiary extinction event. This circumvents requirements for transient microcontinents and island chains to explain discordance between ratite phylogeny and patterns of continental breakup. Ostriches may have dispersed to Africa from Eurasia, putting in question the status of ratites as an iconic Gondwanan relict taxon. [Base composition; flightless; Gondwana; mitochondrial genome; Palaeognathae; phylogeny; ratites.]

Ureaplasma parvum : understanding the complexities of intra-amniotic infection in an ovine model

The human Ureaplasma species are the most frequently isolated bacteria from the upper genital tract of pregnant women and can cause clinically asymptomatic, intra-uterine infections, which are difficult to treat with antimicrobials. Ureaplasma infection of the upper genital tract during pregnancy has been associated with numerous adverse outcomes including preterm birth, chorioamnionitis and neonatal respiratory diseases. The mechanisms by which ureaplasmas are able to chronically colonise the amniotic fluid and avoid eradication by (i) the host immune response and (ii) maternally-administered antimicrobials, remain virtually unexplored. To address this gap within the literature, this study investigated potential mechanisms by which ureaplasmas are able to cause chronic, intra-amniotic infections in an established ovine model. In this PhD program of research the effectiveness of standard, maternal erythromycin for the treatment of chronic, intra-amniotic ureaplasma infections was evaluated. At 55 days of gestation pregnant ewes received an intra-amniotic injection of either: a clinical Ureaplasma parvum serovar 3 isolate that was sensitive to macrolide antibiotics (n = 16); or 10B medium (n = 16). At 100 days of gestation, ewes were then randomised to receive either maternal erythromycin treatment (30 mg/kg/day for four days) or no treatment. Ureaplasmas were isolated from amniotic fluid, chorioamnion, umbilical cord and fetal lung specimens, which were collected at the time of preterm delivery of the fetus (125 days of gestation). Surprisingly, the numbers of ureaplasmas colonising the amniotic fluid and fetal tissues were not different between experimentally-infected animals that received erythromycin treatment or infected animals that did not receive treatment (p > 0.05), nor were there any differences in fetal inflammation and histological chorioamnionitis between these groups (p > 0.05). These data demonstrate the inability of maternal erythromycin to eradicate intra-uterine ureaplasma infections. Erythromycin was detected in the amniotic fluid of animals that received antimicrobial treatment (but not in those that did not receive treatment) by liquid chromatography-mass spectrometry; however, the concentrations were below therapeutic levels (<10 – 76 ng/mL). These findings indicate that the ineffectiveness of standard, maternal erythromycin treatment of intra-amniotic ureaplasma infections may be due to the poor placental transfer of this drug. Subsequently, the phenotypic and genotypic characteristics of ureaplasmas isolated from the amniotic fluid and chorioamnion of pregnant sheep after chronic, intra-amniotic infection and low-level exposure to erythromycin were investigated. At 55 days of gestation twelve pregnant ewes received an intra-amniotic injection of a clinical U. parvum serovar 3 isolate, which was sensitive to macrolide antibiotics. At 100 days of gestation, ewes received standard maternal erythromycin treatment (30 mg/kg/day for four days, n = 6) or saline (n = 6). Preterm fetuses were surgically delivered at 125 days of gestation and ureaplasmas were cultured from the amniotic fluid and the chorioamnion. The minimum inhibitory concentrations (MICs) of erythromycin, azithromycin and roxithromycin were determined for cultured ureaplasma isolates, and antimicrobial susceptibilities were different between ureaplasmas isolated from the amniotic fluid (MIC range = 0.08 – 1.0 mg/L) and chorioamnion (MIC range = 0.06 – 5.33 mg/L). However, the increased resistance to macrolide antibiotics observed in chorioamnion ureaplasma isolates occurred independently of exposure to erythromycin in vivo. Remarkably, domain V of the 23S ribosomal RNA gene (which is the target site of macrolide antimicrobials) of chorioamnion ureaplasmas demonstrated significant variability (125 polymorphisms out of 422 sequenced nucleotides, 29.6%) when compared to the amniotic fluid ureaplasma isolates and the inoculum strain. This sequence variability did not occur as a consequence of exposure to erythromycin, as the nucleotide substitutions were identical between chorioamnion ureaplasmas isolated from different animals, including those that did not receive erythromycin treatment. We propose that these mosaic-like 23S ribosomal RNA gene sequences may represent gene fragments transferred via horizontal gene transfer. The significant differences observed in (i) susceptibility to macrolide antimicrobials and (ii) 23S ribosomal RNA sequences of ureaplasmas isolated from the amniotic fluid and chorioamnion suggests that the anatomical site from which they were isolated may exert selective pressures that alter the socio-microbiological structure of the bacterial population, by selecting for genetic changes and altered antimicrobial susceptibility profiles. The final experiment for this PhD examined antigenic size variation of the multiple banded antigen (MBA, a surface-exposed lipoprotein and predicted ureaplasmal virulence factor) in chronic, intra-amniotic ureaplasma infections. Previously defined ‘virulent-derived’ and ‘avirulent-derived’ clonal U. parvum serovar 6 isolates (each expressing a single MBA protein) were injected into the amniotic fluid of pregnant ewes (n = 20) at 55 days of gestation, and amniotic fluid was collected by amniocentesis every two weeks until the time of near-term delivery of the fetus (at 140 days of gestation). Both the avirulent and virulent clonal ureaplasma strains generated MBA size variants (ranging in size from 32 – 170 kDa) within the amniotic fluid of pregnant ewes. The mean number of MBA size variants produced within the amniotic fluid was not different between the virulent (mean = 4.2 MBA variants) and avirulent (mean = 4.6 MBA variants) ureaplasma strains (p = 0.87). Intra-amniotic infection with the virulent strain was significantly associated with the presence of meconium-stained amniotic fluid (p = 0.01), which is an indicator of fetal distress in utero. However, the severity of histological chorioamnionitis was not different between the avirulent and virulent groups. We demonstrated that ureaplasmas were able to persist within the amniotic fluid of pregnant sheep for 85 days, despite the host mounting an innate and adaptive immune response. Pro-inflammatory cytokines (interleukin (IL)-1â, IL-6 and IL-8) were elevated within the chorioamnion tissue of pregnant sheep from both the avirulent and virulent treatment groups, and this was significantly associated with the production of anti-ureaplasma IgG antibodies within maternal sera (p < 0.05). These findings suggested that the inability of the host immune response to eradicate ureaplasmas from the amniotic cavity may be due to continual size variation of MBA surface-exposed epitopes. Taken together, these data confirm that ureaplasmas are able to cause long-term in utero infections in a sheep model, despite standard antimicrobial treatment and the development of a host immune response. The overall findings of this PhD project suggest that ureaplasmas are able to cause chronic, intra-amniotic infections due to (i) the limited placental transfer of erythromycin, which prevents the accumulation of therapeutic concentrations within the amniotic fluid; (ii) the ability of ureaplasmas to undergo rapid selection and genetic variation in vivo, resulting in ureaplasma isolates with variable MICs to macrolide antimicrobials colonising the amniotic fluid and chorioamnion; and (iii) antigenic size variation of the MBA, which may prevent eradication of ureaplasmas by the host immune response and account for differences in neonatal outcomes. The outcomes of this program of study have improved our understanding of the biology and pathogenesis of this highly adapted microorganism.

Identification of a novel prostate cancer susceptibility variant in the KLK3 gene transcript

Abstract Genome-wide association studies (GWAS) have identified more than 30 prostate cancer (PrCa) susceptibility loci. One of these (rs2735839) is located close to a plausible candidate susceptibility gene, KLK3, which encodes prostate-specific antigen (PSA). PSA is widely used as a biomarker for PrCa detection and disease monitoring. To refine the association between PrCa and variants in this region, we used genotyping data from a two-stage GWAS using samples from the UK and Australia, and the Cancer Genetic Markers of Susceptibility (CGEMS) study. Genotypes were imputed for 197 and 312 single nucleotide polymorphisms (SNPs) from HapMap2 and the 1000 Genome Project, respectively. The most significant association with PrCa was with a previously unidentified SNP, rs17632542 (combined P = 3.9 × 10−22). This association was confirmed by direct genotyping in three stages of the UK/Australian GWAS, involving 10,405 cases and 10,681 controls (combined P = 1.9 × 10−34). rs17632542 is also shown to be associated with PSA levels and it is a non-synonymous coding SNP (Ile179Thr) in KLK3. Using molecular dynamic simulation, we showed evidence that this variant has the potential to introduce alterations in the protein or affect RNA splicing. We propose that rs17632542 may directly influence PrCa risk.

A two-metal ion mechanism operates in the hammerhead ribozyme-mediated cleavage of an RNA substrate

Evidence for a two-metal ion mechanism for cleavage of the HH16 hammerhead ribozyme is provided by monitoring the rate of cleavage of the RNA substrate as a function of La3+ concentration in the presence of a constant concentration of Mg2+. We show that a bell-shaped curve of cleavage activation is obtained as La3+ is added in micromolar concentrations in the presence of 8 mM Mg2+, with a maximal rate of cleavage being attained in the presence of 3 microM La3+. These results show that two-metal ion binding sites on the ribozyme regulate the rate of the cleavage reaction and, on the basis of earlier estimates of the Kd values for Mg2+ of 3.5 mM and > 50 mM, that these sites bind La3+ with estimated Kd values of 0.9 and > 37.5 microM, respectively. Furthermore, given the very different effects of these metal ions at the two binding sites, with displacement of Mg2+ by La3+ at the stronger (relative to Mg2+) binding site activating catalysis and displacement of Mg2+ by La3+ at the weaker (relative to Mg2+) (relative to Mg2+) binding site inhibiting catalysis, we show that the metal ions at these two sites play very different roles. We argue that the metal ion at binding site 1 coordinates the attacking 2'-oxygen species in the reaction and lowers the pKa of the attached proton, thereby increasing the concentration of the attacking alkoxide nucleophile in an equilibrium process. In contrast, the role of the metal ion at binding site 2 is to catalyze the reaction by absorbing the negative charge that accumulates at the leaving 5'-oxygen in the transition state. We suggest structural reasons why the Mg(2+)-La3+ ion combination is particularly suited to demonstrating these different roles of the two-metal ions in the ribozyme cleavage reaction.

The adenine-rich tract in the 5ʹ-end of the hepatitis C virus ORF encodes a peptide regulating the binding of the C protein to RNA

Hepatitis C virus (HCV ) core (C) protein is thought to bind to viral RNA before it undergoes oligomerization leading to RNA encapsidation. Details of these events are so far unknown. The 5ʹ-terminal C protein coding sequence that includes an adenine (A)-rich tract is a part of an internal ribosome entry site(IRES). This nucleotide sequence but not the corresponding protein sequence is needed for proper initiation of translation of viral RNA by an IRES-dependent mechanism. In this study, we examined the importance of this sequence for the ability of the C protein to bind to viral RNA. Serially truncated C proteins with deletions from 10 up to 45 N-terminal amino acids were expressed in Escherichia coli, purified and tested for binding to viral RNA by a gel shift assay. The results showed that truncation of the C protein from its N-terminus by more than 10 amino acids abolished almost completely its expression in E. coli. The latter could be restored by adding a tag to the N-terminus of the protein. The tagged proteins truncated by 15 or more amino acids showed an anomalous migration in SDS-PAGE. Truncation by more than 20 amino acids resulted in a complete loss of ability of tagged C protein to bind to viral RNA. These results provide clues to the early events in the C protein - RNA interactions leading to C protein oligomerization, RNA encapsidation and virion assembly.

Beyond barcoding : a mitochondrial genomics approach to molecular phylogenetics and diagnostics of blowflies (Diptera: Calliphoridae)

Members of the Calliphoridae (blowflies) are significant for medical and veterinary management, due to the ability of some species to consume living flesh as larvae, and for forensic investigations due to the ability of others to develop in corpses. Due to the difficulty of accurately identifying larval blowflies to species there is a need for DNA-based diagnostics for this family, however the widely used DNA-barcoding marker, cox1, has been shown to fail for several groups within this family. Additionally, many phylogenetic relationships within the Calliphoridae are still unresolved, particularly deeper level relationships. Sequencing whole mt genomes has been demonstrated both as an effective method for identifying the most informative diagnostic markers and for resolving phylogenetic relationships. Twenty-seven complete, or nearly so, mt genomes were sequenced representing 13 species, seven genera and four calliphorid subfamilies and a member of the related family Tachinidae. PCR and sequencing primers developed for sequencing one calliphorid species could be reused to sequence related species within the same superfamily with success rates ranging from 61% to 100%, demonstrating the speed and efficiency with which an mt genome dataset can be assembled. Comparison of molecular divergences for each of the 13 protein-coding genes and 2 ribosomal RNA genes, at a range of taxonomic scales identified novel targets for developing as diagnostic markers which were 117–200% more variable than the markers which have been used previously in calliphorids. Phylogenetic analysis of whole mt genome sequences resulted in much stronger support for family and subfamily-level relationships. The Calliphoridae are polyphyletic, with the Polleninae more closely related to the Tachinidae, and the Sarcophagidae are the sister group of the remaining calliphorids. Within the Calliphoridae, there was strong support for the monophyly of the Chrysomyinae and Luciliinae and for the sister-grouping of Luciliinae with Calliphorinae. Relationships within Chrysomya were not well resolved. Whole mt genome data, supported the previously demonstrated paraphyly of Lucilia cuprina with respect to L. sericata and allowed us to conclude that it is due to hybrid introgression prior to the last common ancestor of modern sericata populations, rather than due to recent hybridisation, nuclear pseudogenes or incomplete lineage sorting.

Abrogation of contaminating RNA activity in HIV-1 Gag VLPs

Background: HIV-1 Gag virus like particles (VLPs) used as candidate vaccines are regarded as inert particles as they contain no replicative nucleic acid, although they do encapsidate cellular RNAs. During HIV-1 Gag VLP production in baculovirus-based expression systems, VLPs incorporate the baculovirus Gp64 envelope glycoprotein, which facilitates their entry into mammalian cells. This suggests that HIV-1 Gag VLPs produced using this system facilitate uptake and subsequent expression of encapsidated RNA in mammalian cells - an unfavourable characteristic for a vaccine. Methods. HIV-1 Gag VLPs encapsidating reporter chloramphenicol acetyl transferase (CAT) RNA, were made in insect cells using the baculovirus expression system. The presence of Gp64 on the VLPs was verified by western blotting and RT-PCR used to detect and quantitate encapsidated CAT RNA. VLP samples were heated to inactivate CAT RNA. Unheated and heated VLPs incubated with selected mammalian cell lines and cell lysates tested for the presence of CAT protein by ELISA. Mice were inoculated with heated and unheated VLPs using a DNA prime VLP boost regimen. Results: HIV-1 Gag VLPs produced had significantly high levels of Gp64 (∼1650 Gp64 molecules/VLP) on their surfaces. The amount of encapsidated CAT RNA/g Gag VLPs ranged between 0.1 to 7 ng. CAT protein was detected in 3 of the 4 mammalian cell lines incubated with VLPs. Incubation with heated VLPs resulted in BHK-21 and HeLa cell lysates showing reduced CAT protein levels compared with unheated VLPs and HEK-293 cells. Mice inoculated with a DNA prime VLP boost regimen developed Gag CD8 and CD4 T cell responses to GagCAT VLPs which also boosted a primary DNA response. Heating VLPs did not abrogate these immune responses but enhanced the Gag CD4 T cell responses by two-fold. Conclusions: Baculovirus-produced HIV-1 Gag VLPs encapsidating CAT RNA were taken up by selected mammalian cell lines. The presence of CAT protein indicates that encapsidated RNA was expressed in the mammalian cells. Heat-treatment of the VLPs altered the ability of protein to be expressed in some cell lines tested but did not affect the ability of the VLPs to stimulate an immune response when inoculated into mice. © 2011 Valley-Omar et al; licensee BioMed Central Ltd.

INT6/EIF3E interacts with ATM and is required for proper execution of the DNA damage response in human cells

Altered expression of the INT6 gene, encoding the e subunit of the translational initiation factor eIF3, occurs in human breast cancers, but how INT6 relates to carcinogenesis remains unestablished. Here, we show that INT6 is involved in the DNA damage response. INT6 was required for cell survival following γ-irradiation and G(2)-M checkpoint control. RNA interference-mediated silencing of INT6 reduced phosphorylation of the checkpoint kinases CHK1 and CHK2 after DNA damage. In addition, INT6 silencing prevented sustained accumulation of ataxia telangiectasia mutated (ATM) at DNA damage sites in cells treated with γ-radiation or the radiomimetic drug neocarzinostatin. Mechanistically, this result could be explained by interaction of INT6 with ATM, which together with INT6 was recruited to the sites of DNA damage. Finally, INT6 silencing also reduced ubiquitylation events that promote retention of repair proteins at DNA lesions. Accordingly, accumulation of the repair factor BRCA1 was defective in the absence of INT6. Our findings reveal unexpected and striking connections of INT6 with ATM and BRCA1 and suggest that the protective action of INT6 in the onset of breast cancers relies on its involvement in the DNA damage response.

The complete mitochondrial genome of the flesh fly, Sarcophaga impatiens Walker (Diptera: Sarcophagidae)

Approximately 2500 fly species comprise the Sarcophagidae family worldwide. The complete mitochondrial genome of the carrion-breeding, forensically important Sarcophaga impatiens Walker (Diptera: Sarcophagidae) from Australia was sequenced. The 15,169 bp circular genome contains the 37 genes found in a typical Metazoan genome: 13 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes. It also contains one non-coding A+T-rich region. The arrangement of the genes was the same as that found in the ancestral insect. All the protein initiation codons are ATN, except for cox1 that begins with TCG (encoding S). The 22 tRNA anticodons of S. impatiens are consistent with those observed in Drosophila yakuba, and all form the typical cloverleaf structure, except for tRNA-Ser(AGN) that lacks the DHU arm. The mitochondrial genome of Sarcophaga presented will be valuable for resolving phylogenetic relationships within the family Sarcophagidae and the order Diptera, and could be used to identify favourable genetic markers for species identifications for forensic purposes.

Do RNA viruses require genome cyclisation for replication?

Complementary sequences at the 5′ and 3′ ends of the dengue virus RNA genome are essential for viral replication, and are believed to cyclise the genome through long-range base pairing in cis. Although consistent with evidence in the literature, this view neglects possible biologically active multimeric forms that are equally consistent with the data. Here, we propose alternative multimeric structures, and suggest that multigenome noncovalent concatemers are more likely to exist under cellular conditions than single cyclised monomers. Concatemers provide a plausible mechanism for the dengue virus to overcome the single-stranded (+)-sense RNA virus dilemma, and can potentially assist genome transport from the virus-induced vesicles into the cytosol.

Population Dynamics of RNA viruses

Between 50 and 100 million people are infected with dengue viruses each year and more than 100,000 of these die. Dr Choudhury has demonstrated that populations of dengue viruses in individual patients are genetically and functionally very diverse and that this diversity changes significantly at the time of major outbreaks of disease. The results of his studies may inform strategies which will make dengue vaccines far more effective.

Evolutionary genetics and human assisted movement of a globally invasive pest (Russian wheat aphid : Diuraphis noxia)

This PhD study has examined the population genetics of the Russian wheat aphid (RWA, Diuraphis noxia), one of the world’s most invasive agricultural pests, throughout its native and introduced global range. Firstly, this study investigated the geographic distribution of genetic diversity within and among RWA populations in western China. Analysis of mitochondrial data from 18 sites provided evidence for the long-term existence and expansion of RWAs in western China. The results refute the hypothesis that RWA is an exotic species only present in China since 1975. The estimated date of RWA expansion throughout western China coincides with the debut of wheat domestication and cultivation practices in western Asia in the Holocene. It is concluded that western China represents the limit of the far eastern native range of this species. Analysis of microsatellite data indicated high contemporary gene flow among northern populations in western China, while clear geographic isolation between northern and southern populations was identified across the Tianshan mountain range and extensive desert regions. Secondly, this study analyzed the worldwide pathway of invasion using both microsatellite and endosymbiont genetic data. Individual RWAs were obtained from native populations in Central Asia and the Middle East and invasive populations in Africa and the Americas. Results indicated two pathways of RWA invasion from 1) Syria in the Middle East to North Africa and 2) Turkey to South Africa, Mexico and then North and South America. Very little clone diversity was identified among invasive populations suggesting that a limited founder event occurred together with predominantly asexual reproduction and rapid population expansion. The most likely explanation for the rapid spread (within two years) from South Africa to the New World is by human movement, probably as a result of the transfer of wheat breeding material. Furthermore, the mitochondrial data revealed the presence of a universal haplotype and it is proposed that this haplotype is representative of a wheat associated super-clone that has gained dominance worldwide as a result of the widespread planting of domesticated wheat. Finally, this study examined salivary gland gene diversity to determine whether a functional basis for RWA invasiveness could be identified. Peroxidase DNA sequence data were obtained for a selection of worldwide RWA samples. Results demonstrated that most native populations were polymorphic while invasive populations were monomorphic, supporting previous conclusions relating to demographic founder effects in invasive populations. Purifying selection most likely explains the existence of a universal allele present in Middle Eastern populations, while balancing selection was evident in East Asian populations. Selection acting on the peroxidase gene may provide an allele-dependent advantage linked to the successful establishment of RWAs on wheat, and ultimately their invasion potential. In conclusion, this study is the most comprehensive molecular genetic investigation of RWA population genetics undertaken to date and provides significant insights into the source and pathway of global invasion and the potential existence of a wheat-adapted genotype that has colonised major wheat growing countries worldwide except for Australia. This research has major biosecurity implications for Australia’s grain industry.

Molecular genetics of migraine

Migraine is a common neurological disorder with a significantly heritable component. It is a complex disease and despite numerous molecular genetic studies, the exact pathogenesis causing the neurological disturbance remains poorly understood. Although several known molecular mechanisms have been associated with an increased risk for developing migraine, there remains significant scope for future studies. The majority of studies have investigated the most plausible candidate genes involved in common migraine pathogenesis utilising criteria that takes into account a combination of physiological functionality in conjunction with regions of genomic association. Thus, far genes involved in neurological, vascular or hormonal pathways have been identified and investigated on this basis. Genome-wide association studies (GWAS) studies have helped to identify novel regions that may be associated with migraine and have aided in providing the basis for further molecular investigations. However, further studies utilising sequencing technologies are required to characterise the genetic basis for migraine.