Heterozygote advantage and the maintenance of polymorphism for multilocus traits

Explaining how polymorphism is maintained in the face of selection remains a puzzle since selection tends to erode genetic variation. Provided an infinitely large unsubdivided population and no frequency-dependance of selective values, heterozygote advantage is the text book explanation for the maintenance of polymorphism when selection acts at a diallelic locus. Here, we investigate whether this remains true when selection acts at multiple diallelic loci. We use five different definitions of heterozygote advantage that largely cover this concept for multiple loci. Using extensive numerical simulations, we found no clear associations between the presence of any of the five definitions of heterozygote advantage and the maintenance of polymorphism at all loci. The strength of the association decreases as the number of loci increases or as recombination decreases. We conclude that heterozygote advantage cannot be a general mechanism for the maintenance of genetic polymorphism at multiple loci. These findings suggest that a correlation between the number of heterozygote loci and fitness is not warranted on theoretical ground.     

Development and characterization of eleven polymorphic microsatellite loci from South China field mouse (Apodemus draco)

One population distributed in Yunnan of China was regarded as a new species based on mitochondrial cytochrome b sequences. However, the usefulness of mitochondrial sequence data in determining species boundaries is not universally agreed upon, the frequency data from multiple nuclear gene loci is necessary in determining species boundaries. So, we describe in this paper the isolation and characterization of eleven microsatellite loci in the South China field mouse from genomic DNA-enriched libraries. The eleven loci were tested in 24 individuals from two populations in Southwest China. These loci were highly polymorphic with numbers of alleles per locus ranging from 9 to 24 and expected heterozygosities from 0.898 to 0.967. Eight loci followed Hardy–Weinberg expectations after Bonferroni correction for multiple comparisons. No significant linkage association was found among all these loci. The eleven polymorphic microsatellite loci will be useful in determining species boundaries of the South China field mouse.     

Multiple target loci assembly sequencing (mTAS)

Here we present multiple target loci assembly sequencing (mTAS), a method for examining multiple genomic loci in a single DNA sequencing read. The key to the success of mTAS target sequencing is the uniform amplification of multiple target genomic loci into a single DNA fragment using polymerase cycling assembly (PCA). Using this strategy, we successfully collected multiloci sequence information from a single DNA sequencing run. We applied mTAS to examine 29 different sets of human genomic loci, each containing from 2 to 11 single-nucleotide polymorphisms (SNP) present at different exons. We believe mTAS can be used to reduce the cost of Sanger sequencing-based genetic analysis.     

Presence of multiple independent effects in risk loci of common complex human diseases

Many genetic loci and SNPs associated with many common complex human diseases and traits are now identified. The total genetic variance explained by these loci for a trait or disease, however, has often been very small. Much of the "missing heritability" has been revealed to be hidden in the genome among the large number of variants with small effects. Several recent studies have reported the presence of multiple independent SNPs and genetic heterogeneity in trait-associated loci. It is therefore reasonable to speculate that such a phenomenon could be common among loci known to be associated with a complex trait or disease. For testing this hypothesis, a total of 117 loci known to be associated with rheumatoid arthritis (RA), Crohn disease (CD), type 1 diabetes (T1D), or type 2 diabetes (T2D) were selected. The presence of multiple independent effects was assessed in the case-control samples genotyped by the Wellcome Trust Case Control Consortium study and imputed with SNP genotype information from the HapMap Project and the 1000 Genomes Project. Eleven loci with evidence of multiple independent effects were identified in the study, and the number was expected to increase at larger sample sizes and improved statistical power. The variance explained by the multiple effects in a locus was much higher than the variance explained by the single reported SNP effect. The results thus significantly improve our understanding of the allelic structure of these individual disease-associated loci, as well as our knowledge of the general genetic mechanisms of common complex traits and diseases.     

Induction of gene mutations in mice: The multiple endpoint approach

The multiple endpoint mammalian mutagenesis approach developed in our institute screens in the same animal for recessive specific-locus alleles at 7 loci, approximately 30 loci coding for dominant-cataract mutations, 23 loci controlling protein-charge changes and 12 loci for enzyme-activity alterations. Experiments to screen for the approximately 70 loci in the same offspring of treated male mice were performed with ethylnitrosourea (ENU), procarbazine and X-ray exposure. Mutations were recovered for each genetic endpoint in all treatment groups where a sufficient number of offspring was scored. ENU treatment is highly effective in inducing mutations to all genetic endpoints. The mutations were confirmed by breeding tests. The mutation rates to specific-locus and enzyme-activity alleles were both higher than the mutation rates to either dominant-cataract or protein-charge alleles. The advantages and possibilities of the multiple endpoint approach are discussed in detail.     

