Genomic approaches to analyzing natural variation in Arabidopsis thaliana

The genomics tools available for studying Arabidopsis thaliana are a great resource for researchers trying to characterize and understand the genetic basis of natural variation. Abundant polymorphic markers aid quantitative trait locus (QTL) mapping, the fully sequenced genome provides rapid identification of candidate loci, and extensive knockout collections allow those candidate loci to be tested. Combining QTL mapping of classic phenotypic traits with biochemical or expression analysis is providing mechanistic insight into the traits of interest. Conversely, natural variation studies are now being done on genomic traits such as methylation or chiasma frequency.     

Genomic analysis of synaptotagmin genes.

I used TBLASTn to probe DNA sequence databases with a consensus peptide sequence corresponding to the most highly conserved region of the rodent synaptotagmin (Syt) gene family, which is within the C2B domain. I found human homologues for all known rodent genes, and found six further human genomic loci which encode potential family members. I found eight potential family members in Caenorhabditis elegans, six in Drosophila melanogaster, and four in Arabidopsis thaliana. The C. elegans Syt1 homologue uniquely encodes two alternative C2B exons, one or the other of which is expressed at a time. Comparison of the genomic structures of the Syt genes makes clear the different phylogenies of the different subgroups. Knowledge of the genomic structures will aid the systematic investigation of alternative splicing in Syt genes.     

Genomic instability is associated with natural life span variation in Saccharomyces cerevisiae

Increasing genomic instability is associated with aging in eukaryotes, but the connection between genomic instability and natural variation in life span is unknown. We have quantified chronological life span and loss-of-heterozygosity (LOH) in 11 natural isolates of Saccharomyces cerevisiae. We show that genomic instability increases and mitotic asymmetry breaks down during chronological aging. The age-dependent increase of genomic instability generally lags behind the drop of viability and this delay accounts for approximately 50% of the observed natural variation of replicative life span in these yeast isolates. We conclude that the abilities of yeast strains to tolerate genomic instability co-vary with their replicative life spans. To the best of our knowledge, this is the first quantitative evidence that demonstrates a link between genomic instability and natural variation in life span.     

Genomic variation in neurons

Neurons born during the fetal period have extreme longevity and survive until the death of the individual because the human brain has highly limited tissue regeneration. The brain is comprised of an enormous variety of neurons each exhibiting different morphological and physiological characteristics and recent studies have further reported variations in their genome including chromosomal abnormalities, copy number variations, and single nucleotide mutations. During the early stages of brain development, the increasing number of neurons generated at high speeds has been proposed to lead to chromosomal instability. Additionally, mutations in the neuronal genome can occur in the mature brain. This observed genomic mosaicism in the brain can be produced by multiple endogenous and environmental factors and careful analyses of these observed variations in the neuronal genome remain central for our understanding of the genetic basis of neurological disorders.     

Genomic analysis of detoxification genes in the mosquito Aedes aegypti

Annotation of the recently determined genome sequence of the major dengue vector, Aedes aegypti, reveals an abundance of detoxification genes. Here, we report the presence of 235 members of the cytochrome P450, glutathione transferase and carboxy/cholinesterase families in Ae. aegypti. This gene count represents an increase of 58% and 36% compared with the fruitfly, Drosophila melanogaster, and the malaria mosquito, Anopheles gambiae, respectively. The expansion is not uniform within the gene families. Secure orthologs can be found across the insect species for enzymes that have presumed or proven biosynthetic or housekeeping roles. In contrast, subsets of these gene families that are associated with general xenobiotic detoxification, in particular the CYP6, CYP9 and alpha esterase families, have expanded in Ae. aegypti. In order to identify detoxification genes associated with resistance to insecticides we constructed an array containing unique oligonucleotide probes for these genes and compared their expression level in insecticide resistant and susceptible strains. Several candidate genes were identified with the majority belonging to two gene families, the CYP9 P450s and the Epsilon GSTs. This 'Ae. aegypti Detox Chip' will facilitate the implementation of insecticide resistance management strategies for arboviral control programmes.     

Population history and natural selection shape patterns of genetic variation in 132 genes

Identifying regions of the human genome that have been targets of natural selection will provide important insights into human evolutionary history and may facilitate the identification of complex disease genes. Although the signature that natural selection imparts on DNA sequence variation is difficult to disentangle from the effects of neutral processes such as population demographic history, selective and demographic forces can be distinguished by analyzing multiple loci dispersed throughout the genome. We studied the molecular evolution of 132 genes by comprehensively resequencing them in 24 African-Americans and 23 European-Americans. We developed a rigorous computational approach for taking into account multiple hypothesis tests and demographic history and found that while many apparent selective events can instead be explained by demography, there is also strong evidence for positive or balancing selection at eight genes in the European-American population, but none in the African-American population. Our results suggest that the migration of modern humans out of Africa into new environments was accompanied by genetic adaptations to emergent selective forces. In addition, a region containing four contiguous genes on Chromosome 7 showed striking evidence of a recent selective sweep in European-Americans. More generally, our results have important implications for mapping genes underlying complex human diseases.     

