Making better fertiliser decisions for the cropping regions of Australia

Collating old fertilizer trial data for development of a national database on crop responses.

Multispectral imaging and genomic mapping

Identifying QTLs linked to malting and protein levels in barley using multispectral image data.

Genomic inflation factors under polygenic inheritance

Population structure, including population stratification and cryptic relatedness, can cause spurious associations in genome-wide association studies (GWAS). Usually, the scaled median or mean test statistic for association calculated from multiple single-nucleotide-polymorphisms across the genome is used to assess such effects, and 'genomic control' can be applied subsequently to adjust test statistics at individual loci by a genomic inflation factor. Published GWAS have clearly shown that there are many loci underlying genetic variation for a wide range of complex diseases and traits, implying that a substantial proportion of the genome should show inflation of the test statistic. Here, we show by theory, simulation and analysis of data that in the absence of population structure and other technical artefacts, but in the presence of polygenic inheritance, substantial genomic inflation is expected. Its magnitude depends on sample size, heritability, linkage disequilibrium structure and the number of causal variants. Our predictions are consistent with empirical observations on height in independent samples of ~4000 and ~133,000 individuals.

Molecular mechanisms of cancer predisposition in HNPCC/Lynch syndrome

Hereditary non-polyposis colorectal carcinoma (HNPCC; Lynch syndrome) is among the most common hereditary cancers in man and a model of cancers arising through deficient DNA mismatch repair (MMR). It is inherited in a dominant manner with predisposing germline mutations in the MMR genes, mainly MLH1, MSH2, MSH6 and PMS2. Both copies of the MMR gene need to be inactivated for cancer development. Since Lynch syndrome family members are born with one defective copy of one of the MMR genes in their germline, they only need to acquire a so called second hit to inactivate the MMR gene. Hence, they usually develop cancer at an early age. MMR gene inactivation leads to accumulation of mutations particularly in short repeat tracts, known as microsatellites, causing microsatellite instability (MSI). MSI is the hallmark of Lynch syndrome tumors, but is present in approximately 15% of sporadic tumors as well. There are several possible mechanisms of somatic inactivation (i.e. the second hit ) of MMR genes, for instance deletion of the wild-type copy, leading to loss of heterozygosity (LOH), methylation of promoter regions necessary for gene transcription, or mitotic recombination or gene conversion. In the Lynch syndrome tumors carrying germline mutations in the MMR gene, LOH was found to be the most frequent mechanism of somatic inactivation in the present study. We also studied MLH1/MSH2 deletion carriers and found that somatic mutations identical to the ones in the germline occurred frequently in colorectal cancers and were also present in extracolonic Lynch syndrome-associated tumors. Chromosome-specific marker analysis implied that gene conversion, rather than mitotic recombination or deletion of the respective gene locus accounted for wild-type inactivation. Lynch syndrome patients are predisposed to certain types of cancers, the most common ones being colorectal, endometrial and gastric cancer. Gastric cancer and uroepithelial tumors of bladder and ureter were observed to be true Lynch syndrome tumors with MMR deficiency as the driving force of tumorigenesis. Brain tumors and kidney carcinoma, on the other hand, were mostly MSS, implying the possibility of alternative routes of tumor development. These results present possible implications in clinical cancer surveillance. In about one-third of families suspected of Lynch syndrome, mutations in MMR genes are not found, and we therefore looked for alternative mechanisms of predisposition. According to our results, large genomic deletions, mainly in MSH2, and germline epimutations in MLH1, together explain a significant fraction of point mutation-negative families suspected of Lynch syndrome and are associated with characteristic clinical and family features. Our findings have important implications in the diagnosis and management of Lynch syndrome families.

Methicillin Resistance In Staphylococci : Horizontal Transfer of Mobile Genetic Element (SCCmec) Between Staphylococcal Species

