The GP problem: quantifying gene-to-phenotype relationships

In this paper we refer to the gene-to-phenotype modeling challenge as the GP problem. Integrating information across levels of organization within a genotype-environment system is a major challenge in computational biology. However, resolving the GP problem is a fundamental requirement if we are to understand and predict phenotypes given knowledge of the genome and model dynamic properties of biological systems. Organisms are consequences of this integration, and it is a major property of biological systems that underlies the responses we observe. We discuss the E(NK) model as a framework for investigation of the GP problem and the prediction of system properties at different levels of organization. We apply this quantitative framework to an investigation of the processes involved in genetic improvement of plants for agriculture. In our analysis, N genes determine the genetic variation for a set of traits that are responsible for plant adaptation to E environment-types within a target population of environments. The N genes can interact in epistatic NK gene-networks through the way that they influence plant growth and development processes within a dynamic crop growth model. We use a sorghum crop growth model, available within the APSIM agricultural production systems simulation model, to integrate the gene-environment interactions that occur during growth and development and to predict genotype-to-phenotype relationships for a given E(NK) model. Directional selection is then applied to the population of genotypes, based on their predicted phenotypes, to simulate the dynamic aspects of genetic improvement by a plant-breeding program. The outcomes of the simulated breeding are evaluated across cycles of selection in terms of the changes in allele frequencies for the N genes and the genotypic and phenotypic values of the populations of genotypes.

Towards a taxonomic revision of the genus Echinococcus

Echinococcus remains a significant public health problem worldwide and, in several regions, the aetiological agents of cystic hydatid disease/echinococcosis are extending their range. The taxonomy of Echinococcus has been a controversial issue for decades, but the outcome of recent molecular epidemiological studies has served to reinforce proposals made ten years ago to revise the taxonomy of Echinococcus. A formal nomenclature is essential for effective communication, and provides the stability that underpins epidemiological investigations. It will also serve to recognize the contribution of early taxonomists.

Streptococcus pyogenes prtFII, but not sfbI, sfbII or fbp54, is represented more frequently among invasive-disease isolates of tropical Australia

Streptococcus pyogenes (group A streptococcus) strains may express several distinct fibronectin-binding proteins (FBPs) which are considered as major streptococcal adhesins. Of the FBPs, SfbI was shown in vitro to promote internalization of the bacterium into host cells and has been implicated in persistence. In the tropical Northern Territory, where group A streptococcal infection is common, multiple genotypes of the organism were found among isolates from invasive disease cases and no dominant strains were observed. To determine whether any FBPs is associated with invasive disease propensity of S. pyogenes, we have screened streptococcal isolates from bacteraemic and necrotizing fasciitis patients and isolates from uncomplicated infections for genetic endowment of 4 FBPs. No difference was observed in the distribution of sfbII, fbp54 and sfbI between the blood isolates' and isolates from uncomplicated infection. We conclude that the presence of sfbI does not appear to promote invasive diseases, despite its association with persistence. We also show a higher proportion of group A streptococcus strains isolated from invasive disease cases possess prtFII when compared to strains isolated from non-invasive disease cases. We suggest that S. pyogenes may recruit different FBPs for different purposes.

Complete mitochondrial genomes confirm the distinctiveness of the horse-dog and sheep-dog strains of Echinococcus granulosus

Unlike other members of the genus, Echinococcus granulosus is known to exhibit considerable levels of variation in biology, physiology and molecular genetics. Indeed, some of the taxa regarded as 'genotypes' within E. granulosus might be sufficiently distinct as to merit specific status. Here, complete mitochondrial genomes are presented of 2 genotypes of E. granulosus (G1-sheep-dog strain: G4-horse-dog strain) and of another taeniid cestode, Taenia crassiceps. These genomes are characterized and compared with those of Echinococcus multilocularis and Hymenolepis diminuta. Genomes of all the species are very similar in structure, length and base-composition. Pairwise comparisons of concatenated protein-coding genes indicate that the G1 and G4 genotypes of E. granulosus are almost as distant from each other as each is from a distinct species, E. multilocularis. Sequences for the variable genes atp6 and nad3 were obtained from additional genotypes of E. granulosus, from E. vogeli and E. oligarthrus. Again, pairwise comparisons showed the distinctiveness of the G1 and G4 genotypes. Phylogenetic analyses of concatenated atp6, nad1 (partial) and cox1 (partial) genes from E. multilocularis, E. vogeli, E. oligarthrus, 5 genotypes of E. granulosus, and using T. crassiceps as an outgroup, yielded the same results. We conclude that the sheep-dog and horse-dog strains of E. granulosus should be regarded as distinct at the specific level.