Genomic analysis of cyclic-di-GMP-related genes in rhizobial type strains and functional analysis in Rhizobium etli

Rhizobia are soil bacteria that can fix nitrogen in symbiosis with leguminous plants or exist free living in the rhizosphere. Crucial to their complex lifestyle is the ability to sense and respond to diverse environmental stimuli, requiring elaborate signaling pathways. In the majority of bacteria, the nucleotide-based second messenger cyclic diguanosine monophosphate (c-di-GMP) is involved in signal transduction. Surprisingly, little is known about the importance of c-di-GMP signaling in rhizobia. We have analyzed the genome sequences of six well-studied type species (Bradyrhizobium japonicum, Mesorhizobium loti, Rhizobium etli, Rhizobium leguminosarum, Sinorhizobium fredii, and Sinorhizobium meliloti) for proteins possibly involved in c-di-GMP signaling based on the presence of four domains: GGDEF (diguanylate cyclase), EAL and HD-GYP (phosphodiesterase), and PilZ (c-di-GMP sensor). We find that rhizobia possess a high number of these proteins. Conservation analysis suggests that c-di-GMP signaling proteins modulate species-specific pathways rather than ancient rhizobia-specific processes. Two hybrid GGDEF-EAL proteins were selected for functional analysis, R. etli RHE_PD00105 (CdgA) and RHE_PD00137 (CdgB). Expression of cdgA and cdgB is repressed by the alarmone (p)ppGpp. cdgB is significantly expressed on plant roots and free living. Mutation of cdgA, cdgB, or both does not affect plant root colonization, nitrogen fixation capacity, biofilm formation, motility, and exopolysaccharide production. However, heterologous expression of the individual GGDEF and EAL domains of each protein in Escherichia coli strongly suggests that CdgA and CdgB are bifunctional proteins, possessing both diguanylate cyclase and phosphodiesterase activities. Taken together, our results provide a platform for future studies of c-di-GMP signaling in rhizobia.     

Genomic analysis of the major bovine milk protein genes.

The genomic arrangement of the major bovine milk protein genes has been determined using a combination of physical mapping techniques. The major milk proteins consist of the four caseins, alpha s1 (CASAS1), alpha s2 (CASAS2), beta (CASB), and kappa (CASK), as well as the two major whey proteins, alpha-lactalbumin (LALBA) and beta-lactoglobulin (LGB). A panel of bovine X hamster hybrid somatic cells analyzed for the presence or absence of bovine specific restriction fragments revealed the genes coding for the major milk proteins to reside on three chromosomes. The four caseins were assigned to syntenic group U15 and localized to bovine chromosome 6 at q31-33 by in situ hybridization. LALBA segregated with syntenic group U3, while LGB segregated with U16. Pulsed-field gel electrophoresis confirmed genetic mapping results indicating tight linkage of the casein genes. The four genes reside on less than 200 kb of DNA in the order CASAS1-CASB-CASAS2-CASK. Multiple restriction fragment length polymorphisms were also found at the six loci in three breeds of cattle.     

Regulatory components involved in altering lateral root development in response to localized iron: Evidence for natural genetic variation

On the search for sparingly available nutrients, plants may alter their root architecture to improve soil exploration. So far, the examples for root system modifications induced by a heterogeneous availability of nutrients have been reported for the macronutrients nitrogen (N) and phosphorous (P). In an attempt to extend this type of knowledge to other nutrients, we recently provided evidence that Arabidopsis roots are able to sense a local availability of the micronutrient iron (Fe) and to respond with lateral root elongation into the Fe-containing patch. This specific root response was caused by enhanced elongation of cells leaving the root meristem and was dependent on an AUX1-mediated auxin accumulation in the lateral root apices. In this report, we compare mechanisms underlying this response with those known for other nutrients and show that a substantial genotypic variation exists among accessions of A. thaliana in the responsiveness of lateral roots toward localized Fe supplies.     

Genomic convergence and network analysis approach to identify candidate genes in Alzheimer's disease

Background

