Weighted gene coexpression network analysis strategies applied to mouse weight

Systems-oriented genetic approaches that incorporate gene expression and genotype data are valuable in the quest for genetic regulatory loci underlying complex traits. Gene coexpression network analysis lends itself to identification of entire groups of differentially regulated genes—a highly relevant endeavor in finding the underpinnings of complex traits that are, by definition, polygenic in nature. Here we describe one such approach based on liver gene expression and genotype data from an F₂ mouse intercross utilizing weighted gene coexpression network analysis (WGCNA) of gene expression data to identify physiologically relevant modules. We describe two strategies: single-network analysis and differential network analysis. Single-network analysis reveals the presence of a physiologically interesting module that can be found in two distinct mouse crosses. Module quantitative trait loci (mQTLs) that perturb this module were discovered. In addition, we report a list of genetic drivers for this module. Differential network analysis reveals differences in connectivity and module structure between two networks based on the liver expression data of lean and obese mice. Functional annotation of these genes suggests a biological pathway involving epidermal growth factor (EGF). Our results demonstrate the utility of WGCNA in identifying genetic drivers and in finding genetic pathways represented by gene modules. These examples provide evidence that integration of network properties may well help chart the path across the gene–trait chasm. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users. Tova F. Fuller, Anatole Ghazalpour contributed equally to this work.     

EDISA: extracting biclusters from multiple time-series of gene expression profiles

Background  

Cells dynamically adapt their gene expression patterns in response to various stimuli. This response is orchestrated into a number of gene expression modules consisting of co-regulated genes. A growing pool of publicly available microarray datasets allows the identification of modules by monitoring expression changes over time. These time-series datasets can be searched for gene expression modules by one of the many clustering methods published to date. For an integrative analysis, several time-series datasets can be joined into a three-dimensional gene-condition-time dataset, to which standard clustering or biclustering methods are, however, not applicable. We thus devise a probabilistic clustering algorithm for gene-condition-time datasets.     

Mobilization of pdif modules in Acinetobacter: A novel mechanism for antibiotic resistance gene shuffling?

XerCD-dif site-specific recombination is a well characterized system, found in most bacteria and archaea. Its role is resolution of chromosomal dimers that arise from homologous recombination. Xer-mediated recombination is also used by several plasmids for multimer resolution to enhance stability and by some phage for integration into the chromosome. In the past decade, it has been hypothesized that an alternate and novel function exists for this system in the dissemination of genetic elements, notably antibiotic resistance genes, in Acinetobacter species. Currently the mechanism underlying this apparent genetic mobility is unknown. Multidrug resistant Acinetobacter baumannii is an increasingly problematic pathogen that can cause recurring infections. Sequencing of numerous plasmids from clinical isolates of A. baumannii revealed the presence of possible mobile modules: genes were found flanked by pairs of Xer recombination sites, called plasmid-dif (pdif) sites. These modules have been identified in multiple otherwise unrelated plasmids and in different genetic contexts suggesting they are mobile elements. In most cases, the pairs of sites flanking a gene (or genes) are in inverted repeat, but there can be multiple modules per plasmid providing pairs of recombination sites that can be used for inversion or fusion/deletion reactions; as many as 16 pdif sites have been seen in a single plasmid. Similar modules including genes for surviving environmental toxins have also been found in strains of Acinetobacter Iwoffi isolated from permafrost cores; this suggests that these mobile modules are an ancient adaptation and not a novel response to antibiotic pressure. These modules bear all the hallmarks of mobile genetic elements, yet, their movement has never been directly observed to date. This review gives an overview of the current state of this novel research field.     

Expression of inositol 1,3,4-trisphosphate 5/6-kinase (ITPK1) and its role in neural tube defects