Regardless of the existence of antibiotics, infectious diseases are the leading causes of death in the world. Staphylococci cause many infections of varying severity, although they can also exist peacefully in many parts of the human body. Most often Staphylococcus aureus colonises the nose, and that colonisation is considered to be a risk factor for spread of this bacterium. S. aureus is considered to be the most important Staphylococcus species. It poses a challenge to the field of medicine, and one of the most problematic aspects is the drastic increase of the methicillin-resistant S. aureus (MRSA) strains in hospitals and community world-wide, including Finland. In addition, most of the clinical coagulase-negative staphylococcus (CNS) isolates express resistance to methicillin. Methicillin-resistance in S. aureus is caused by the mecA gene that encodes an extra penicillin-binding protein (PBP) 2a. The mecA gene is found in a mobile genomic island called staphylococcal chromosome cassette mec (SCCmec). The SCCmec consists of the mec gene and cassette chromosome recombinase (ccr)gene complexes. The areas of the SCCmec element outside the ccr and mec complex are known as the junkyard J regions. So far, eight types of SCCmec(SCCmec I- SCCmec VIII) and a number of variants have been described. The SCCmec island is an acquired element in S. aureus. Lately, it appears that CNS might be the storage place of the SCCmec that aid the S. aureus by providing it with the resistant elements. The SCCmec is known to exist only in the staphylococci. The aim of the present study was to investigate the horizontal transfer of SCCmec between the S. aureus and CNS. One specific aim was to study whether or not some methicillin-sensitive S. aureus (MSSA) strains are more inclined to receive the SCCmec than others. This was done by comparing the genetic background of clinical MSSA isolates in the health care facilities of the Helsinki and Uusimaa Hospital District in 2001 to the representatives of the epidemic MRSA (EMRSA) genotypes, which have been encountered in Finland during 1992-2004. Majority of the clinical MSSA strains were related to the EMRSA strains. This finding suggests that horizontal transfer of SCCmec from unknown donor(s) to several MSSA background genotypes has occurred in Finland. The molecular characteristics of representative clinical methicillin-resistant S. epidermidis (MRSE) isolates recovered in Finnish hospitals between 1990 and 1998 were also studied, examining their genetic relation to each other and to the internationally recognised MRSE clones as well, so as to ascertain the common traits between the SCCmec elements in MRSE and MRSA. The clinical MRSE strains were genetically related to each other; eleven PFGE types were associated with sequence type ST2 that has been identified world-wide. A single MRSE strain may possess two SCCmec types III and IV, which were recognised among the MRSA strains. Moreover, six months after the onset of an outbreak of MRSA possessing a SCCmec type V in a long-term care facility in Northern Finland (LTCF) in 2003, the SCCmec element of nasally carried methicillin-resistant staphylococci was studied. Among the residents of a LTCF, nasal carriage of MR-CNS was common with extreme diversity of SCCmec types. MRSE was the most prevalent CNS species. Horizontal transfer of SCCmec elements is speculated to be based on the sharing of SCCmec type V between MRSA and MRSE in the same person. Additionally, the SCCmec element of the clinical human S. sciuri isolates was studied. Some of the SCCmec regions were present in S. sciuri and the pls gene was common in it. This finding supports the hypothesis of genetic exchange happening between staphylococcal species. Evaluation of the epidemiology of methicillin-resistant staphylococcal colonisation is necessary in order to understand the apparent emergence of these strains and to develop appropriate control strategies. SCCmec typing is essential for understanding the emergence of MRSA strains from CNS, considering that the MR-CNS may represent the gene pool for the continuous creation of new SCCmec types from which MRSA might originate.

Hundreds of variants clustered in genomic loci and biological pathways affect human height

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

Reliable Means of Diagnosis and Serovar Determination of Blood-Borne Salmonella Strains: Quick PCR Amplification of Unique Genomic Loci by Novel Primer Sets

Typhoid fever is becoming an ever increasing threat in the developing countries. We have improved considerably upon the existing PCR-based diagnosis method by designing primers against a region that is unique to Salmonella enterica subsp. enterica serovar Typhi and Salmonella enterica subsp. enterica serovar Paratyphi A, corresponding to the STY0312 gene in S. Typhi and its homolog SPA2476 in S. Paratyphi A. An additional set of primers amplify another region in S. Typhi CT18 and S. Typhi Ty2 corresponding to the region between genes STY0313 to STY0316 but which is absent in S. Paratyphi A. The possibility of a false-negative result arising due to mutation in hypervariable genes has been reduced by targeting a gene unique to typhoidal Salmonella serovars as a diagnostic marker. The amplified region has been tested for genomic stability by amplifying the region from clinical isolates of patients from various geographical locations in India, thereby showing that this region is potentially stable. These set of primers can also differentiate between S. Typhi CT18, S. Typhi Ty2, and S. Paratyphi A, which have stable deletions in this specific locus. The PCR assay designed in this study has a sensitivity of 95% compared to the Widal test which has a sensitivity of only 63%. As observed, in certain cases, the PCR assay was more sensitive than the blood culture test was, as the PCR-based detection could also detect dead bacteria.

