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.

Genomic evidence for the parallel evolution of coastal forms in the Senecio lautus complex

Instances of parallel ecotypic divergence where adaptation to similar conditions repeatedly cause similar phenotypic changes in closely related organisms are useful for studying the role of ecological selection in speciation. Here we used a combination of traditional and next generation genotyping techniques to test for the parallel divergence of plants from the Senecio lautus complex, a phenotypically variable groundsel that has adapted to disparate environments in the South Pacific. Phylogenetic analysis of a broad selection of Senecio species showed that members of the S. lautus complex form a distinct lineage that has diversified recently in Australasia. An inspection of thousands of polymorphisms in the genome of 27 natural populations from the S. lautus complex in Australia revealed a signal of strong genetic structure independent of habitat and phenotype. Additionally, genetic differentiation between populations was correlated with the geographical distance separating them, and the genetic diversity of populations strongly depended on geographical location. Importantly, coastal forms appeared in several independent phylogenetic clades, a pattern that is consistent with the parallel evolution of these forms. Analyses of the patterns of genomic differentiation between populations further revealed that adjacent populations displayed greater genomic heterogeneity than allopatric populations and are differentiated according to variation in soil composition. These results are consistent with a process of parallel ecotypic divergence in face of gene flow.

Population structure and patterns of genetic variation in a pearl oyster (Pinctada radiata) native to the Arabian Gulf

The study assessed natural levels and patterns of genetic variation in Arabian Gulf populations of a native pearl oyster to define wild population structure considering potential intrinsic and extrinsic factors that could influence any wild structure detected. The study was also the first attempt to develop microsatellite markers and to generate a genome survey sequence (GSS) dataset for the target species using next generation sequencing technology. The partial genome dataset generated has potential biotechnological applications and for pearl oyster farming in the future.

A population genomics overview of the neuronal nitric oxide synthase (nNOS) gene and its relationship to migraine susceptibility

The ubiquitous chemical messenger molecule nitric oxide (NO) has been implicated in a diverse range of biological activities including neurotransmission, smooth muscle motility and mediation of nociception. Endogenous synthesis of NO by the neuronal isoform of the nitric oxide synthase gene family has an essential role within the central and peripheral nervous systems in addition to the autonomic innervation of cerebral blood vessels. To investigate the potential role of NO and more specifically the neuronal nitric oxide synthase (nNOS) gene in migraine susceptibility, we investigated two microsatellite repeat variants residing within the 5′ and 3′ regions of the nNOS gene. Population genomic evaluation of the two nNOS repeat variants indicated significant linkage disequilibrium between the two loci. Z-DNA conformational sequence structures within the 5′ region of the nNOS gene have the potential to enhance or repress gene promoter activity. We suggest that genetic analysis of this 5′ repeat variant is the more functional variant expressing gene wide information that could affect endogenous NO synthesis and potentially result in diseased states. However, no association with migraine (with or without aura) was seen in our extensive case-control cohort (n = 579 affected with matched controls), when both the 5′ and 3′ genetic variants were investigated.

A genetic variant located in miR-423 is associated with reduced breast cancer risk

Background/Aim: Since microRNAs (miRNAs) act as translational regulators of multiple genes, single nucleotide polymorphisms (SNP) in them can have potentially wide-ranging effects. Using an association approach, this research examined the effects of the rs6505162 SNP, an A>C polymorphism located in the premiRNA region of miR-423, on breast cancer development. Materials and Methods: Caucasian Australian women with breast cancer and controls matched for age and ethnicity were genotyped for rs6505162 and their genotypic and allelic frequencies analysed for significant differences. Results: Analysis indicated that there were significant differences between the case and control populations (χ 2=6.70, p=0.035), with the CC genotype conferring reduced risk of breast cancer development (odds ratio=0.50 95% confidence interval=0.27-0.92, p=0.03). Conclusion: Further functional research is required to determine the mechanism of action of this SNP on miRNA function.

Analysis of expressed sequence tags from Actinidia : Applications of a cross species EST database for gene discovery in the areas of flavor, health, color and ripening