Coevolution of self-fertilization and inbreeding depression. III. Homozygous lethal mutations at multiple loci.

We study the evolution of the rate of self-fertilization in response to deleterious mutations at multiple loci. Although partial selfing induces associations among loci even in the absence of linkage, associations among mutations at different loci are of a smaller order of magnitude than the mutation rate. Genotypes that carry homozygous lethal mutations in heterozygous form at i loci occur in frequencies of the order (Ti) mu i, in which T denotes the number of viability loci and mu the mutation rate. While associations between mutations at different loci remain small even under inbreeding, each viability locus develops an association with the modifier of the rate of self-fertilization that substantially affects the evolution of the breeding system. Positive associations between enhancers of selfing and haplotypes carrying multiple wild-type alleles and positive associations in heterozygosity between the modifier locus and the viability loci promote evolutionary increases in the rate of self-fertilization.     

Biochemical genetics of mosquitofish III. Incidence and significance of multiple insemination

Evidence for multiple insemination in the mosquitofish, Gambusia affinis, was detected by statistical analysis of the enzyme phenotypes in mother-offspring combinations for three loci. The probability of nondetection of multiple insemination events for diallelic loci with two or more males mating with a female was presented. The frequency of multiple insemination for G. affinis is at least 56% but is probably nearer to 100%. The high incidence of multiple insemination together with sperm storage and delayed fertilization could have important effects on allele frequencies through time and space. These effects are due to the increase in effective population size and overlap among generations.     

Multiplexed single-cell profiling of chromatin states at genomic loci by expansion microscopy

Proper regulation of genome architecture and activity is essential for the development and function of multicellular organisms. Histone modifications, acting in combination, specify these activity states at individual genomic loci. However, the methods used to study these modifications often require either a large number of cells or are limited to targeting one histone mark at a time. Here, we developed a new method called Single Cell Evaluation of Post-TRanslational Epigenetic Encoding (SCEPTRE) that uses Expansion Microscopy (ExM) to visualize and quantify multiple histone modifications at non-repetitive genomic regions in single cells at a spatial resolution of ∼75 nm. Using SCEPTRE, we distinguished multiple histone modifications at a single housekeeping gene, quantified histone modification levels at multiple developmentally-regulated genes in individual cells, and evaluated the relationship between histone modifications and RNA polymerase II loading at individual loci. We find extensive variability in epigenetic states between individual gene loci hidden from current population-averaged measurements. These findings establish SCEPTRE as a new technique for multiplexed detection of combinatorial chromatin states at single genomic loci in single cells.     

Genetic control of hexokinase variation in insects

Hexokinase variation in insects appears to be under the control of a single locus in some species and under multiple-locus control in others. It is often difficult to distinguish the number of loci controlling hexokinase expression. Analysis of hexokinase electrophoretic patterns in six species of mosquitoes and five species of crickets, as well as a review of hexokinase variation in other insect species, is used to emphasize the importance of interspecific comparisons when making genetic inferences. Evidence is provided which adds support for multiple hexokinase loci in dipterans. Hexokinase control by multiple loci may be difficult to determine in some species because of tight linkage, disequilibrium, and/or posttranslational modification.This work was supported by NIH Grant AI 11132-09, NSF Grant DEB-7811188, and Yale Sesil Anonymous Funds.     

Likelihood Analysis of Recombinational Disequilibrium in Multiple-Locus Gametic Frequencies

Likelihood methods are developed for the estimation and testing of multiple-locus gametic disequilibria, using log-linear models of parametric effects. The estimates of disquilibrium are related to Kimura's Z-measure, and may be extended to multiple alleles and multiple loci. Likelihood ratio test criteria are constructed, which are asymptotically distributed as chi(2). The analysis is partitioned into various components corresponding to two-locus, residual three-locus, and higher order disequilibria. A four-locus example from Hordeum vulgare L. is utilized to illustrate the analysis. Most of the multiple-locus disequilibrium is accounted for by two-locus effects, and closely linked loci show considerably more disequilibrium than unlinked loci. It is shown that all possible pairwise comparisons are not statistically independent.     

Coincidence in map positions between pathogen-induced defense-responsive genes and quantitative resistance loci in rice

Quantitative disease resistance conferred by quantitative trait loci (QTLs) is presumably of wider spectrum and durable. Forty-four cDNA clones, representing 44 defense-responsive genes, were fine mapped to 56 loci distributed on 9 of the 12 rice chromosomes. The locations of 32 loci detected by 27 cDNA clones were associated with previously identified resistance QTLs for different rice diseases, including blast, bacterial blight, sheath blight and yellow mottle virus. The loci detected by the same multiple-copy cDNA clones were frequently located on similar locations of different chromosomes. Some of the multiple loci detected by the same clones were all associated with resistance QTLs. These results suggest that some of the genes may be important components in regulation of defense responses against pathogen invasion and they may be the candidates for studying the mechanism of quantitative disease resistance in rice.     

