Genetic Basis Of Natural Variation Research Articles

Genetic architecture of natural variation in Drosophila melanogaster aggressive behavior.

Aggression is an evolutionarily conserved complex behavior essential for survival and the organization of social hierarchies. With the exception of genetic variants associated with bioamine signaling, which have been implicated in aggression in many species, the genetic basis of natural variation in aggression is largely unknown. Drosophila melanogaster is a favorable model system for exploring the genetic basis of natural variation in aggression. Here, we performed genome-wide association analyses using the inbred, sequenced lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and replicate advanced intercross populations derived from the most and least aggressive DGRP lines. We identified genes that have been previously implicated in aggressive behavior as well as many novel loci, including gustatory receptor 63a (Gr63a), which encodes a subunit of the receptor for CO2, and genes associated with development and function of the nervous system. Although genes from the two association analyses were largely nonoverlapping, they mapped onto a genetic interaction network inferred from an analysis of pairwise epistasis in the DGRP. We used mutations and RNAi knock-down alleles to functionally validate 79% of the candidate genes and 75% of the candidate epistatic interactions tested. Epistasis for aggressive behavior causes cryptic genetic variation in the DGRP that is revealed by changing allele frequencies in the outbred populations derived from extreme DGRP lines. This phenomenon may pertain to other fitness traits and species, with implications for evolution, applied breeding, and human genetics.

Longevity GWAS Using the Drosophila Genetic Reference Panel.

We used 197 Drosophila melanogaster Genetic Reference Panel (DGRP) lines to perform a genome-wide association analysis for virgin female lifespan, using ~2M common single nucleotide polymorphisms (SNPs). We found considerable genetic variation in lifespan in the DGRP, with a broad-sense heritability of 0.413. There was little power to detect signals at a genome-wide level in single-SNP and gene-based analyses. Polygenic score analysis revealed that a small proportion of the variation in lifespan (~4.7%) was explicable in terms of additive effects of common SNPs (≥2% minor allele frequency). However, several of the top associated genes are involved in the processes previously shown to impact ageing (eg, carbohydrate-related metabolism, regulation of cell death, proteolysis). Other top-ranked genes are of unknown function and provide promising candidates for experimental examination. Genes in the target of rapamycin pathway (TOR; Chrb, slif, mipp2, dredd, RpS9, dm) contributed to the significant enrichment of this pathway among the top-ranked 100 genes (p = 4.79×10−06). Gene Ontology analysis suggested that genes involved in carbohydrate metabolism are important for lifespan; including the InterPro term DUF227, which has been previously associated with lifespan determination. This analysis suggests that our understanding of the genetic basis of natural variation in lifespan from induced mutations is incomplete.

High Levels of Diversity in the Phytochemistry, Ploidy and Genetics of the Medicinal Plant Acorus calamus L.

Acorus calamus L. is an emergent aquatic macrophyte that is used in several medicines and food supplements. A. calamus ecotypes with low concentrations of β-asarone can be used for medicinal purposes. We investigated 50 populations of this species, representing the entire natural distribution range in Indian Himalayan Region (IHR), to understand natural variations in β-asarone concentration, ploidy level and genetic divergence. Our results showed that: 1) all four cytotypes (diploid, triploid, tetraploid and hexaploid) of A. calamus exist in natural populations in the IHR; 2) the β-asarone concentration is highly correlated with the ploidy level (diploid and hexaploid populations exhibit low β-asarone concentrations, whereas triploids and tetraploids show high β-asarone concentrations); 3) there is a significant correlation between the β-asarone concentration and its distribution along latitudinal and longitudinal ranges, but not ecological factors (mean annual temperature and mean annual precipitation); and 4) geographic boundaries play an important role in the population divergence of A. calamus. This study provides a foundation for understanding the genetic basis of natural variation in complex traits of A. calamus and for devising conservation strategies for this medically important plant.