Background Kiwifruit (Actinidia spp.) are a relatively new, but economically important crop grown in many different parts of the world. Commercial success is driven by the development of new cultivars with novel consumer traits including flavor, appearance, healthful components and convenience. To increase our understanding of the genetic diversity and gene-based control of these key traits in Actinidia, we have produced a collection of 132,577 expressed sequence tags (ESTs). Results The ESTs were derived mainly from four Actinidia species (A. chinensis, A. deliciosa, A. arguta and A. eriantha) and fell into 41,858 non redundant clusters (18,070 tentative consensus sequences and 23,788 EST singletons). Analysis of flavor and fragrance-related gene families (acyltransferases and carboxylesterases) and pathways (terpenoid biosynthesis) is presented in comparison with a chemical analysis of the compounds present in Actinidia including esters, acids, alcohols and terpenes. ESTs are identified for most genes in color pathways controlling chlorophyll degradation and carotenoid biosynthesis. In the health area, data are presented on the ESTs involved in ascorbic acid and quinic acid biosynthesis showing not only that genes for many of the steps in these pathways are represented in the database, but that genes encoding some critical steps are absent. In the convenience area, genes related to different stages of fruit softening are identified. Conclusion This large EST resource will allow researchers to undertake the tremendous challenge of understanding the molecular basis of genetic diversity in the Actinidia genus as well as provide an EST resource for comparative fruit genomics. The various bioinformatics analyses we have undertaken demonstrates the extent of coverage of ESTs for genes encoding different biochemical pathways in Actinidia.

Characterisation of a glutathione reductase gene and its genetic locus from pea (Pisum sativum L.)

A cDNA encoding the chloroplast/mitochondrial form of glutathione reductase (GR:EC 1,6,4,2) from pea (Pisum sativum L.) was used to map a single GR locus, named GORI. In two domesticated genotypes of pea (cv, Birte and JI 399) it is likely that the GORI locus contains a single gene. However, in a semi-domesticated land race of pea sequences were detected but closely related sets of GR gene sequences were in JI 281 represent either a second intact gene or a partial or pseudogene copy. A GR gene was cloned from ev. Birte, sequenced and its structure analysed. No features of the transcription or structure of the gene suggested a mechanism for generating any more than one form of . From these data plus previously published biochemical evidence was suggested a second, distinct gene encoding for the cytosolic form of GR should be present in peas. The GORI-encoded GR mRNA can be detected in all main organs of the plant and no alternative spliced species was present which could perhaps account for the generation of multiple isoforms of GR. The mismatch between the number of charge-separable isoforms in pea and the proposed number suggests that different GR isoforms arise by some form of post-transnational modification.

Pattern analysis of the diversity of morphological plant attributes and herbage yield in a world collection of white clover (Trifolium repens L.) germplasm characterised in a summer moisture stress environment of Australia

Information on the variation available for different plant attributes has enabled germplasm collections to be effectively utilised in plant breeding. A world sourced collection of white clover germplasm has been developed at the White Clover Resource Centre at Glen Innes, New South Wales. This collection of 439 accessions was characterised under field conditions as a preliminary study of the genotypic variation for morphological attributes; stolon density, stolon branching, number of nodes. number of rooted nodes, stolon thickness, internode length, leaf length, plant height and plant spread, together with seasonal herbage yield. Characterisation was conducted on different batches of germplasm (subsets of accessions taken from the complete collection) over a period of five years. Inclusion of two check cultivars, Haifa and Huia, in each batch enabled adjustment of the characterisation data for year effects and attribute-by-year interaction effects. The component of variance for seasonal herbage yield among batches was large relative to that for accessions. Accession-by-experiment and accession-by-season interactions for herbage yield were not detected. Accession mean repeatability for herbage yield across seasons was intermediate (0.453). The components of genotypic variance among accessions for all attributes, except plant height, were larger than their respective standard errors. The estimates of accession mean repeatability for the attributes ranged from low (0.277 for plant height) to intermediate (0.544 for internode length). Multivariate techniques of clustering and ordination were used to investigate the diversity present among the accessions in the collection. Both cluster analysis and principal component analysis suggested that seven groups of accessions existed. It was also proposed from the pattern analysis results that accessions from a group characterised by large leaves, tall plants and thick stolons could be crossed with accessions from a group that had above average stolon density and stolon branching. This material could produce breeding populations to be used in recurrent selection for the development of white clover cultivars for dryland summer moisture stress environments in Australia. The germplasm collection was also found to be deficient in genotypes with high stolon density, high number of branches high number of rooted nodes and large leaves. This warrants addition of new germplasm accessions possessing these characteristics to the present germplasm collection.