ITPK1 is the rate-limiting enzyme in the pathway leading to formation of the highly phosphorylated inositol phosphates including IP₆ and the inositol pyrophosphates. One or more of these metabolites are essential for life as deletion of either of the kinases that form IP₅ or IP₆ in mice results in embryonic lethality. We have produced mice harboring a hypomorphic allele for Itpk1, and mice homozygous for this gene trap allele produce low but detectable levels of active enzyme. We have studied the expression of Itpk1 in various tissues and found that the enzyme is highly expressed in smooth muscle of vessels and other tissues. In addition, these mice have neural tube defects in 12% of homozygous embryos. Since the levels of enzyme expression vary greatly in homozygous animals, we speculate that relative deficiency of one or more inositol phosphates accounts for these defects. We plan to feed an inositol deficient diet or one with supplemental inositol to animals to demonstrate altered prevalence of neural tube defects.     

Mapping of genetic modifiers affecting the eye phenotype of ocular retardation (Chx10 or-J ) mice

Ocular retardation is a recessive murine mutation whose phenotypic expression is greatly affected by genetic background effects. Mice of the inbred 129/SvJ background that are homozygous for the Chx10^or-J mutation are blind and have a thin, poorly differentiated retina and no optic nerve. A backcross between 129/SvJ and Mus musculus castaneus (CASA/Rk) produced animals that were homozygous for the Chx10^or-J mutation, yet showed a much milder phenotype. Such animals, when brother-sister mated and selected for mild phenotype for several generations, resulted in partial recovery of visual function, including presence of an optic nerve and pupillary response. In this article we report a genome scan of phenotypic extremes of the backcross to identify the genetic loci affecting this phenotype modification. Our scan revealed significant loci on Chromosomes 6 and 14 where the CASA/Rk alleles are maintained selectively. Markers were developed near candidate genes, but no candidate gene could be identified unequivocally. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users.     

Search for Drosophila genes based on patterned expression of mini-white reporter gene of a P lacW vector in adult eyes

Developmental expression of transduced mini-white(w) gene of Drosophila is sensitive to its flanking genomic enhancers. Taking advantage of this phenomenon, we mobilized a P lacW transposon and screened for new transposant lines which showed patterned expression of the mini-w gene in adult eyes. From a screen of about 1,000 independent P lacW transposant lines on the second chromosome, we identified 7 lines which showed patterned w expression in adult eyes. These P insertions were assigned to engrailed, wingless and teashirt genes based on their chromosomal locations, developmental expression of the lacZ reporter gene, lethal embryonic mutant phenotypes and, finally, their failure to complement the lethal alleles of the respective genetic loci. Our results show that although only a small fraction of the total transposant lines displayed patterned w expression, the genetic loci thus identified are those which play essential roles in pattern formation. Scopes of screens for genetic loci based on w reporter gene expression in adult eyes are discussed.     

Co-transformation of gene expression cassettes via particle bombardment to generate safe transgenic plant without any unwanted DNA

The presence of resistant selectable marker genes and other added DNAs such as the vector backbone sequence in transgenic plant might be an unpredictable hazard to the ecosystem as well as to human health, which have affected the safe assessment of transgenic plants seriously. Using minimal gene expression cassette (containing the promoter, coding region, and terminator) without vector backbone sequence for particle bombardment is the new trend of plant genetic transformation. In the present paper, we co-transformed the selectable marker bar gene cassette and non-selected cecropinB gene cassette into rice (Oryza sativa L.) by particle bombardment, then eliminated the selectable marker bar gene in R₁ generation applying the hereditary segregation strategy and attained two safe transgenic plants only harboring cecropinB gene cassettes without any superfluous DNA. This is the fist report indicating that the combination of minimal gene cassettes transformation with the co-transformation and segregation strategy can generate selectable marker-free transgenic plants, which will promote the advancement in plant genetic engineering greatly.     

A reporter gene system used to study developmental expressionof alternative oxidase and isolate mitochondrial retrograde regulationmutants in Arabidopsis