The genetic basis of natural variation in seed size and seed number and their trade-off using Arabidopsis thaliana MAGIC lines.

Offspring number and size are key traits determining an individual's fitness and a crop's yield. Yet, extensive natural variation within species is observed for these traits. Such variation is typically explained by trade-offs between fecundity and quality, for which an optimal solution is environmentally dependent. Understanding the genetic basis of seed size and number, as well as any possible genetic constraints preventing the maximization of both, is crucial from both an evolutionary and applied perspective. We investigated the genetic basis of natural variation in seed size and number using a set of Arabidopsis thaliana multiparent advanced generation intercross (MAGIC) lines. We also tested whether life history affects seed size, number, and their trade-off. We found that both seed size and seed number are affected by a large number of mostly nonoverlapping QTL, suggesting that seed size and seed number can evolve independently. The allele that increases seed size at most identified QTL is from the same natural accession, indicating past occurrence of directional selection for seed size. Although a significant trade-off between seed size and number is observed, its expression depends on life-history characteristics, and generally explains little variance. We conclude that the trade-off between seed size and number might have a minor role in explaining the maintenance of variation in seed size and number, and that seed size could be a valid target for selection.

Genome-wide delineation of natural variation for pod shatter resistance in Brassica napus.

Resistance to pod shattering (shatter resistance) is a target trait for global rapeseed (canola, Brassica napus L.), improvement programs to minimise grain loss in the mature standing crop, and during windrowing and mechanical harvest. We describe the genetic basis of natural variation for shatter resistance in B. napus and show that several quantitative trait loci (QTL) control this trait. To identify loci underlying shatter resistance, we used a novel genotyping-by-sequencing approach DArT-Seq. QTL analysis detected a total of 12 significant QTL on chromosomes A03, A07, A09, C03, C04, C06, and C08; which jointly account for approximately 57% of the genotypic variation in shatter resistance. Through Genome-Wide Association Studies, we show that a large number of loci, including those that are involved in shattering in Arabidopsis, account for variation in shatter resistance in diverse B. napus germplasm. Our results indicate that genetic diversity for shatter resistance genes in B. napus is limited; many of the genes that might control this trait were not included during the natural creation of this species, or were not retained during the domestication and selection process. We speculate that valuable diversity for this trait was lost during the natural creation of B. napus. To improve shatter resistance, breeders will need to target the introduction of useful alleles especially from genotypes of other related species of Brassica, such as those that we have identified.

Genome-wide analysis in Drosophila reveals age-specific effects of SNPs on fitness traits.

Most organisms exhibit senescence; a decline in physiological function with age. In nature, rates of senescence vary extensively among individuals and this variation has a significant genetic component; however, we know little about the genes underlying senescence. Here we show the first evidence that individual alleles influence fecundity in an age-specific manner and so the genetic basis of natural variation in fecundity changes dramatically with age. We complete a genome-wide association to identify single-nucleotide polymorphisms (SNPs) affecting lifespan and age-specific fecundity using the Drosophila melanogaster Genetic Reference Panel. We identify 1,031 SNPs affecting fecundity and 52 influencing lifespan. Only one SNP is associated with both early- and late-age fecundity. The age-specific effect of candidate genes on fecundity is validated using RNA interference. In addition, there is a dramatic increase in the number of SNPs influencing fecundity with age. This result provides support for the mutation accumulation theory of aging.

Genetic control of heterochrony in Eucalyptus globulus.

A change in the timing or rate of developmental events throughout ontogeny is referred to as heterochrony, and it is a major evolutionary process in plants and animals. We investigated the genetic basis for natural variation in the timing of vegetative phase change in the tree Eucalyptus globulus, which undergoes a dramatic change in vegetative morphology during the juvenile-to-adult transition. Quantitative trait loci analysis in an outcross F2 family derived from crosses between individuals from a coastal population of E. globulus with precocious vegetative phase change and individuals from populations in which vegetative phase change occurs several years later implicated the microRNA EglMIR156.5 as a potential contributor to this heterochronic difference. Additional evidence for the involvement of EglMIR156.5 was provided by its differential expression in trees with early and late phase change. Our findings suggest that changes in the expression of miR156 underlie natural variation in vegetative phase change in E. globulus, and may also explain interspecific differences in the timing of this developmental transition.

