A linear programming approach for estimating the structure of a sparse linear genetic network from transcript profiling data

Background: A genetic network can be represented as a directed graph in which a node corresponds to a gene and a directed edge specifies the direction of influence of one gene on another. The reconstruction of such networks from transcript profiling data remains an important yet challenging endeavor. A transcript profile specifies the abundances of many genes in a biological sample of interest. Prevailing strategies for learning the structure of a genetic network from high-dimensional transcript profiling data assume sparsity and linearity. Many methods consider relatively small directed graphs, inferring graphs with up to a few hundred nodes. This work examines large undirected graphs representations of genetic networks, graphs with many thousands of nodes where an undirected edge between two nodes does not indicate the direction of influence, and the problem of estimating the structure of such a sparse linear genetic network (SLGN) from transcript profiling data. Results: The structure learning task is cast as a sparse linear regression problem which is then posed as a LASSO (l1-constrained fitting) problem and solved finally by formulating a Linear Program (LP). A bound on the Generalization Error of this approach is given in terms of the Leave-One-Out Error. The accuracy and utility of LP-SLGNs is assessed quantitatively and qualitatively using simulated and real data. The Dialogue for Reverse Engineering Assessments and Methods (DREAM) initiative provides gold standard data sets and evaluation metrics that enable and facilitate the comparison of algorithms for deducing the structure of networks. The structures of LP-SLGNs estimated from the INSILICO1, INSILICO2 and INSILICO3 simulated DREAM2 data sets are comparable to those proposed by the first and/or second ranked teams in the DREAM2 competition. The structures of LP-SLGNs estimated from two published Saccharomyces cerevisae cell cycle transcript profiling data sets capture known regulatory associations. In each S. cerevisiae LP-SLGN, the number of nodes with a particular degree follows an approximate power law suggesting that its degree distributions is similar to that observed in real-world networks. Inspection of these LP-SLGNs suggests biological hypotheses amenable to experimental verification. Conclusion: A statistically robust and computationally efficient LP-based method for estimating the topology of a large sparse undirected graph from high-dimensional data yields representations of genetic networks that are biologically plausible and useful abstractions of the structures of real genetic networks. Analysis of the statistical and topological properties of learned LP-SLGNs may have practical value; for example, genes with high random walk betweenness, a measure of the centrality of a node in a graph, are good candidates for intervention studies and hence integrated computational – experimental investigations designed to infer more realistic and sophisticated probabilistic directed graphical model representations of genetic networks. The LP-based solutions of the sparse linear regression problem described here may provide a method for learning the structure of transcription factor networks from transcript profiling and transcription factor binding motif data.

A large-scale analysis of genetic variants within putative miRNA binding sites in prostate cancer

Prostate cancer is the second most common malignancy among men worldwide. Genome-wide association studies have identified 100 risk variants for prostate cancer, which can explain approximately 33% of the familial risk of the disease. We hypothesized that a comprehensive analysis of genetic variations found within the 3' untranslated region of genes predicted to affect miRNA binding (miRSNP) can identify additional prostate cancer risk variants. We investigated the association between 2,169 miRSNPs and prostate cancer risk in a large-scale analysis of 22,301 cases and 22,320 controls of European ancestry from 23 participating studies. Twenty-two miRSNPs were associated (P<2.3×10(-5)) with risk of prostate cancer, 10 of which were within 7 genes previously not mapped by GWAS studies. Further, using miRNA mimics and reporter gene assays, we showed that miR-3162-5p has specific affinity for the KLK3 rs1058205 miRSNP T-allele, whereas miR-370 has greater affinity for the VAMP8 rs1010 miRSNP A-allele, validating their functional role. SIGNIFICANCE Findings from this large association study suggest that a focus on miRSNPs, including functional evaluation, can identify candidate risk loci below currently accepted statistical levels of genome-wide significance. Studies of miRNAs and their interactions with SNPs could provide further insights into the mechanisms of prostate cancer risk.