Leaf area development of ABA-deficient and wild-type peas at two levels of nitrogen supply

The ABA-deficient wilty pea (Pisum sativum L.) and its wild-type (WT) were grown at two levels of nitrogen supply (0.5 and 5.0 mM) for 5-6 weeks from sowing, to determine whether leaf ABA status altered the leaf growth response to N deprivation. Plants were grown at high relative humidity to prevent wilting of the wilty peas. Irrespective of N supply, expanding wilty leaflets had ca 50% less ABA than WT leaflets but similar ethylene evolution rates. Fully expanded wilty leaflets had lower relative water contents (RWC) and were 10-60% smaller in area (according to the node of measurement) than WT leaflets. However, there were no genotypic differences in plant relative leaf expansion rate (RLER). Growth of both genotypes at 0.5 mM N increased the RWC of fully expanded leaflets, but did not alter ethylene evolution or ABA concentration of expanding leaflets. Plants grown at 0.5 mM N showed a 20-30% reduction in RLER, which was similar in magnitude in both wilty and WT peas. Thus, leaf ABA status did not alter the leaf growth response to N deprivation.

Osteoarthritis of the hands, hips and knees in an Australian twin sample - Evidence of association with the aggrecan VNTR polymorphism

Age-related changes in the composition of the cartilage matrix may be associated with the development of osteoarthritis, a relatively late-onset disease characterised by the destruction of joint cartilage. In order to investigate whether differences in the VNTR polymorphic region of aggrecan affect cartilage functionality and therefore the development of osteoarthritis, we examined the aggrecan polymorphic genotypes of a sample of 134 Australian twins aged over 50 (including 34 monozygotic and 27 dizygotic twin pairs). Clinical measures of hand, hip and knee osteoarthritis, as well as self-reported bone and joint pain, were tested for association with the aggrecan polymorphism. The results were consistent with either a deleterious effect of allele 27, or a protective effect of alleles 25 and 28, providing some additional evidence for an association between the aggrecan VNTR polymorphism and osteoarthritis of the hands, hips and knees.

Pathogenic specialisation within Colletotrichum trifolii in Australia, and lucerne cultivar reactions to all known Australian pathotypes

Anthracnose and crown rot, caused by Colletotrichum trifolii, are serious diseases of lucerne (Medicago saliva L.) in humid regions of the world. A race survey was conducted by inoculating individual lucerne clones (genotypes) with C. trifolii isolates collected from a range of Medicago hosts, locations, and years in south-eastern Queensland. This survey revealed for the first time in Australia the presence of race 2 (virulence on anthracnose resistance gene An I) and the first world report of race 4 (virulence on An(2)). A collection of North American race I and race 2 C. trifolii isolates, when inoculated onto the Australian differential clones, gave responses that were in agreement with their North American reactions. A RAPD analysis was conducted on 9 Australian C. trifolii isolates including races 1, 2, and 4; two C. destructivum and one C. gloeosporioides isolate were included as known outliers. For the C. trifolii isolates, 94.6% similarity was found regardless of host origin or race, compared with 2.2% similarity between this group and the C. gloeosporioides and C. destructivum isolates, confirming that the new races belong to C. trifolii. Currently, it is hypothesised that only plants carrying genes An, and An2 are resistant to the 3 races. Of 22 cultivars screened against the 3 races, only UQL-1, Hallmark, and Pioneer 54Q53 had >30% of plants resistant to the 3 races in separate screenings. The research highlights the need to find new sources of resistance to C. trifolii in lucerne.

eMelanoBase: An online locus-specific variant database for familial melanoma