Integrating sorghum whole genome sequence information with a compendium of sorghum QTL studies reveals uneven distribution of QTL and of gene-rich regions with significant implications for crop improvement.

A comprehensive analysis was conducted using 48 sorghum QTL studies published from 1995 to 2010 to make information from historical sorghum QTL experiments available in a form that could be more readily used by sorghum researchers and plant breeders. In total, 771 QTL relating to 161 unique traits from 44 studies were projected onto a sorghum consensus map. Confidence intervals (CI) of QTL were estimated so that valid comparisons could be made between studies. The method accounted for the number of lines used and the phenotypic variation explained by individual QTL from each study. In addition, estimated centimorgan (cM) locations were calculated for the predicted sorghum gene models identified in Phytozome (JGI GeneModels SBI v1.4) and compared with QTL distribution genome-wide, both on genetic linkage (cM) and physical (base-pair/bp) map scales. QTL and genes were distributed unevenly across the genome. Heterochromatic enrichment for QTL was observed, with approximately 22% of QTL either entirely or partially located in the heterochromatic regions. Heterochromatic gene enrichment was also observed based on their predicted cM locations on the sorghum consensus map, due to suppressed recombination in heterochromatic regions, in contrast to the euchromatic gene enrichment observed on the physical, sequence-based map. The finding of high gene density in recombination-poor regions, coupled with the association with increased QTL density, has implications for the development of more efficient breeding systems in sorghum to better exploit heterosis. The projected QTL information described, combined with the physical locations of sorghum sequence-based markers and predicted gene models, provides sorghum researchers with a useful resource for more detailed analysis of traits and development of efficient marker-assisted breeding strategies.

White shark (Carcharodon carcharias) microsatellite genotypes (six loci) from regions along Australia's southern and eastern coastline.

Despite international protection of white sharks (Carcharodon carcharias), important conservation parameters such as abundance, population structure and genetic diversity are largely unknown. The tissue of 97 predominately juvenile white sharks sampled from spatially distant eastern and southwestern Australian coastlines was sequenced for the mitochondrial DNA (mtDNA) control region and genotyped with six nuclear-encoded microsatellite loci. MtDNA population structure was found between the eastern and southwestern coasts (FST = 0.142, p < 0.001), implying female natal philopatry. This concords with recent satellite and acoustic tracking findings which suggest the sustained presence of discrete east coast nursery areas. Furthermore, population subdivision was found between the same regions with biparentally inherited microsatellite markers (FST = 0.009, p <0.05), suggesting that males may also exhibit some degree of reproductive philopatry. Five sharks captured along the east coast had mtDNA haplotypes that resembled western Indian Ocean sharks more closely than Australian/New Zealand sharks, suggesting that transoceanic dispersal or migration resulting in breeding may occur sporadically. Our most robust estimate of contemporary genetic effective population size was low and below the threshold at which adaptive potential may be lost. For a variety of reasons, these contemporary estimates were at least one, possibly two orders of magnitude below our historical effective size estimates. Further population decline could expose these genetically isolated populations to detrimental genetic effects. Regional Australian white shark conservation management units should be implemented until genetic population structure, size and diversity can be investigated in more detail. Reference: Blower, D. C., Pandolfi, J. M., Gomez-Cabrera, M. del C., Bruce, B. D. & Ovenden, J. R. (In press - April 2012). Population genetics of Australian white sharks reveals fine-scale spatial structure, transoceanic dispersal events and low effective population sizes. Marine Ecology Progress Series.

Estimating winter crop area across seasons and regions using time-sequential MODIS imagery.