Alzheimer’s disease (AD) is one of the leading genetically complex and heterogeneous disorder that is influenced by both genetic and environmental factors. The underlying risk factors remain largely unclear for this heterogeneous disorder. In recent years, high throughput methodologies, such as genome-wide linkage analysis (GWL), genome-wide association (GWA) studies, and genome-wide expression profiling (GWE), have led to the identification of several candidate genes associated with AD. However, due to lack of consistency within their findings, an integrative approach is warranted. Here, we have designed a rank based gene prioritization approach involving convergent analysis of multi-dimensional data and protein-protein interaction (PPI) network modelling.

Results

Our approach employs integration of three different AD datasets- GWL,GWA and GWE to identify overlapping candidate genes ranked using a novel cumulative rank score (S_R) based method followed by prioritization using clusters derived from PPI network. S_R for each gene is calculated by addition of rank assigned to individual gene based on either p value or score in three datasets. This analysis yielded 108 plausible AD genes. Network modelling by creating PPI using proteins encoded by these genes and their direct interactors resulted in a layered network of 640 proteins. Clustering of these proteins further helped us in identifying 6 significant clusters with 7 proteins (EGFR, ACTB, CDC2, IRAK1, APOE, ABCA1 and AMPH) forming the central hub nodes. Functional annotation of 108 genes revealed their role in several biological activities such as neurogenesis, regulation of MAP kinase activity, response to calcium ion, endocytosis paralleling the AD specific attributes. Finally, 3 potential biochemical biomarkers were found from the overlap of 108 AD proteins with proteins from CSF and plasma proteome. EGFR and ACTB were found to be the two most significant AD risk genes.

Conclusions

With the assumption that common genetic signals obtained from different methodological platforms might serve as robust AD risk markers than candidates identified using single dimension approach, here we demonstrated an integrated genomic convergence approach for disease candidate gene prioritization from heterogeneous data sources linked to AD.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-199) contains supplementary material, which is available to authorized users.     

Candidate genes and QTLs for fruit ripening and softening in melon

Different factors affect the quality of melon fruit and among them long shelf life is critical from the consumer’s point of view. In melon, cultivars showing both climacteric and non-climacteric ripening types are found. In this study we have investigated climacteric ripening and fruit softening using a collection of near-isogenic lines (NILs) derived from the non-climacteric melon parental lines PI 161375 (SC) and “Piel de Sapo” (PS). Surprisingly, we found that QTL eth3.5 in NIL SC3-5b induced a climacteric-ripening phenotype with increased respiration and ethylene levels. Data suggest that the non-climacteric phenotypes from PI 161375 and “Piel de Sapo” may be the result of mutations in different genes. Several QTLs for fruit flesh firmness were also detected. Candidate genes putatively involved in ethylene regulation, biosynthesis and perception and cell wall degradation were mapped and some colocations with QTLs were observed. These results may provide additional data towards understanding of non-climacteric ripening in melon.     

Genomic sequence around butterfly wing development genes: annotation and comparative analysis

Background

Analysis of genomic sequence allows characterization of genome content and organization, and access beyond gene-coding regions for identification of functional elements. BAC libraries, where relatively large genomic regions are made readily available, are especially useful for species without a fully sequenced genome and can increase genomic coverage of phylogenetic and biological diversity. For example, no butterfly genome is yet available despite the unique genetic and biological properties of this group, such as diversified wing color patterns. The evolution and development of these patterns is being studied in a few target species, including Bicyclus anynana, where a whole-genome BAC library allows targeted access to large genomic regions.

Methodology/Principal Findings

We characterize ∼1.3 Mb of genomic sequence around 11 selected genes expressed in B. anynana developing wings. Extensive manual curation of in silico predictions, also making use of a large dataset of expressed genes for this species, identified repetitive elements and protein coding sequence, and highlighted an expansion of Alcohol dehydrogenase genes. Comparative analysis with orthologous regions of the lepidopteran reference genome allowed assessment of conservation of fine-scale synteny (with detection of new inversions and translocations) and of DNA sequence (with detection of high levels of conservation of non-coding regions around some, but not all, developmental genes).

Conclusions

The general properties and organization of the available B. anynana genomic sequence are similar to the lepidopteran reference, despite the more than 140 MY divergence. Our results lay the groundwork for further studies of new interesting findings in relation to both coding and non-coding sequence: 1) the Alcohol dehydrogenase expansion with higher similarity between the five tandemly-repeated B. anynana paralogs than with the corresponding B. mori orthologs, and 2) the high conservation of non-coding sequence around the genes wingless and Ecdysone receptor, both involved in multiple developmental processes including wing pattern formation.     

Genomic analysis of the relationship between gene expression variation and DNA polymorphism in Drosophila simulans

Background

Understanding how DNA sequence polymorphism relates to variation in gene expression is essential to connecting genotypic differences with phenotypic differences among individuals. Addressing this question requires linking population genomic data with gene expression variation.