Common variants in the CYP2C19 gene are associated with susceptibility to endometriosis

Objective: To follow-up previous studies highlighting a possible role for cytochrome P450, family 2, subfamily C, 19 (CYP2C19) in susceptibility to endometriosis by searching for additional variants in the CYP2C19 gene that may be associated with the disease. Design Case-control study. Setting Academic research. Subject(s) The cases comprised 2,271 women with surgically confirmed endometriosis; the controls comprised 939 women with self-report of no endometriosis and 1,770 unscreened population samples. Intervention(s) Sequencing of the CYP2C19 region and follow-up of 80 single nucleotide polymorphisms (SNPs) in two case-control samples. Main Outcome Measure(s) Allele frequency differences between cases and controls. Result(s) Sequencing of the CYP2C19 gene region resulted in the detection of a large number of known and novel SNPs. Genotyping of 80 polymorphic SNPs in 901 endometriosis cases and 939 controls resulted in study-wide significant association signals for SNPs in moderate or complete linkage disequilibrium with rs4244285, a functional SNP in exon 5 that abrogates CYP2C19 function through the creation of an alternative splice site. Evidence of association was also detected for another functional SNP in the CYP2C19 promoter, rs12248560, which was highlighted in our previous study. Conclusion(s) Functional variants in CYP2C19 may contribute to endometriosis susceptibility in both familial and sporadic cases. © 2014 by American Society for Reproductive Medicine.

Genetic structure of the Korean black scraper Thamnaconus modestus inferred from microsatellite marker analysis

The Korean black scraper, Thamnaconus modestus, is one of the most economically important maricultural fish species in Korea. However, the annual catch of this fish has been continuously declining over the past several decades. In this study, the genetic diversity and relationships among four wild populations and two hatchery stocks of Korean black scraper were assessed based on 16 microsatellite (MS) markers. A total of 319 different alleles were detected over all loci with an average of 19.94 alleles per locus. The hatchery stocks [mean number of alleles (N A) = 12, allelic richness (A R) = 12, expected heterozygosity (He) = 0.834] showed a slight reduction (P > 0.05) in genetic variability in comparison with wild populations (mean N A = 13.86, A R = 12.35, He = 0.844), suggesting a sufficient level of genetic variation in the hatchery populations. Similarly low levels of inbreeding and significant Hardy–Weinberg equilibrium deviations were detected in both wild and hatchery populations. The genetic subdivision among all six populations was low but significant (overall F ST = 0.008, P < 0.01). Pairwise F ST, a phylogenetic tree, and multidimensional scaling analysis suggested the existence of three geographically structured populations based on different sea basin origins, although the isolation-by-distance model was rejected. This result was corroborated by an analysis of molecular variance. This genetic differentiation may result from the co-effects of various factors, such as historical dispersal, local environment and ocean currents. These three geographical groups can be considered as independent management units. Our results show that MS markers may be suitable not only for the genetic monitoring of hatchery stocks but also for revealing the population structure of Korean black scraper populations. These results will provide critical information for breeding programs, the management of cultured stocks and the conservation of this species.

Association of variants at 1q32 and STAT3 with ankylosing spondylitis suggests genetic overlap with Crohn's disease