The wheat grain industry is Australia's second largest agricultural export commodity. There is an increasing demand for accurate, objective and near real-time crop production information by industry. The advent of the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite platform has augmented the capability of satellite-based applications to capture reflectance over large areas at acceptable pixel scale, cost and accuracy. The use of multi-temporal MODIS-enhanced vegetation index (EVI) imagery to determine crop area was investigated in this article. Here the rigour of the harmonic analysis of time-series (HANTS) and early-season metric approaches was assessed when extrapolating over the entire Queensland (QLD) cropping region for the 2005 and 2006 seasons. Early-season crop area estimates, at least 4 months before harvest, produced high accuracy at pixel and regional scales with percent errors of -8.6% and -26% for the 2005 and 2006 seasons, respectively. In discriminating among crops at pixel and regional scale, the HANTS approach showed high accuracy. The errors for specific area estimates for wheat, barley and chickpea were 9.9%, -5.2% and 10.9% (for 2005) and -2.8%, -78% and 64% (for 2006), respectively. Area estimates of total winter crop, wheat, barley and chickpea resulted in coefficient of determination (R(2)) values of 0.92, 0.89, 0.82 and 0.52, when contrasted against the actual shire-scale data. A significantly high coefficient of determination (0.87) was achieved for total winter crop area estimates in August across all shires for the 2006 season. Furthermore, the HANTS approach showed high accuracy in discriminating cropping area from non-cropping area and highlighted the need for accurate and up-to-date land use maps. The extrapolability of these approaches to determine total and specific winter crop area estimates, well before flowering, showed good utility across larger areas and seasons. Hence, it is envisaged that this technology might be transferable to different regions across Australia.

Potato virus A as a heterologous protein expression tool in plants

Viruksien käyttö tuotekehityksen ja tutkimuksen vaatimien proteiinien tuottamiseen, syötävien rokotteiden kehittämiseen ja geeniterapiaan edustavat kasvavia biotekniikan sovellusalueita. Perunan A-virus (PVA) kuuluu potyviruksiin, joiden proteiinit tuotetaan aluksi yhtenä suurena molekyylinä, joka pilkotaan yksittäisiksi proteiineiksi viruksen itsensä tuottamilla entsyymeillä. Siten virusgenomiin lisätty vieras geeni käännetään proteiiniksi virusproteiinien mukana. Lopputuloksena kaikkia proteiineja tuotetaan kasvisoluissa samansuuruinen määrä. Lisäksi, viruksen proteiinikuoren koontimekanismi sallii perintöaineksen merkittävän lisäyksen ilman että viruksen tartutuskyky merkittävästi heikkenee. Koska virus monistuu ja leviää koko kasviin, jo melko pieni määrä kasveja riittää huomattavan proteiinimäärän tuottamiseen esimerkiksi säännösten mukaisessa kasvihuoneessa. Tämän työn tarkoituksena oli muuntaa PVA:n genomia siten, että virus soveltuisi yhden vieraan proteiinin tai useiden erilaisten proteiinien samanaikaiseen tuottamiseen kasveissa. Aluksi kokeiltiin viruksen replikaasia ja kuoriproteiinia koodaavien genomialueiden välistä kohtaa ja ihmisestä peräisi olevaa geeniä, joka tuotti S-COMT-entsyymiä (katekoli-O-metyylitransferaasi). Sen aktiivisuuden rajoittaminen auttaa Parkinsonintaudin hoidossa. Kasvissa tuotettua S-COMT:ia voitaisiin käyttää lääkekehityksessä estolääkkeiden testaukseen. Kahden viikon kuluttua tartutuksesta tupakan lehdissä oli entsymaattisesti aktiivista S-COMT:ia n. 1 % lehden liukoisista proteiineista. PVA:n P1-proteiinia koodaavalta alueelta oli paikannettu kohta, johon ehkä voitaisiin siirtää vieras geeni. Asia varmistettiin siirtämällä tähän kohtaan meduusan geeni, joka tuottaa UV-valossa vihreänä fluoresoivaa proteiinia (GFP). GFP-geeniä kantava PVA levisi kasvissa ja lisääntyi n. 30-50 %:iin viruksen normaalista pitoisuudesta. Koko kasvi fluoresoi vihreänä UV-valossa. Vieras geeni voidaan sijoittaa myös potyviruksen P1- ja HCpro-proteiineja koodaavien alueiden väliin. Samaan PVA-genomiin siirrettiin kolme geeniä, yksi kuhunkin kolmesta kloonauskohdasta: GFP-geeni P1:n sisälle, merivuokon lusiferaasigeeni P1/HCpro-kohtaan ja bakteerin beta-glukuronidaasigeeni (GUS) replikaasi/kuoriproteiini-kohtaan. Virusgenomin ja itse viruksen pituudet kasvoivat 38 %, mutta virus säilytti tartutuskykynsä. Se levisi kasveissa saavuttaen n. 15 % viruksen normaalista pitoisuudesta. Kaikki kolme vierasta proteiinia esiintyivät lehdissä aktiivisina.

Sweet potato chlorotic stunt virus: Studies on viral synergism and suppression of RNA silencing