A genetic variant of MDM4 influences regulation by multiple microRNAs in prostate cancer

The oncogene MDM4, also known as MDMX or HDMX, contributes to cancer susceptibility and progression through its capacity to negatively regulate a range of genes with tumour-suppressive functions. As part of a recent genome-wide association study it was determined that the A-allele of the rs4245739 SNP (A>C), located in the 3'-UTR of MDM4, is associated with an increased risk of prostate cancer. Computational predictions revealed that the rs4245739 SNP is located within a predicted binding site for three microRNAs (miRNAs): miR-191-5p, miR-887 and miR-3669. Herein, we show using reporter gene assays and endogenous MDM4 expression analyses that miR-191-5p and miR-887 have a specific affinity for the rs4245739 SNP C-allele in prostate cancer. These miRNAs do not affect MDM4 mRNA levels, rather they inhibit its translation in C-allele-containing PC3 cells but not in LNCaP cells homozygous for the A-allele. By analysing gene expression datasets from patient cohorts, we found that MDM4 is associated with metastasis and prostate cancer progression and that targeting this gene with miR-191-5p or miR-887 decreases in PC3 cell viability. This study is the first, to our knowledge, to demonstrate regulation of the MDM4 rs4245739 SNP C-allele by two miRNAs in prostate cancer, and thereby to identify a mechanism by which the MDM4 rs4245739 SNP A-allele may be associated with an increased risk for prostate cancer.

Genetic Options for Drought Management in Groundnut

Groundnut is one of the principal oilseeds in the world. It is cultivated on 24.8 million ha with a total production of 32.8 million t and an average productivity of 1.32 t ha-'. Developing countries account for 96.9% of the world groundnut area and 93.8% of total production. Production is concentrated in Asia (56.8% area and 66.5% production of the world) and Africa (38.0% area and 24.7% production). The groundnut productivity in Africa is only 0.86 t ha-' compared with 1.55 t hx1 of Asia. The world groundnut economy-facts, trends and outlook are desaibed in detail by Freeman et al., 1999. Briefly, in medium-term (i.e. up to 2010), 'groundnut production and consumption is likely to shift increasingly to developing countries; production will grow in all regions but most rapidly in Asia, slowly in sub-Saharan Africa and decline in Latin America; and utilizationwill continue to shift away from groundnut oil toward groundnut meal, specially confectionery products'.

A robust approach to design a single facility layout plan in dynamic manufacturing environments using a permutation-based genetic algorithm

During the past few decades, developing efficient methods to solve dynamic facility layout problems has been focused on significantly by practitioners and researchers. More specifically meta-heuristic algorithms, especially genetic algorithm, have been proven to be increasingly helpful to generate sub-optimal solutions for large-scale dynamic facility layout problems. Nevertheless, the uncertainty of the manufacturing factors in addition to the scale of the layout problem calls for a mixed genetic algorithm–robust approach that could provide a single unlimited layout design. The present research aims to devise a customized permutation-based robust genetic algorithm in dynamic manufacturing environments that is expected to be generating a unique robust layout for all the manufacturing periods. The numerical outcomes of the proposed robust genetic algorithm indicate significant cost improvements compared to the conventional genetic algorithm methods and a selective number of other heuristic and meta-heuristic techniques.

Genetic analysis of resistance to root-lesion nematode (Pratylenchus thornei) in wheat

The inheritance of resistance to root-lesion nematode was investigated in five synthetic hexaploid wheat lines and two bread wheat lines using a half-diallel design of F1 and F2 crosses. The combining ability of resistance genes in the synthetic hexaploid wheat lines was compared with the performance of the bread wheat line 'GS50a', the source of resistance to Pratylenchus thornei used in Australian wheat breeding programmes. Replicated glasshouse trials identified P. thornei resistance as polygenic and additive in gene action. General combining ability (GCA) of the parents was more important than specific combining ability (SCA) effects in the inheritance of P. thornei resistance in both F1 and F2 populations. The synthetic hexaploid wheat line 'CPI133872' was identified as the best general combiner, however, all five synthetic hexaploid wheat lines possessed better GCA than 'GS50a'. The synthetic hexaploid wheat lines contain novel sources of P. thornei resistance that will provide alternative and more effective sources of resistance to be utilized in wheat breeding programmes