The genetic basis of natural variation for iron homeostasis in the maize IBM population.

BackgroundIron (Fe) deficiency symptoms in maize (Zea mays subsp. mays) express as leaf chlorosis, growth retardation, as well as yield reduction and are typically observed when plants grow in calcareous soils at alkaline pH. To improve our understanding of genotypical variability in the tolerance to Fe deficiency-induced chlorosis, the objectives of this study were to (i) determine the natural genetic variation of traits related to Fe homeostasis in the maize intermated B73 × Mo17 (IBM) population, (ii) to identify quantitative trait loci (QTLs) for these traits, and (iii) to analyze expression levels of genes known to be involved in Fe homeostasis as well as of candidate genes obtained from the QTL analysis.ResultsIn hydroponically-grown maize, a total of 47 and 39 QTLs were detected for the traits recorded under limited and adequate supply of Fe, respectively.ConclusionsFrom the QTL results, we were able to identify new putative candidate genes involved in Fe homeostasis under a deficient or adequate Fe nutritional status, like Ferredoxin class gene, putative ferredoxin PETF, metal tolerance protein MTP4, and MTP8. Furthermore, our expression analysis of candidate genes suggested the importance of trans-acting regulation for 2’-deoxymugineic acid synthase 1 (DMAS1), nicotianamine synthase (NAS3, NAS1), formate dehydrogenase 1 (FDH1), methylthioribose-1-phosphate isomerase (IDI2), aspartate/tyrosine/aromatic aminotransferase (IDI4), and methylthioribose kinase (MTK).

Exploring the effects of gene dosage on mandible shape in mice as a model for studying the genetic basis of natural variation

Mandible shape in the mouse is a complex trait that is influenced by many genetic factors. However, little is known about the action of single genes on adult mandible shape so far, since most developmentally relevant genes are already required during embryogenesis, i.e., knockouts lead to embryonic death or severe deformations, before the mandible is fully formed. We employ here a geometric morphometric approach to identify subtle phenotypic differences caused by dosage effects of candidate genes. We use mouse strains with specific gene modifications (knockouts and knockins) to compare heterozygous animals with controls from the same stock, which is expected to be equivalent to a change of gene expression of the respective locus. Such differences in expression level are also likely to occur as part of the natural variation. We focus on Bmp pathway genes (Bmp4, its antagonist Noggin, and combinations of Bmp5-7 genotypes), but include also two other developmental control genes suspected to affect mandible development in some way (Egfr and Irf6). In addition, we study the effects of Hoxd13, as well as an extracellular matrix constituent (Col2a1). We find that subtle but significant shape differences are caused by differences in gene dosage of several of these genes. The changes seen for Bmp4 and Noggin are partially compatible with the action of these genes known from birds and fish. We find significant shape changes also for Hoxd13, although this gene has so far only been implicated in skeletal patterning processes of the limbs. Comparing the effect sizes of gene dosage changes to the variation found in natural populations of mice as well as quantitative trait loci (QTL) effects on mandible shape, we find that the effect sizes caused by gene dosage changes are at the lower end of the spectrum of natural variation, but larger than the average additive effects found in QTL studies. We conclude that studying gene dosage effects have the potential to provide new insights into aspects of craniofacial development, variation, and evolution.

Genome-wide association study of sleep in Drosophila melanogaster.