Results

Using whole genome expression data and recent light shotgun genome sequencing of six Drosophila simulans genotypes, we assessed the relationship between expression variation in males and females and nucleotide polymorphism across thousands of loci. By examining sequence polymorphism in gene features, such as untranslated regions and introns, we find that genes showing greater variation in gene expression between genotypes also have higher levels of sequence polymorphism in many gene features. Accordingly, X-linked genes, which have lower sequence polymorphism levels than autosomal genes, also show less expression variation than autosomal genes. We also find that sex-specifically expressed genes show higher local levels of polymorphism and divergence than both sex-biased and unbiased genes, and that they appear to have simpler regulatory regions.

Conclusion

The gene-feature-based analyses and the X-to-autosome comparisons suggest that sequence polymorphism in cis-acting elements is an important determinant of expression variation. However, this relationship varies among the different categories of sex-biased expression, and trans factors might contribute more to male-specific gene expression than cis effects. Our analysis of sex-specific gene expression also shows that female-specific genes have been overlooked in analyses that only point to male-biased genes as having unusual patterns of evolution and that studies of sexually dimorphic traits need to recognize that the relationship between genetic and expression variation at these traits is different from the genome as a whole.     

Towards the cloning of genes underlying murine QTLs

Many quantitative trait loci (QTLs), including those for ethanol-related traits, have been mapped in the mouse. In light of rapidly developing tools and resources, we briefly review the strategy for identifying the genes underlying these QTLs. We note that positional cloning will soon be a matter of testing candidate genes rather than discovering genes; therefore, we describe a ``congenic test' to support that a candidate gene is indeed a QTL. Considering the rapid development of congenics and mutants, we also identify four areas of investigation—phenotypes, ethanol specificity, environment, and gene interactions—that might be exploited during the course of positional cloning to gain insights into QTL pathways. In particular, we note that multiple mutants of nearly every major neurotransmitter pathway have now been made. These mutants are not only useful for phenotypic tests, but also could be used to conduct ``gene dependence' tests of QTLs. We also consider potential applications for the very recently developed ability to clone mice. Received: 15 September 1998 / Accepted: 8 October 1998     

Genetic and behavioral analysis of natural variation in Drosophila melanogaster pupation position

Drosophila melanogaster pupae are exposed to many biotic and abiotic dangers while immobilized during several days of metamorphosis. As a passive defense mechanism, appropriate pupation site selection represents an important mitigation of these threats. Pupation site selection is sensitive to genetic and environmental influences, but the specific mechanisms of the behavior are largely unknown. Using a set of 76 recombinant inbred strains we identify a single quantitative trait locus, at polytene position 56A01-C11, associated with pupation site variation. We furthermore present a detailed investigation into the wandering behaviors of two strains expressing different pupation position tendencies, and identify behavioral differences. Larvae from a strain that tends to pupate relatively far from the food also tend to travel significantly farther from the media during wandering. We did not observe consistent differences in either the number or duration of wandering forays made by near or far pupating strains. The ability of larvae to integrate several internal and external environmental cues while choosing a contextually appropriate pupation site, and specifically, the variation in this ability, presents a very interesting behavioral phenotype in this highly tractable genetic model organism.     

Modeling stochastic noise in gene regulatory systems

The Master equation is considered the gold standard for modeling the stochastic mechanisms of gene regulation in molecular detail, but it is too complex to solve exactly in most cases, so approximation and simulation methods are essential. However, there is still a lack of consensus about the best way to carry these out. To help clarify the situation, we review Master equation models of gene regulation, theoretical approximations based on an expansion method due to N.G. van Kampen and R. Kubo, and simulation algorithms due to D.T. Gillespie and P. Langevin. Expansion of the Master equation shows that for systems with a single stable steady-state, the stochastic model reduces to a deterministic model in a first-order approximation. Additional theory, also due to van Kampen, describes the asymptotic behavior of multistable systems. To support and illustrate the theory and provide further insight into the complex behavior of multistable systems, we perform a detailed simulation study comparing the various approximation and simulation methods applied to synthetic gene regulatory systems with various qualitative characteristics. The simulation studies show that for large stochastic systems with a single steady-state, deterministic models are quite accurate, since the probability distribution of the solution has a single peak tracking the deterministic trajectory whose variance is inversely proportional to the system size. In multistable stochastic systems, large fluctuations can cause individual trajectories to escape from the domain of attraction of one steady-state and be attracted to another, so the system eventually reaches a multimodal probability distribution in which all stable steady-states are represented proportional to their relative stability. However, since the escape time scales exponentially with system size, this process can take a very long time in large systems.