A proportion of melanoma,prone individuals in both familial and non,familial contexts has been shown to carry inactivating mutations in either CDKN2A or, rarely, CDK4. CDKN2A is a complex locus that encodes two unrelated proteins from alternately spliced transcripts that are read in different frames. The alpha transcript (exons 1a, 2, and 3) produces the p16INK4A cyclin-dependent kinase inhibitor, while the beta transcript (exons 1beta and 2) is translated as p14ARF, a stabilizing factor of p53 levels through binding to MDM2. Mutations in exon 2 can impair both polypeptides and insertions and deletions in exons 1alpha, 1beta, and 2, which can theoretically generate p16INK4A,p14ARF fusion proteins. No online database currently takes into account all the consequences of these genotypes, a situation compounded by some problematic previous annotations of CDKN2A related sequences and descriptions of their mutations. As an initiative of the international Melanoma Genetics Consortium, we have therefore established a database of germline variants observed in all loci implicated in familial melanoma susceptibility. Such a comprehensive, publicly accessible database is an essential foundation for research on melanoma susceptibility and its clinical application. Our database serves two types of data as defined by HUGO. The core dataset includes the nucleotide variants on the genomic and transcript levels, amino acid variants, and citation. The ancillary dataset includes keyword description of events at the transcription and translation levels and epidemiological data. The application that handles users' queries was designed in the model,view. controller architecture and was implemented in Java. The object-relational database schema was deduced using functional dependency analysis. We hereby present our first functional prototype of eMelanoBase. The service is accessible via the URL www.wmi.usyd.e, du.au:8080/melanoma.html.

Dizygotic twinning is not associated with methylenetetrahydrofolate reductase haplotypes