BackgroundSleep is a highly conserved behavior, yet its duration and pattern vary extensively among species and between individuals within species. The genetic basis of natural variation in sleep remains unknown.ResultsWe used the Drosophila Genetic Reference Panel (DGRP) to perform a genome-wide association (GWA) study of sleep in D. melanogaster. We identified candidate single nucleotide polymorphisms (SNPs) associated with differences in the mean as well as the environmental sensitivity of sleep traits; these SNPs typically had sex-specific or sex-biased effects, and were generally located in non-coding regions. The majority of SNPs (80.3%) affecting sleep were at low frequency and had moderately large effects. Additive models incorporating multiple SNPs explained as much as 55% of the genetic variance for sleep in males and females. Many of these loci are known to interact physically and/or genetically, enabling us to place them in candidate genetic networks. We confirmed the role of seven novel loci on sleep using insertional mutagenesis and RNA interference.ConclusionsWe identified many SNPs in novel loci that are potentially associated with natural variation in sleep, as well as SNPs within genes previously known to affect Drosophila sleep. Several of the candidate genes have human homologues that were identified in studies of human sleep, suggesting that genes affecting variation in sleep are conserved across species. Our discovery of genetic variants that influence environmental sensitivity to sleep may have a wider application to all GWA studies, because individuals with highly plastic genotypes will not have consistent phenotypes.

The genetic basis for natural variation in alcohol sensitivity in drosophila

The genetic basis for natural variation in alcohol sensitivity in drosophila

Rapid creation of Arabidopsis doubled haploid lines for quantitative trait locus mapping.

Quantitative trait loci (QTL) mapping is a powerful tool for investigating the genetic basis of natural variation. QTL can be mapped using a number of different population designs, but recombinant inbred lines (RILs) are among the most effective. Unfortunately, homozygous RIL populations are time consuming to construct, typically requiring at least six generations of selfing starting from a heterozygous F(1). Haploid plants produced from an F(1) combine the two parental genomes and have only one allele at every locus. Converting these sterile haploids into fertile diploids (termed "doubled haploids," DHs) produces immortal homozygous lines in only two steps. Here we describe a unique technique for rapidly creating recombinant doubled haploid populations in Arabidopsis thaliana: centromere-mediated genome elimination. We generated a population of 238 doubled haploid lines that combine two parental genomes and genotyped them by reduced representation Illumina sequencing. The recombination rate and parental allele frequencies in our population are similar to those found in existing RIL sets. We phenotyped this population for traits related to flowering time and for petiole length and successfully mapped QTL controlling each trait. Our work demonstrates that doubled haploid populations offer a rapid, easy alternative to RILs for Arabidopsis genetic analysis.

Mapping of ionomic traits in Mimulus guttatus reveals Mo and Cd QTLs that colocalize with MOT1 homologues.

Natural variation in the regulation of the accumulation of mineral nutrients and trace elements in plant tissues is crucial to plant metabolism, development, and survival across different habitats. Studies of the genetic basis of natural variation in nutrient metabolism have been facilitated by the development of ionomics. Ionomics is a functional genomic approach for the identification of the genes and gene networks that regulate the elemental composition, or ionome, of an organism. In this study, we evaluated the genetic basis of divergence in elemental composition between an inland annual and a coastal perennial accession of Mimulus guttatus using a recombinant inbred line (RIL) mapping population. Out of 20 elements evaluated, Mo and Cd were the most divergent in accumulation between the two accessions and were highly genetically correlated in the RILs across two replicated experiments. We discovered two major quantitative trait loci (QTL) for Mo accumulation, the largest of which consistently colocalized with a QTL for Cd accumulation. Interestingly, both Mo QTLs also colocalized with the two M. guttatus homologues of MOT1, the only known plant transporter to be involved in natural variation in molybdate uptake.

Multiple quantitative trait loci influence intra-specific variation in genital morphology between phylogenetically distinct lines of Drosophila montana