Pronounced genetic population structure in a potentially vagile fish species (Pristipomoides multidens, Teleostei; Perciformes; Lutjanidae) from the East Indies triangle

The East Indies triangle, bordered by the Phillipines, Malay Peninsula and New Guinea, has a high level of tropical marine species biodiversity. Pristipomoides multidens is a large, long-lived, fecund snapper species that is distributed throughout the East Indies and Indo-Pacific. Samples were analysed from central and eastern Indonesia and northern Australia to test for genetic discontinuities in population structure. Fish (n = 377) were collected from the Indonesian islands of Bali, Sumbawa, Flores, West Timor, Tanimbar and Tual along with 131 fish from two northern Australian locations (Arafura and Timor Seas) from a previous study. Genetic variation in the control region of the mitochondrial genome was assayed using restriction fragment length polymorphism and direct sequencing. Haplotype diversity was high (0.67-0.82), as was intraspecific sequence divergence (range 0-5.8%). FST between pairs of populations ranged from 0 to 0.2753. Genetic subdivision was apparent on a small spatial scale; FST was 0.16 over 191 km (Bali/Sumbawa) and 0.17 over 491 km (Bali/Flores). Constraints to dispersal that contribute to, and maintain, the observed degree of genetic subdivision are experienced presumably by all life history stages of this tropical marine finfish. The constraints may include (1) little or no movement of eggs or larvae, (2) little or no home range or migratory movement of adults and (3) loss of larval cohorts due to transport of larvae away from suitable habitat by prevailing currents

Genetic response to growth rate selection in pigs measured on two feeding levels

From the findings of McPhee et al. (1988), there is an expectation that selection in the growing pig for bodyweight gain measured on restricted feeding will result in favourable responses in the rate and efficiency of growth of lean pork on different levels of feeding. This paper examines this in two lines of Australian Large White pigs which have undergone 3 years of selection for high and for low growth rate over a 6-week period starting at 50 kg liveweight. Over this test period, pigs of both lines are all fed the same total amount of grower food, restricted to an estimated 80% of average ad libitum intake. 'Animal production for a consuming world': proceedings of 9th Congress of the AAAAP Societies and 23rd Biennial Conference of the ASAP and 17th Annual Symposium of the University of Sydney, Dairy Research Foundation, (DRF). Sydney, Australia.

Functionally associated molecular genetic marker map construction in perennial ryegrass (Lolium perenne L.)

A molecular marker-based map of perennial ryegrass (Lolium perenne L.) has been constructed through the use of polymorphisms associated with expressed sequence tags (ESTs). A pair-cross between genotypes from a North African ecotype and the cultivar Aurora was used to generate a two-way pseudo-testcross population. A selection of 157 cDNAs assigned to eight different functional categories associated with agronomically important biological processes was used to detect polymorphic EST–RFLP loci in the F1(NA6 × AU6) population. A comprehensive set of EST–SSR markers was developed from the analysis of 14,767 unigenes, with 310 primer pairs showing efficient amplification and detecting 113 polymorphic loci. Two parental genetic maps were produced: the NA6 genetic map contains 88 EST–RFLP and 71 EST–SSR loci with a total map length of 963 cM, while the AU6 genetic map contains 67 EST–RFLP and 58 EST–SSR loci with a total map length of 757 cM. Bridging loci permitted the alignment of homologous chromosomes between the parental maps, and a sub-set of genomic DNA-derived SSRs was used to relate linkage groups to the perennial ryegrass reference map. Regions of segregation distortion were identified, in some instances in common with other perennial ryegrass maps. The EST-derived marker-based map provides the basis for in silico comparative genetic mapping, as well as the evaluation of co-location between QTLs and functionally associated genetic loci.