Stochastic search gene suggestion: a Bayesian hierarchical model for gene mapping

Mapping the genes for a complex disease, such as diabetes or rheumatoid arthritis (RA), involves finding multiple genetic loci that may contribute to the onset of the disease. Pairwise testing of the loci leads to the problem of multiple testing. Looking at haplotypes, or linear sets of loci, avoids multiple tests but results in a contingency table with sparse counts, especially when using marker loci with multiple alleles. We propose a hierarchical Bayesian model for case-parent triad data that uses a conditional logistic regression likelihood to model the probability of transmission to a diseased child. We define hierarchical prior distributions on the allele main effects to model the genetic dependencies present in the human leukocyte antigen (HLA) region of chromosome 6. First, we add a hierarchical level for model selection that accounts for both locus and allele selection. This allows us to cast the problem of identifying genetic loci relevant to the disease into a problem of Bayesian variable selection. Second, we attempt to include linkage disequilibrium as a covariance structure in the prior for model coefficients. We evaluate the performance of the procedure with some simulated examples and then apply our procedure to identifying genetic markers in the HLA region that influence risk for RA. Our software is available on the website http://www.epigenetic.org/Linkage/ssgs-public/.     

Coincidence in map positions between pathogen-induced defense-responsive genes and quantitative resistance loci in rice

Quantitative disease resistance conferred by quantitative trait loci (QTLs) is presumably of wider spectrum and durable. Forty-four cDNA clones, representing 44 defense-responsive genes, were fine mapped to 56 loci distributed on 9 of the 12 rice chromosomes. The locations of 32 loci detected by 27 cDNA clones were associated with previously identified resistance QTLs for different rice diseases, including blast, bacterial blight, sheath blight and yellow mottle virus. The loci detected by the same multiple-copy cDNA clones were frequently located on similar locations of different chromosomes. Some of the multiple loci detected by the same clones were all associated with resistance QTLs. These results suggest that some of the genes may be important components in regulation of defense responses against pathogen invasion and they may be the candidates for studying the mechanism of quantitative disease resistance in rice.     

A comprehensive study of multiple mapping and feature selection for correction strategy in the analysis of small RNAs from SOLiD sequencing

High-throughput sequencing is a powerful tool for discovering and profiling microRNAs (miRNAs) to gain further insights into their biogenesis and function. Due to shorter size, short RNAs from deep sequencing dataset are prone to map to multiple loci with an equal number of mismatches, especially among multicopy miRNA precursors and homologous miRNA genes. Systematic analysis of SOLiD sequencing dataset showed that 37.94% short RNAs could simultaneously map to more than one miRNA precursor, and more short RNAs were found to have multiple genomic loci. Improper selection from candidate loci might lose some mapping information, influence miRNA expression profile or even mislead to identify novel miRNAs. A comprehensive study indicated several potential features for correction strategy: location and distribution of mismatches, quality values, expression profiles of multiple isomiRs (miRNA variants), miRNA* and moRs (miRNA-offset-RNAs) at candidate locus and in its flank sequence. Further studies should develop an approach to correct the widespread phenomenon of multiple mapping based on these features, and improve accuracy of profiling and discovering miRNAs.     

Identification of selective sweeps using a dynamically adjusted number of linked microsatellites

There is currently large interest in distinguishing the signatures of genetic variation produced by demographic events from those produced by natural selection. We propose a simple multilocus statistical test to identify candidate sites of selective sweeps with high power. The test is based on the variability profile measured in an array of linked microsatellites. We also show that the analysis of flanking markers drastically reduces the number of false positives among the candidates that are identified in a genomewide survey of unlinked loci and find that this property is maintained in many population-bottleneck scenarios. However, for a certain range of intermediately severe population bottlenecks we find genomic signatures that are very similar to those produced by a selective sweep. While in these worst-case scenarios the power of the proposed test remains high, the false-positive rate reaches values close to 50%. Hence, selective sweeps may be hard to identify even if multiple linked loci are analyzed. Nevertheless, the integration of information from multiple linked loci always leads to a considerable reduction of the false-positive rate compared to a genome scan of unlinked loci. We discuss the application of this test to experimental data from Drosophila melanogaster.     

A Bayesian heterogeneous analysis of variance approach to inferring recent selective sweeps