The evolution of animal genitalia has gained renewed interest because of their potential roles during sexual selection and early stages of species formation. Although central to understanding the evolutionary process, knowledge of the genetic basis of natural variation in genital morphology is limited to a very few species. Using an outbred cross between phylogenetically distinct lines of Drosophila montana, we characterized quantitative trait loci (QTLs) affecting the size and shape of the distiphallus, a prominent part of the male intromittent organ. Our microsatellite-based linkage analysis shows that intra-specific variation in the distiphallus involves several QTLs of largely additive effect and that a highly significant QTL co-localizes with the same inversion where we have earlier localized a large QTL for a sexually selected courtship song trait. The latter indicates that inversions can play an important role in shaping the evolution of rapidly evolving traits with a potential influence on speciation.

Population-Based Resequencing of Experimentally Evolved Populations Reveals the Genetic Basis of Body Size Variation in Drosophila melanogaster

Body size is a classic quantitative trait with evolutionarily significant variation within many species. Locating the alleles responsible for this variation would help understand the maintenance of variation in body size in particular, as well as quantitative traits in general. However, successful genome-wide association of genotype and phenotype may require very large sample sizes if alleles have low population frequencies or modest effects. As a complementary approach, we propose that population-based resequencing of experimentally evolved populations allows for considerable power to map functional variation. Here, we use this technique to investigate the genetic basis of natural variation in body size in Drosophila melanogaster. Significant differentiation of hundreds of loci in replicate selection populations supports the hypothesis that the genetic basis of body size variation is very polygenic in D. melanogaster. Significantly differentiated variants are limited to single genes at some loci, allowing precise hypotheses to be formed regarding causal polymorphisms, while other significant regions are large and contain many genes. By using significantly associated polymorphisms as a priori candidates in follow-up studies, these data are expected to provide considerable power to determine the genetic basis of natural variation in body size.

The genetic basis of divergent pigment patterns in juvenile threespine sticklebacks

Animal pigment patterns are important for a range of functions, including camouflage and communication. Repeating pigment patterns, such as stripes, bars and spots have been of particular interest to developmental and theoretical biologists, but the genetic basis of natural variation in such patterns is largely unexplored. In this study, we identify a difference in a periodic pigment pattern among juvenile threespine sticklebacks (Gasterosteus aculeatus) from different environments. Freshwater sticklebacks exhibit prominent vertical bars that visually break up the body shape, but sticklebacks from marine populations do not. We hypothesize that these distinct pigment patterns are tuned to provide crypsis in different habitats. This phenotypic difference is widespread and appears in most of the freshwater populations that we sampled. We used quantitative trait locus (QTL) mapping in freshwater-marine F2 hybrids to elucidate the genetic architecture underlying divergence in this pigmentation pattern. We identified two QTL that were significantly associated with variation in barring. Interestingly, these QTL were associated with two distinct aspects of the pigment pattern: melanophore number and overall pigment level. We compared the QTL locations with positions of known pigment candidate genes in the stickleback genome. We also identified two major QTL for juvenile body size, providing new insights into the genetic basis of juvenile growth rates in natural populations. In summary, although there is a growing literature describing simple genetic bases for adaptive coloration differences, this study emphasizes that pigment patterns can also possess a more complex genetic architecture.

Novel insights into the genomics of rice root adaptive development

Deciphering the genetic and molecular mechanisms controlling the development of the root system and its adaptive plasticity under adverse environments is of primary importance for the sustainable establishment of the rice crop. Rice displays a complex root structure comprising several root types mostly of postembryonic origin. The large natural variation in root architecture among cultivars reflects their adaptation to contrasting agro-environmental conditions. This article reviews the current knowledge on the organization and anatomy of the various types of roots of the fibrous root system of rice, the diversity and genetic basis of natural variation of root system architecture and performance, and the molecular mechanisms underlying constitutive and adaptive root development. This paper also throws light on how the integrated approach of new tools in high-resolution microscopy imaging, expression profiling, mutant screening, and reverse genetics could facilitate the rapid discovery and analysis of the key genes and regulatory networks involved in root architectural traits affecting plant performance under field conditions.

Sleuthing the difference a nucleotide can make