Timing the Emergence of Resistance to Anti-HIV Drugs with Large Genetic Barriers

New antiretroviral drugs that offer large genetic barriers to resistance, such as the recently approved inhibitors of HIV-1 protease, tipranavir and darunavir, present promising weapons to avert the failure of current therapies for HIV infection. Optimal treatment strategies with the new drugs, however, are yet to be established. A key limitation is the poor understanding of the process by which HIV surmounts large genetic barriers to resistance. Extant models of HIV dynamics are predicated on the predominance of deterministic forces underlying the emergence of resistant genomes. In contrast, stochastic forces may dominate, especially when the genetic barrier is large, and delay the emergence of resistant genomes. We develop a mathematical model of HIV dynamics under the influence of an antiretroviral drug to predict the waiting time for the emergence of genomes that carry the requisite mutations to overcome the genetic barrier of the drug. We apply our model to describe the development of resistance to tipranavir in in vitro serial passage experiments. Model predictions of the times of emergence of different mutant genomes with increasing resistance to tipranavir are in quantitative agreement with experiments, indicating that our model captures the dynamics of the development of resistance to antiretroviral drugs accurately. Further, model predictions provide insights into the influence of underlying evolutionary processes such as recombination on the development of resistance, and suggest guidelines for drug design: drugs that offer large genetic barriers to resistance with resistance sites tightly localized on the viral genome and exhibiting positive epistatic interactions maximally inhibit the emergence of resistant genomes.

Structural Basis of Drug Resistance by Genetic Variants of HIV Type 1 Clade C Protease from India

Using computer modeling of three-dimensional structures and structural information available on the crystal structures of HIV-1 protease, we investigated the structural effects of mutations, in treatment-naive and treatment-exposed individuals from India and postulated mechanisms of resistance in clade C variants. A large number of models (14) have been generated by computational mutation of the available crystal structures of drug bound proteases. Localized energy minimization was carried out in and around the sites of mutation in order to optimize the geometry of interactions present. Most of the mutations result in structural differences at the flap that favors the semiopen state of the enzyme. Some of the mutations were also found to confer resistance by affecting the geometry of the active site. The E35D mutation affects the flap structure in clade B strains and E35N and E35K mutation, seen in our modeled strains, have a more profound effect. Common polymorphisms at positions 36 and 63 in clade C also affected flap structure. Apart from a few other residues Gln-58, Asn-83, Asn-88, and Gln-92 and their interactions are important for the transition from the closed to the open state. Development of protease inhibitors by structure-based design requires investigation of mechanisms operative for clade C to improve the efficacy of therapy.

Genetic Parameters for Production Traits in Queensland Merinos

Genetic and phenotypic parameters are presented for production traits, greasy fleece weight (GFW), yield (YLD), clean fleece weight (CFW), average fibre diameter (DIAM) and liveweight (LWT), in 15 month old medium Peppin Merino sheep at Longreach and Julia Creek, Queensland. Heritabilities for GFW, YLD, CFW, DIAM and LWT were respectively 0.35, 0.62, 0.34, 0.74, and 0.37 for Longreach and 0.23, 0.52, 0.20, 0.67 and 0.56 for Julia Creek. Most estimates were consistent with other reported values. AAABG 13th Conference; Proceedings of the Association for the Advancement of Animal Breeding and Genetics.

Genetic association analysis of miRNA SNPs implicates MIR145 in breast cancer susceptibility