Ankylosing spondylitis (AS) is a common inflammatory arthritic condition. Overt inflammatory bowel disease (IBD) occurs in about 10% of AS patients, and in addition 70% of AS cases may have subclinical terminal ileitis. Spondyloarthritis is also common in IBD patients. We therefore tested Crohn's disease susceptibility genes for association with AS, aiming to identify pleiotropic genetic associations with both diseases. Genotyping was carried out using Sequenom and Applied Biosystems TaqMan and OpenArray technologies on 53 markers selected from 30 Crohn's disease associated genomic regions. We tested genotypes in a population of unrelated individual cases (n = 2,773) and controls (n = 2,215) of white European ancestry for association with AS. Statistical analysis was carried out using a Cochran-Armitage test for trend in PLINK. Strong association was detected at chr1q32 near KIF21B (rs11584383, P = 1.66 x 10-10, odds ratio (OR) = 0.74, 95% CI:0.68-0.82). Association with disease was also detected for 2 variants within STAT3 (rs6503695, P = 4.6×10-4. OR = 0.86 (95% CI:0.79-0.93); rs744166, P = 2.6×10-5, OR = 0.84 (95% CI:0.77-0.91)). Association was confirmed for IL23R (rs11465804, P = 1.2×10-5, OR = 0.65 (95% CI:0.54-0.79)), and further associations were detected for IL12B (rs10045431, P = 5.261025, OR = 0.83 (95% CI:0.76-0.91)), CDKAL1 (rs6908425, P = 1.1×10-4, OR = 0.82 (95% CI:0.74-0.91)), LRRK2/MUC19 (rs11175593, P = 9.9×10-5, OR = 1.92 (95% CI: 1.38-2.67)), and chr13q14 (rs3764147, P = 5.9×10-4, OR = 1.19 (95% CI: 1.08-1.31)). Excluding cases with clinical IBD did not significantly affect these findings. This study identifies chr1q32 and STAT3 as ankylosing spondylitis susceptibility loci. It also further confirms association for IL23R and detects suggestive association with another 4 loci. STAT3 is a key signaling molecule within the Th17 lymphocyte differentiation pathway and further enhances the case for a major role of this T-lymphocyte subset in ankylosing spondylitis. Finally these findings suggest common aetiopathogenic pathways for AS and Crohn's disease and further highlight the involvement of common risk variants across multiple diseases.

Ankylosing spondylitis is associated with the anthrax toxin receptor 2 gene (ANTXR2)

Objectives ANTXR2 variants have been associated with ankylosing spondylitis (AS) in two previous genome-wide association studies (GWAS) (p∼9×10-8). However, a genome-wide significant association (p<5×10-8) was not observed. We conducted a more comprehensive analysis of ANTXR2 in an independent UK sample to confirm and refine this association. Methods A replication study was carried out with 2978 cases and 8365 controls. Then, these were combined with non-overlapping samples from the two previous GWAS in a meta-analysis. Human leukocyte antigen (HLA)-B27 stratification was also performed to test for ANTXR2-HLA-B27 interaction. Results Out of nine single nucleotide polymorphisms (SNP) in the study, five SNPs were nominally associated (p<0.05) with AS in the replication dataset. In the meta-analysis, eight SNPs showed evidence of association, the strongest being with rs12504282 (OR=0.88, p=6.7×10-9). Seven of these SNPs showed evidence for association in the HLA-B27-positive subgroup, but none was associated with HLA-B27-negative AS. However, no statistically significant interaction was detected between HLA-B27 and ANTXR2 variants. Conclusions ANTXR2 variants are clearly associated with AS. The top SNPs from two previous GWAS (rs4333130 and rs4389526) and this study (rs12504282) are in strong linkage disequilibrium (r2≥0.76). All are located near a putative regulatory region. Further studies are required to clarify the role played by these ANTXR2 variants in AS.

Genetic studies in osteoporosis: The end of the beginning

Osteoporosis and disorders of bone fragility are highly heritable, but despite much effort the identities of few of the genes involved has been established. Recent developments in genetics such as genome-wide association studies are revolutionizing research in this field, and it is likely that further contributions will be made through application of next-generation sequencing technologies, analysis of copy number variation polymorphisms, and high-throughput mouse mutagenesis programs. This article outlines what we know about osteoporosis genetics to date and the probable future directions of research in this field.

Evidence of genetic association between TNFRSF1A encoding the p55 tumour necrosis factor receptor, and ankylosing spondylitis in UK Caucasians