Identifying the genetic basis of natural variation in genome-wide gene expression is a daunting task, yet central to characterizing the potential evolutionary impact of gene expression polymorphism. In this issue, Brown et al. reveal the genetic variation that underlies a significant portion of the gene expression polymorphism that segregates in a natural isolate of the wine and bread yeast Saccharomyces cerevisiae. That polymorphism is shown to be associated with a single nucleotide insertion. Interestingly, the major expression polymorphism is not caused by an insertion in a cis-regulatory region or in a transcription factor, but in a homonucleotide repeat within the coding sequence of SSY1, which encodes an amino acid sensor. The expression polymorphism, perhaps unsurprisingly, is also environment-dependent and leads to differential fitness across environments. The result calls for a deeper consideration within molecular ecological genomics of the gene-by-environment interactions that lead to differential gene expression. The finding also calls for further data, clarifying the kinds of genetic variation that constitute gene expression polymorphism between individuals in natural populations.

Differential Expression of Genes Important for Adaptation inCapsella bursa-pastoris(Brassicaceae)

Understanding the genetic basis of natural variation is of primary interest for evolutionary studies of adaptation. In Capsella bursa-pastoris, a close relative of Arabidopsis (Arabidopsis thaliana), variation in flowering time is correlated with latitude, suggestive of an adaptation to photoperiod. To identify pathways regulating natural flowering time variation in C. bursa-pastoris, we have studied gene expression differences between two pairs of early- and late-flowering C. bursa-pastoris accessions and compared their response to vernalization. Using Arabidopsis microarrays, we found a large number of significant differences in gene expression between flowering ecotypes. The key flowering time gene FLOWERING LOCUS C (FLC) was not differentially expressed prior to vernalization. This result is in contrast to those in Arabidopsis, where most natural flowering time variation acts through FLC. However, the gibberellin and photoperiodic flowering pathways were significantly enriched for gene expression differences between early- and late-flowering C. bursa-pastoris. Gibberellin biosynthesis genes were down-regulated in late-flowering accessions, whereas circadian core genes in the photoperiodic pathway were differentially expressed between early- and late-flowering accessions. Detailed time-series experiments clearly demonstrated that the diurnal rhythm of CIRCADIAN CLOCK-ASSOCIATED1 (CCA1) and TIMING OF CAB EXPRESSION1 (TOC1) expression differed between flowering ecotypes, both under constant light and long-day conditions. Differential expression of flowering time genes was biologically validated in an independent pair of flowering ecotypes, suggesting a shared genetic basis or parallel evolution of similar regulatory differences. We conclude that genes involved in regulation of the circadian clock, such as CCA1 and TOC1, are strong candidates for the evolution of adaptive flowering time variation in C. bursa-pastoris.

Quantitative Trait Loci for Locomotor Behavior inDrosophila melanogaster

Locomotion is an integral component of most animal behaviors and many human diseases and disorders are associated with locomotor deficits, but little is known about the genetic basis of natural variation in locomotor behavior. Locomotion is a complex trait, with variation attributable to the joint segregation of multiple interacting quantitative trait loci (QTL), with effects that are sensitive to the environment. We assessed variation in a component of locomotor behavior (locomotor reactivity) in a population of 98 recombinant inbred lines of Drosophila melanogaster and mapped four QTL affecting locomotor reactivity by linkage to polymorphic roo transposable element insertion sites. We used complementation tests of deficiencies to fine map these QTL to 12 chromosomal regions and complementation tests of mutations to identify 13 positional candidate genes affecting locomotor reactivity, including Dopa decarboxylase (Ddc), which catalyzes the final step in the synthesis of serotonin and dopamine. Linkage disequilibrium mapping in a population of 164 second chromosome substitution lines derived from a single natural population showed that polymorphisms at Ddc were associated with naturally occurring genetic variation in locomotor behavior. These data implicate variation in the synthesis of bioamines as a factor contributing to natural variation in locomotor reactivity.