Perturbation of mitochondrial function causes altered nuclear gene expression in plants. To study this response, called mitochondrial retrograde regulation, and developmental gene expression, a transgenic Arabidopsis thaliana (Col-0) line containing a firefly luciferase gene controlled by a promoter region of the Arabidopsis alternative oxidase 1a gene (AtAOX1a) was created. The transgene and the endogenous gene were developmentally induced in young cotyledons to a level higher than in older cotyledons and leaves. Analysis of transgene expression suggests that this is true for emerging leaves as well. Antimycin A (AA), a mitochondrial electron transport chain inhibitor, and monofluroacetate (MFA), a TCA cycle inhibitor, induced expression of the transgene and the endogenous gene in parallel. The following comparative responses of the transgene to inhibitors were observed: (a) the response in cotyledons to AA treatment differed greatly in magnitude from the response in leaves; (b) the induction kinetics in cotyledons following MFA treatment differed greatly from the kinetics in leaves; and (c) the induction kinetics following MFA treatment differed from the kinetics of AA in both leaves and cotyledons. The transgenic line was used in a genetic screen to isolate mutants with greatly decreased transgene and AtAOX1a induction in response to AA. Some of these mutant lines showed greatly decreased induction by MFA, but one did not. Taken altogether, the data provide genetic evidence that suggests that induction of the AtAOX1a gene by distinct mitochondrial perturbations are via distinct, but overlapping signaling pathways that are tissue specific.     

Native xylose-inducible promoter expands the genetic tools for the biomass-degrading,extremely thermophilic bacterium Caldicellulosiruptor bescii

Regulated control of both homologous and heterologous gene expression is essential for precise genetic manipulation and metabolic engineering of target microorganisms. However, there are often no options available for inducible promoters when working with non-model microorganisms. These include extremely thermophilic, cellulolytic bacteria that are of interest for renewable lignocellulosic conversion to biofuels and chemicals. In fact, improvements to the genetic systems in these organisms often cease once transformation is achieved. This present study expands the tools available for genetically engineering Caldicellulosiruptor bescii, the most thermophilic cellulose-degrader known growing up to 90 °C on unpretreated plant biomass. A native xylose-inducible (P _xi ) promoter was utilized to control the expression of the reporter gene (ldh) encoding lactate dehydrogenase. The P _xi -ldh construct resulted in a both increased ldh expression (20-fold higher) and lactate dehydrogenase activity (32-fold higher) in the presence of xylose compared to when glucose was used as a substrate. Finally, lactate production during growth of the recombinant C. bescii strain was proportional to the initial xylose concentration, showing that tunable expression of genes is now possible using this xylose-inducible system. This study represents a major step in the use of C. bescii as a potential platform microorganism for biotechnological applications using renewable biomass.

     

Effects of selection for ethanol preference on gene expression in the nucleus accumbens of HS‐CC mice

Previous studies on changes in murine brain gene expression associated with the selection for ethanol preference have used F₂ intercross or heterogeneous stock (HS) founders, derived from standard laboratory strains. However, these populations represent only a small proportion of the genetic variance available in Mus musculus. To investigate a wider range of genetic diversity, we selected mice for ethanol preference using an HS derived from the eight strains of the collaborative cross. These HS mice were selectively bred (four generations) for high and low ethanol preference. The nucleus accumbens shell of naive S₄ mice was interrogated using RNA sequencing (RNA‐Seq). Gene networks were constructed using the weighted gene coexpression network analysis assessing both coexpression and cosplicing. Selection targeted one of the network coexpression modules (greenyellow) that was significantly enriched in genes associated with receptor signaling activity including Chrna7, Grin2a, Htr2a and Oprd1. Connectivity in the module as measured by changes in the hub nodes was significantly reduced in the low preference line. Of particular interest was the observation that selection had marked effects on a large number of cell adhesion molecules, including cadherins and protocadherins. In addition, the coexpression data showed that selection had marked effects on long non‐coding RNA hub nodes. Analysis of the cosplicing network data showed a significant effect of selection on a large cluster of Ras GTPase‐binding genes including Cdkl5, Cyfip1, Ndrg1, Sod1 and Stxbp5. These data in part support the earlier observation that preference is linked to Ras/Mapk pathways.     

Expression profiles of growth-related genes in two Nile tilapia strains and their crossbred provide insights into introgressive breeding effects