The distribution of microsatellite allele sizes in populations aids in understanding the genetic diversity of species and the evolutionary history of recent selective sweeps. We propose a heterogeneous Bayesian analysis of variance model for inferring loci involved in recent selective sweeps by analyzing the distribution of allele sizes at multiple loci in multiple populations. Our model is shown to be consistent with a multilocus test statistic, ln RV, proposed for identifying microsatellite loci involved in recent selective sweeps. Our methodology differs in that it accepts original allele size data rather than summary statistics and allows the incorporation of prior knowledge about allele frequencies using a hierarchical prior distribution consisting of log normal and gamma probability distributions. Interesting features of the model are its ability to simultaneously analyze allele size data for any number of populations and to cope with the presence of any number of selected loci. The utility of the method is illustrated by application to two sets of microsatellite allele size data for a group of West African Anopheles gambiae populations. The results are consistent with the suppressed-recombination model of speciation, and additional candidate loci on chromosomes 2 (079 and 175) and 3 (088) are discovered that escaped former analysis.     

High prevalence of multiple paternity within fruits in natural populations of Silene latifolia, as revealed by microsatellite DNA analysis

Data on multiple paternity within broods has been gathered in several animal species, and comparable data in plants would be of great importance to understand the evolution of reproductive traits in a common framework. In this study, we first isolated and characterized six microsatellite loci from the dioecious plant Silene latifolia (Caryophyllaceae). The polymorphism of the loci was assessed in 60 individual females from four different populations. Two of the investigated loci showed a pattern of inheritance consistent with X-linkage. These microsatellite loci were highly polymorphic and therefore useful tools for parentage analysis. We then used four of the markers to determine paternity within naturally pollinated fruits in four European populations. This study revealed widespread multiple paternity in all populations investigated. The minimum number of fathers per fruit varied from one to nine, with population means ranging from 3.4 to 4.9. The number of fathers per fruit was not significantly correlated with offspring sex ratios. High prevalence of multiple paternity within fruits strongly suggest that pollen competition is likely to occur in this species. This may substantially impact male reproductive success and possibly contribute to increase female and offspring fitness, either through postpollination selection or increased genetic diversity. Wide variation in outcrossing rates may be an overlooked aspect of plant mating systems.     

Genetic dissection of drought and heat‐responsive agronomic traits in wheat

High yield and wide adaptation are principal targets of wheat breeding but are hindered by limited knowledge on genetic basis of agronomic traits and abiotic stress tolerances. In this study, 277 wheat accessions were phenotyped across 30 environments with non‐stress, drought‐stressed, heat‐stressed, and drought‐heat‐stressed treatments and were subjected to genome‐wide association study using 395 681 single nucleotide polymorphisms. We detected 295 associated loci including consistent loci for agronomic traits across different treatments and eurytopic loci for multiple abiotic stress tolerances. A total of 22 loci overlapped with quantitative trait loci identified by biparental quantitative trait loci mapping. Six loci were simultaneously associated with agronomic traits and abiotic stress tolerance, four of which fell within selective sweep regions. Selection in Chinese wheat has increased the frequency of superior marker alleles controlling yield‐related traits in the four loci during past decades, which conversely diminished favourable genetic variation controlling abiotic stress tolerance in the same loci; two promising candidate paralogous genes colocalized with such loci, thereby providing potential targets for studying the molecular mechanism of stress tolerance–productivity trade‐off. These results uncovering promising alleles controlling agronomic traits and/or multiple abiotic stress tolerances, providing insights into heritable covariation between yield and abiotic stress tolerance, will accelerate future efforts for wheat improvement.     

Physical molecular maps of wheat chromosomes

In bread wheat, a set of 527 simple sequence repeats (SSRs) were tried on 164 deletion lines, leading to a successful mapping of 270 SSRs on 313 loci covering all 21 chromosomes. A maximum of 119 loci (38%) were located on B subgenome, and a minimum of 90 loci (29%) mapped on D subgenome. Similarly, homoeologous group 7 carried a maximum of 61 loci (19%), and group 4 carried a minimum of 22 loci (7%). Of the cited 270 SSRs, 39 had multiple loci, but only eight of these detected homoeologous loci. Linear order of loci in physical maps largely corresponded with those in the genetic maps. Apparently, distances between each of only 26 pairs of loci significantly differed from the corresponding distances on genetic maps. Some loci, which were genetically mapped close to the centromere, were physically located distally, while other loci that were mapped distally in the genetic maps were located in the proximal bins in the physical maps. This suggested that although the linear order of the loci was largely conserved, variation does exist between genetic and physical distances.Electronic Supplementary Material Supplementary material is available for this article at .