Background: Folate metabolism is critical to embryonic development, influencing neural tube defects (NTD) and recurrent early pregnancy loss. Polymorphisms in 5,10-methylenetetrahydrofolate reductase (MTHFR) have been associated with dizygotic (DZ) twinning through pregnancy loss. Methods: The C677T and A1298C polymorphisms in MTHFR were genotyped in 258 Australasian families (1016 individuals) and 118 Dutch families (462 individuals) of mothers of DZ twins and a population sample of 462 adolescent twin families (1861 individuals). Haplotypes were constructed from the alleles, and transmission of the MTHFR haplotypes to mothers of DZ twins and from parents to twins in the adolescent twin families analysed. Results: The C677T and A1298C were common in all three populations (frequencies > 0.29). There was strong linkage disequilibrium (D'=1) between the variants, showing that specific combinations of alleles (haplotypes) were transmitted together. Three haplotypes accounted for nearly all the variation. There was no evidence of any association between MTHFR genotype and twinning in mothers of twins, or of the loss of specific MTHFR genotypes during twin pregnancies. Conclusions: It is concluded that variation in twinning frequency is not associated with MTHFR genotype.

ADH genotype does not modify the effects of alcohol on high-density lipoprotein

Background: Alcohol consumption has beneficial effects on mortality which are mainly due to reduction in cardiovascular disease. These are believed to be due, at least in part, to the increase in plasma high-density lipoprotein (HDL) which is associated with alcohol consumption. It has been proposed that ADH3 genotype modifies the relationships between alcohol intake and cardiovascular disease by altering the HDL response to alcohol. The aim of this paper was to test for effects of ADH2 and ADH3 genotypes on the response of HDL components to habitual alcohol consumption. Methods: Adult male and female subjects were genotyped for ADH2 and ADH3; and plasma HDL cholesterol, apolipoprotein A-I, and apolipoprotein A-II were measured. Nine hundred one subjects had both ADH2 and ADH3 genotypes and HDL cholesterol results, while 753 had both genotypes and all three lipid results. The effect of alcohol intake on the three measured HDL components, and a factor score derived from them, was estimated for each of the ADH2 and ADH3 genotype groups. Results: All the measured components of HDL increased with increasing alcohol consumption over the range of intakes studied, 0-4 drinks per day. There were no significant interactions between alcohol consumption and ADH2 or ADH3 genotypes. Conclusions: The concept that alcohol dehydrogenase genotype and alcohol metabolic rate modify the effects of alcohol on plasma HDL concentration is not supported by our results.

Screening of human primary melanocytes of defined melanocortin-1 receptor genotype: Pigmentation marker, ultrastructural and UV-survival studies

Recent population studies have demonstrated an association with the red-hair and fair-skin phenotype with variant alleles of the melanocortin-1 receptor (MC1R) which result in amino acid substitutions within the coding region leading to an altered receptor activity. In particular, Arg151Cys, Arg160Trp and Asp294His were the most commonly associated variants seen in the south-east Queensland population with at least one of these alleles found in 93% of those with red hair. In order to study the individual effects of these variants on melanocyte biology and melanocytic pigmentation, we established a series of human melanocyte strains genotyped for the MC1R receptor which included wild-type consensus, variant heterozygotes, compound heterozygotes and homozygotes for Arg151Cys, Arg160Trp, Val60Leu and Val92Met alleles. These strains ranged from darkly pigmented to amelanotic, with all strains of consensus sequence having dark pigmentation. UV sensitivity was found not to be associated with either MC1R genotype or the level of pigmentation with a range of sensitivities seen across all genotypes. Ultrastructural analysis demonstrated that while consensus strains contained stage IV melanosomes in their terminal dendrites, Arg151Cys and Arg160Trp homozygote strains contained only stage II melanosomes. This was despite being able to show expression of tyrosinase and tyrosinase-related protein-1 markers, although at reduced levels and an ability to convert exogenous 3,4-dihydroxyphenyl-alanine (DOPA) to melanin in these strains.

Complete genomic sequence of the Australian south-west genotype of Sindbis virus: Comparisons with other Sindbis strains and identification of a unique deletion in the 3 '-untranslated region

Our previous studies have shown that two distinct genotypes of Sindbis (SIN) virus occur in Australia. One of these, the Oriental/Australian type, circulates throughout most of the Australian continent, whereas the recently identified south-west (SW) genetic type appears to be restricted to a distinct geographic region located in the temperate south-west of Australia. We have now determined the complete nucleotide and translated amino acid sequences of a SW isolate of SIN virus (SW6562) and performed comparative analyses with other SIN viruses at the genomic level. The genome of SW6562 is 11,569 nucleotides in length, excluding the cap nucleotide and poly (A) tail. Overall this virus differs from the prototype SIN virus (strain AR339) by 23% in nucleotide sequence and 12.5% in amino acid sequence. Partial sequences of four regions of the genome of four SW isolates were determined and compared with the corresponding sequences from a number of SIN isolates from different regions of the World. These regions are the non-structural protein (nsP3), the E2 gene, the capsid gene, and the repeated sequence elements (RSE) of the 3'UTR. These comparisons revealed that the SW SIN viruses were more closely related to South African and European strains than to other Australian isolates of SIN virus. Thus the SW genotype of SIN virus may have been introduced into this region of Australia by viremic humans or migratory birds and subsequently evolved independently in the region. The sequence data also revealed that the SW genotype contains a unique deletion in the RSE of the 3'UTR region of the genome. Previous studies have shown that deletions in this region of the SIN genome can have significant effects on virus replication in mosquito and avian cells, which may explain the restricted distribution of this genotype of SIN virus.

Codon usage bias and A+T content variation in human papillomavirus genomes

We analyzed the codon usage bias of eight open reading frames (ORFs) across up to 79 human papillomavirus (HPV) genotypes from three distinct phylogenetic groups. All eight ORFs across HPV genotypes show a strong codon usage bias, amongst degenerately encoded amino acids, toward 18 codons mainly with T at the 3rd position. For all 18 degenerately encoded amino acids, codon preferences amongst human and animal PV ORFs are significantly different from those averaged across mammalian genes. Across the HPV types, the L2 ORFs show the highest codon usage bias (73.2 +/- 1.6% and the E4 ORFs the lowest (51.1 +/- 0.5%), reflecting as similar bias in codon 3rd position A + T content (L2: 76.1 +/- 4.2%; E4: 58.6 +/- 4.5%). The E4 ORF, uniquely amongst the HPV ORFs, is G + C rich, while the other ORFs are A + T rich. Codon usage bias correlates positively with A + T content at the codon 3rd position in the E2, E6, L1 and L2 ORFs, but negatively in the E4 ORFs. A general conservation of preferred codon usage across human and non-human PV genotypes whether they originate from a same supergroup or not, together with observed difference between the preferred codon usage for HPV ORFs and for genes of the cells they infect, suggests that specific codon usage bias and A + T content variation may somehow increase the replicational fitness of HPVs in mammalian epithelial cells, and have practical implications for gene therapy of HPV infection. (C) 2003 Elsevier B.V. All rights reserved.

Part 1. Growth, carcass and meat quality parameters of male goats: effects of genotype and liveweight at slaughter

Male kids (110) from six goat genotypes, i.e. Boer x Angora (BA), Boer x Feral (1317), Boer x Saanen (BS), Feral x Feral (FF), Saanen x Angora (SA) and Saanen x Feral (SF) and two slaughter weight groups, i.e. Capretto and Chevon (liveweight at slaughter 14-22 and 30-35 kg, respectively) were compared for growth, carcass and meat quality characteristics. Due to their better growth rate, kids from BS and SF genotypes reached the required liveweight for slaughter earlier than kids from other Genotypes used in the study. Chevon kids had a significantly (P < 0.05) lower average daily gain (119 g per day) compared to Capretto kids (171 g per day). SA, SF and FF kids deposited more internal fat in comparison to kids from other genotypes. The dressing percentage of kids ranged from 51 to 54%, with significant differences between genotypes. BS and SF kids had longer carcasses. while BF kids had larger eye muscle area compared to other genotypes. Goat carcasses had a thin subcutaneous fat cover (1.6-2.2 mm). Genotype had a significant (P < 0.05) influence on cooking loss, pigment concentration and muscle colour parameters (CIE L*, a* and b* values). As denoted by the higher V and fibre optic probe values and lower subjective muscle score, the longissimus muscle colour was lighter for BS kids than other genotypes. Cooked meat from the BF kids had lower shear force values and better sensory scores compared to other genotypes. A significant (P < 0.05) decrease in muscle tenderness was observed from Capretto to Chevon carcasses, whereas cooked meat from these two slaughter weight groups was equally accepted (P > 0.05) by the panellists. (C) 2003 Elsevier Science B.V. All rights reserved.

Part 2. Carcass composition and fatty acid profiles of adipose tissue of male goats: effects of genotype and liveweight at slaughter

The dissected carcass composition and fatty acid profiles of intermuscular fat from 110 male goat kids from six genotypes i.e. Boer x Angora (BA), Boer x Feral (BF), Boer x Saanen (BS), Feral x Feral (1717), Saanen x Angora (SA) and Saanen x Feral (SF) and two slaughter weight groups i.e. Capretto and Chevon (liveweight at slaughter 14-22 and 30-35 kg, respectively) were compared. Carcass tissue distribution for various genotypes was: muscle (63-66%), fat (10-13%) and bone (21-24%). Genotype significantly (P < 0.05) influenced the carcass composition; BA and FF carcasses had significantly higher muscle to bone ratio, while carcasses from BS kids were leaner compared to other genotypes. However, the two slaughter weight groups did not differ significantly (P > 0.05) in terms of carcass composition, when compared at the same carcass weight. In the present study, significant (P < 0.01) correlations were observed between percentage of muscle, fat and bone in most of the primal cuts and that in the carcass side. The main saturated fatty acids (SFAs) identified were palmitic (16:0) and stearic acid (18:0), while oleic acid (18: 1, omega9) was the main unsaturated fatty acid (UFA) in the intermuscular fat from goat kids. There were significant (P < 0.05) differences between genotypes in the proportions of individual fatty acids. Adipose tissue from BS kids had significantly higher UFAs (mainly oleic acid) and thus had a significantly lower melting point compared to other genotypes. There were significantly higher proportions of palmitic acid (35%) in the adipose tissue from Capretto kids compared to that from Chevon kids (22%). The concentration of UFAs increased in the adipose tissue from Capretto to Chevon carcasses. (C) 2003 Elsevier Science B.V. All rights reserved.