Objectives: To replicate the possible genetic association between ankylosing spondylitis (AS) and TNFRSF1A. Methods: TNFRSF1A was re-sequenced in 48 individuals with AS to identify novel polymorphisms. Nine single nucleotide polymorphisms (SNPs) in TNFRSF1A and 5 SNPs in the neighbouring gene SCNN1A were genotyped in 1604 UK Caucasian individuals with AS and 1019 matched controls. An extended study was implemented using additional genotype data on 8 of these SNPs from 1400 historical controls from the 1958 British Birth Cohort. A meta-analysis of previously published results was also undertaken. Results: One novel variant in intron 6 was identified but no new coding variants. No definite associations were seen in the initial study but in the extended study there were weak associations with rs4149576 (p=0.04) and rs4149577 (p=0.007). In the metaanalysis consistent, somewhat stronger associations were seen with rs4149577 (p=0.002) and rs4149578 (p=0.006). Conclusions: These studies confirm the weak genetic associations between AS and TNFRSF1A. In view of the previously reported associations of TNFRSF1A with AS, in Caucasians and Chinese, and the biological plausibility of this candidate gene, replication of this finding in well powered studies is clearly indicated.

Cardiometabolic effects of genetic upregulation of the interleukin 1 receptor antagonist: A Mendelian randomisation analysis

Background To investigate potential cardiovascular and other effects of long-term pharmacological interleukin 1 (IL-1) inhibition, we studied genetic variants that produce inhibition of IL-1, a master regulator of inflammation. Methods We created a genetic score combining the effects of alleles of two common variants (rs6743376 and rs1542176) that are located upstream of IL1RN, the gene encoding the IL-1 receptor antagonist (IL-1Ra; an endogenous inhibitor of both IL-1α and IL-1β); both alleles increase soluble IL-1Ra protein concentration. We compared effects on inflammation biomarkers of this genetic score with those of anakinra, the recombinant form of IL-1Ra, which has previously been studied in randomised trials of rheumatoid arthritis and other inflammatory disorders. In primary analyses, we investigated the score in relation to rheumatoid arthritis and four cardiometabolic diseases (type 2 diabetes, coronary heart disease, ischaemic stroke, and abdominal aortic aneurysm; 453 411 total participants). In exploratory analyses, we studied the relation of the score to many disease traits and to 24 other disorders of proposed relevance to IL-1 signalling (746 171 total participants). Findings For each IL1RN minor allele inherited, serum concentrations of IL-1Ra increased by 0·22 SD (95% CI 0·18–0·25; 12·5%; p=9·3 × 10−33), concentrations of interleukin 6 decreased by 0·02 SD (−0·04 to −0·01; −1·7%; p=3·5 × 10−3), and concentrations of C-reactive protein decreased by 0·03 SD (−0·04 to −0·02; −3·4%; p=7·7 × 10−14). We noted the effects of the genetic score on these inflammation biomarkers to be directionally concordant with those of anakinra. The allele count of the genetic score had roughly log-linear, dose-dependent associations with both IL-1Ra concentration and risk of coronary heart disease. For people who carried four IL-1Ra-raising alleles, the odds ratio for coronary heart disease was 1·15 (1·08–1·22; p=1·8 × 10−6) compared with people who carried no IL-1Ra-raising alleles; the per-allele odds ratio for coronary heart disease was 1·03 (1·02–1·04; p=3·9 × 10−10). Per-allele odds ratios were 0·97 (0·95–0·99; p=9·9 × 10−4) for rheumatoid arthritis, 0·99 (0·97–1·01; p=0·47) for type 2 diabetes, 1·00 (0·98–1·02; p=0·92) for ischaemic stroke, and 1·08 (1·04–1·12; p=1·8 × 10−5) for abdominal aortic aneurysm. In exploratory analyses, we observed per-allele increases in concentrations of proatherogenic lipids, including LDL-cholesterol, but no clear evidence of association for blood pressure, glycaemic traits, or any of the 24 other disorders studied. Modelling suggested that the observed increase in LDL-cholesterol could account for about a third of the association observed between the genetic score and increased coronary risk. Interpretation Human genetic data suggest that long-term dual IL-1α/β inhibition could increase cardiovascular risk and, conversely, reduce the risk of development of rheumatoid arthritis. The cardiovascular risk might, in part, be mediated through an increase in proatherogenic lipid concentrations. Funding UK Medical Research Council, British Heart Foundation, UK National Institute for Health Research, National Institute for Health Research Cambridge Biomedical Research Centre, European Research Council, and European Commission Framework Programme 7.