The Nile tilapia (Oreochromis niloticus) is a prominent farmed fish in aquaculture worldwide. Crossbreeding has recently been carried out between the Red-Stirling and the wt Chitralada strains of Nile tilapia, producing a heterotic hybrid (7/8 Chitralada and 1/8 Red-Stirling) that combines the superior growth performance of the Chitralada with the reddish coloration of the Red-Stirling strain. While classical selective breeding and crossbreeding strategies are well known, the molecular mechanisms underlying the phenotypic expression of economically advantageous traits in tilapia remain largely unknown. Molecular investigations have shown that variable expression of growth hormone (gh), insulin-like growth factors (igf1 and 2) and somatolactin (smtla) – components of the growth hormone/insulin-like growth factor (GH/IGF) axis – and myostatin (mstn) genes can affect traits of economic relevance in farmed animals. The aim of this study was to assess and compare the gene expression signature among Chitralada, Red-Stirling and their backcross hybrid in order to gain insights into the effects of introgressive breeding in modulation of the GH/IGF axis. Gene expression analyses in distinct tissues showed that most genes of the GH/IGF axis were up-regulated and mstn was down-regulated in backcross animals in comparison with Red-Stirling and Chitralada animals. These gene expression profiles revealed that backcross animals displayed a distinctive expression signature, which attests to the effectiveness of the introgressive breeding technique. Our findings also suggest that the GH/IGF axis and mstn genes might be candidate markers for fish performance and prove useful within genetic improvement programs aimed at the production of superior-quality tilapia strains using introgressive breeding.     

The number of breeders explains genetic connectivity in an endangered bird

Connectivity is central to ecology and evolution as it focuses on the movement of individuals or genes across landscapes. Genetic connectivity approaches aim to understand gene flow but often estimate it indirectly based on metrics of genetic differentiation, which can also be affected by other evolutionary forces such as genetic drift. Gene flow and genetic drift are driven by separate ecological mechanisms with potentially differing effects on genetic differentiation and interpretations of genetic connectivity. The ecological mechanisms contributing to gene flow and genetic drift are primarily effective dispersal, or movement followed by successful reproduction, and the number of breeders in a local population, N_b, respectively. Yet, rarely are these ecological mechanisms and genetic connectivity measured simultaneously across landscapes. We examine the roles of effective dispersal and N_b on genetic connectivity across the entire range of the endangered snail kite (Rostrhamus sociabilis plumbeus), between 2006–2015. We find that both N_b and effective dispersal are important predictors of genetic connectivity across this landscape, but that N_b has a 3 × stronger effect on genetic connectivity. Furthermore, N_b is positively correlated with heterozygosity and allelic richness within patches, suggesting a potentially important role of genetic drift, in addition to gene flow, on genetic connectivity. These results emphasize that conservation efforts should focus on not only between‐patch processes of movement but also within‐patch processes regarding habitat quality and local population size for increasing genetic connectivity.     

Estimating number of transgene copies in transgenic rapeseed by real-time PCR assay with<Emphasis Type="Italic">HMG I/Y</Emphasis> as an endogenous reference gene

In transgenic plants, the number of transgene copies can greatly influence the level of expression and genetic stability of the target gene. Transgene copy numbers are estimated by Southern blot analysis, which is laborious and time-consuming, requires relatively large amounts of plant materials, and may involve hazardous radioisotopes. Here we report the development of a sensitive, convenient real-time PCR technique for estimating the number of transgene copies in transgenic rapeseed. This system uses TaqMan quantitative real-time PCR and comparison with a novel, confirmed single-copy endogenous reference gene, high-mobile-group protein I/Y (HMG I/Y), to determine the numbers of copies of exogenous β-glucuronidase (GUS) and neomycin phosphotransferase II (nptII) genes. TheGUS andnptII copy numbers in primary transformants (T₀) were calculated by comparing threshold cycle (C _T) values of theGUS andnptII genes with those of the internal standard,HMG I/Y. This method is more convenient and accurate than Southern blotting because the number of copies of the exogenous gene could be directly deduced by comparing itsC _T value to that of the single-copy endogenous gene in each sample. Unlike other similar procedures of real-time PCR assay, this method does not require identical amplification efficiencies between the PCR systems for target gene and endogenous reference gene, which can avoid the bias that may result from slight variations in amplification efficiencies between PCR systems of the target and endogenous reference genes.