The studies presented in this thesis aimed to a better understanding of the molecular biology of Sweet potato chlorotic stunt virus (SPCSV, Crinivirus, Closteroviridae) and its role in the development of synergistic viral diseases. The emphasis was on the severe sweet potato virus disease (SPVD) that results from a synergistic interaction of SPCSV and Sweet potato feathery mottle virus (SPFMV, Potyvirus, Potyviridae). SPVD is the most important disease affecting sweetpotato. It is manifested as a significant increase in symptom severity and SPFMV titres. This is accompanied by a dramatic sweetpotato yield reduction. SPCSV titres remain little affected in the diseased plants. Viral synergistic interactions have been associated with the suppression of an adaptive general defence mechanism discovered in plants and known as RNA silencing. In the studies of this thesis two novel proteins (RNase3 and p22) identified in the genome of a Ugandan SPCSV isolate were shown to be involved in suppression of RNA silencing. RNase3 displayed a dsRNA-specific endonuclease activity that enhanced the RNA-silencing suppression activity of p22. Comparative analyses of criniviral genomes revealed variability in the gene content at the 3´end of the genomic RNA1. Molecular analyses of different isolates of SPCSV indicated a marked intraspecific heterogeneity in this region where the p22 and RNase3 genes are located. Isolates of the East African strain of SPCSV from Tanzania and Peru and an isolate from Israel were missing a 767-nt fragment that included the p22 gene. However, regardless of the absence of p22, all SPCSV isolates acted synergistically with SPFMV in co-infected sweetpotato, enhanced SPFMV titres and caused SPVD. These results showed that p22 is dispensable for development of SPVD. The role of RNase3 in SPVD was then studied by generating transgenic plants expressing the RNase3 protein. These plants had increased titres of SPFMV (ca. 600-fold higher in comparison with nontransgenic plants) 2-3 weeks after graft inoculation and displayed the characteristic SPVD symptoms. RNA silencing suppression (RSS) activity of RNase3 was detected in agroinfiltrated leaves of Nicotiana bethamiana. In vitro studies showed that RNase3 was able to cleave small interferring RNAs (siRNA) to products of ~14-nt. The data thus identified RNase3 as a suppressor of RNA silencing able to cleave siRNAs. RNase3 expression alone was sufficient for breaking down resistance to SPFMV in sweetpotato and for the development of SPVD. Similar RNase III-like genes exist in animal viruses which points out a novel and possibly more general mechanism of RSS by viruses. A reproducible method of sweetpotato transformation was used to target RNA silencing against the SPCSV polymerase region (RdRp) with an intron-spliced hairpin construct. Hence, engineered resistance to SPCSV was obtained. Ten out of 20 transgenic events challenged with SPCSV alone showed significantly reduced virus titres. This was however not sufficient to prevent SPVD upon coinfection with SPFMV. Immunity to SPCSV seems to be required to control SPVD and targeting of different SPCSV regions need to be assessed in further studies. Based on the identified key role of RNase3 in SPVD the possibility to design constructs that target this gene might prove more efficient in future studies.

A comparison of DNA pools constructed following whole genome amplification for two-stage SNP genotyping designs

Genotyping in DNA pools reduces the cost and the time required to complete large genotyping projects. The aim of the present study was to evaluate pooling as part of a strategy for fine mapping in regions of significant linkage. Thirty-nine single nucleotide polymorphisms (SNPs) were analyzed in two genomic DNA pools of 384 individuals each and results compared with data after typing all individuals used in the pools. There were no significant differences using data from either 2 or 8 heterozygous individuals to correct frequency estimates for unequal allelic amplification. After correction, the mean difference between estimates from the genomic pool and individual allele frequencies was .033. A major limitation of the use of DNA pools is the time and effort required to carefully adjust the concentration of each individual DNA sample before mixing aliquots. Pools were also constructed by combining DNA after Multiple Displacement Amplification (MDA). The MDA pools gave similar results to pools constructed after careful DNA quantitation (mean difference from individual genotyping .040) and MDA provides a rapid method to generate pools suitable for some applications. Pools provide a rapid and cost-effective screen to eliminate SNPs that are not polymorphic in a test population and can detect minor allele frequencies as low as 1% in the pooled samples. With current levels of accuracy, pooling is best suited to an initial screen in the SNP validation process that can provide high-throughput comparisons between cases and controls to prioritize SNPs for subsequent individual genotyping.