Background MicroRNAs (miRNAs) are important small non-coding RNA molecules that regulate gene expression in cellular processes related to the pathogenesis of cancer. Genetic variation in miRNA genes could impact their synthesis and cellular effects and single nucleotide polymorphisms (SNPs) are one example of genetic variants studied in relation to breast cancer. Studies aimed at identifying miRNA SNPs (miR-SNPs) associated with breast malignancies could lead towards further understanding of the disease and to develop clinical applications for early diagnosis and treatment. Methods We genotyped a panel of 24 miR-SNPs using multiplex PCR and chip-based matrix assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry (MS) analysis in two Caucasian breast cancer case control populations (Primary population: 173 cases and 187 controls and secondary population: 679 cases and 301 controls). Association to breast cancer susceptibility was determined using chi-square (X 2 ) and odds ratio (OR) analysis. Results Statistical analysis showed six miR-SNPs to be non-polymorphic and twelve of our selected miR-SNPs to have no association with breast cancer risk. However, we were able to show association between rs353291 (located in MIR145) and the risk of developing breast cancer in two independent case control cohorts (p = 0.041 and p = 0.023). Conclusions Our study is the first to report an association between a miR-SNP in MIR145 and breast cancer risk in individuals of Caucasian background. This finding requires further validation through genotyping of larger cohorts or in individuals of different ethnicities to determine the potential significance of this finding as well as studies aimed to determine functional significance. Keywords: Association analysis; Breast cancer; microRNA; miR-SNPs; MIR145

Genetic variability of Tomato spotted wilt virus in Australia and validation of real time RT-PCR for its detection in single and bulked leaf samples

The potential for large-scale use of a sensitive real time reverse transcription polymerase chain reaction (RT-PCR) assay was evaluated for the detection of Tomato spotted wilt virus (TSWV) in single and bulked leaf samples by comparing its sensitivity with that of DAS-ELISA. Using total RNA extracted with RNeasy® or leaf soak methods, real time RT-PCR detected TSWV in all infected samples collected from 16 horticultural crop species (including flowers, herbs and vegetables), two arable crop species, and four weed species by both assays. In samples in which DAS-ELISA had previously detected TSWV, real time RT-PCR was effective at detecting it in leaf tissues of all 22 plant species tested at a wide range of concentrations. Bulk samples required more robust and extensive extraction methods with real time RT-PCR, but it generally detected one infected sample in 1000 uninfected ones. By contrast, ELISA was less sensitive when used to test bulked samples, once detecting up to 1 infected in 800 samples with pepper but never detecting more than 1 infected in 200 samples in tomato and lettuce. It was also less reliable than real time RT-PCR when used to test samples from parts of the leaf where the virus concentration was low. The genetic variability among Australian isolates of TSWV was small. Direct sequencing of a 587 bp region of the nucleoprotein gene (S RNA) of 29 isolates from diverse crops and geographical locations yielded a maximum of only 4.3% nucleotide sequence difference. Phylogenetic analysis revealed no obvious groupings of isolates according to geographic origin or host species. TSWV isolates, that break TSWV resistance genes in tomato or pepper did not differ significantly in the N gene region studied, indicating that a different region of the virus genome is responsible for this trait.

Rapid and efficient screening of a Representational Difference Analysis library using reverse Southern hybridisation: Identification of genetic differences between Haemophilus parasuis isolates

Representational Difference Analysis (RDA) is an established technique used for isolation of specific genetic differences between or within bacterial species. This method was used to investigate the genetic basis of serovar-specificity and the relationship between serovar and virulence in Haemophilus parasuis. An RDA clone library of 96 isolates was constructed using H. parasuis strains H425(P) (serovar 12) and HS1967 (serovar 4). To screen such a large clone library to determine which clones are strain-specific would typically involved separately labelling each clone for use in Southern hybridisation against genomic DNA from each of the strains. In this study, a novel application of reverse Southern hybridisation was used to screen the RDA library: genomic DNA from each strain was labelled and used to probe the library to identify strain-specific clones. This novel approach represents a significant improvement in methodology that is rapid and efficient.