Unlinked Genomic Regions Research Articles

Evidence that genomic incompatibilities and other multilocus processes impact hybrid fitness in a rattlesnake hybrid zone.

Hybrid zones provide valuable opportunities to understand the genomic mechanisms that promote speciation by providing insight into factors involved in intermediate stages of speciation. Here, we investigate introgression in a hybrid zone between two rattlesnake species (Crotalus viridis and Crotalus oreganus concolor) that have undergone historical allopatric divergence and recent range expansion and secondary contact. We use Bayesian genomic cline models to characterize genomic patterns of introgression between these lineages and identify loci potentially subject to selection in hybrids. We find evidence for a large number of genomic regions with biased ancestry that deviate from the genomic background in hybrids (i.e., excess ancestry loci), which tend to be associated with genomic regions with higher recombination rates. We also identify suites of excess ancestry loci that show highly correlated allele frequencies (including conspecific and heterospecific combinations) across physically unlinked genomic regions in hybrids. Our findings provide evidence for multiple multilocus evolutionary processes impacting hybrid fitness in this system.

The genetic basis for panicle trait variation in switchgrass (Panicum virgatum)

Key messageWe investigate the genetic basis of panicle architecture in switchgrass in two mapping populations across a latitudinal gradient, and find many stable, repeatable genetic effects and limited genetic interactions with the environment.Grass species exhibit large diversity in panicle architecture influenced by genes, the environment, and their interaction. The genetic study of panicle architecture in perennial grasses is limited. In this study, we evaluate the genetic basis of panicle architecture including panicle length, primary branching number, and secondary branching number in an outcrossed switchgrass QTL population grown across ten field sites in the central USA through multi-environment mixed QTL analysis. We also evaluate genetic effects in a diversity panel of switchgrass grown at three of the ten field sites using genome-wide association (GWAS) and multivariate adaptive shrinkage. Furthermore, we search for candidate genes underlying panicle traits in both of these independent mapping populations. Overall, 18 QTL were detected in the QTL mapping population for the three panicle traits, and 146 unlinked genomic regions in the diversity panel affected one or more panicle trait. Twelve of the QTL exhibited consistent effects (i.e., no QTL by environment interactions or no QTL × E), and most (four of six) of the effects with QTL × E exhibited site-specific effects. Most (59.3%) significant partially linked diversity panel SNPs had significant effects in all panicle traits and all field sites and showed pervasive pleiotropy and limited environment interactions. Panicle QTL co-localized with significant SNPs found using GWAS, providing additional power to distinguish between true and false associations in the diversity panel.

Genetic architecture of adaptive radiation across two trophic levels.

Evolution of trophic diversity is a hallmark of adaptive radiation. Yet, transitions between carnivory and herbivory are rare in young adaptive radiations. Haplochromine cichlid fish of the African Great Lakes are exceptional in this regard. Lake Victoria was colonized by an insectivorous generalist and in less than 20 000 years, several clades of specialized herbivores evolved. Carnivorous versus herbivorous lifestyles in cichlids require many different adaptations in functional morphology, physiology and behaviour. Ecological transitions in either direction thus require many traits to change in a concerted fashion, which could be facilitated if genomic regions underlying these traits were physically linked or pleiotropic. However, linkage/pleiotropy could also constrain evolvability. To investigate components of the genetic architecture of a suite of traits that distinguish invertivores from algae scrapers, we performed quantitative trait locus (QTL) mapping using a second-generation hybrid cross. While we found indications of linkage/pleiotropy within trait complexes, QTLs for distinct traits were distributed across several unlinked genomic regions. Thus, a mixture of independently segregating variation and some pleiotropy may underpin the rapid trophic transitions. We argue that the emergence and maintenance of associations between the different genomic regions underpinning co-adapted traits that evolved and persist against some gene flow required reproductive isolation.

Genetic coupling of life-history and aerobic performance in Atlantic salmon.

A better understanding of the genetic and phenotypic architecture underlying life-history variation is a longstanding aim in biology. Theories suggest energy metabolism determines life-history variation by modulating resource acquisition and allocation trade-offs, but the genetic underpinnings of the relationship and its dependence on ecological conditions have rarely been demonstrated. The strong genetic determination of age-at-maturity by two unlinked genomic regions (vgll3 and six6) makes Atlantic salmon (Salmo salar) an ideal model to address these questions. Using more than 250 juveniles in common garden conditions, we quantified the covariation between metabolic phenotypes—standard and maximum metabolic rates (SMR and MMR), and aerobic scope (AS)—and the life-history genomic regions, and tested if food availability modulates the relationships. We found that the early maturation genotype in vgll3 was associated with higher MMR and consequently AS. Additionally, MMR exhibited physiological epistasis; it was decreased when late maturation genotypes co-occurred in both genomic regions. Contrary to our expectation, the life-history genotypes had no effects on SMR. Furthermore, food availability had no effect on the genetic covariation, suggesting a lack of genotype-by-environment interactions. Our results provide insights on the key organismal processes that link energy use at the juvenile stage to age-at-maturity, indicating potential mechanisms by which metabolism and life-history can coevolve.

Next Generation Sequencing and Bioinformatics Analysis of Family Genetic Inheritance.

Mendelian and complex genetic trait diseases continue to burden and affect society both socially and economically. The lack of effective tests has hampered diagnosis thus, the affected lack proper prognosis. Mendelian diseases are caused by genetic mutations in a singular gene while complex trait diseases are caused by the accumulation of mutations in either linked or unlinked genomic regions. Significant advances have been made in identifying novel diseases associated mutations especially with the introduction of next generation and third generation sequencing. Regardless, some diseases are still without diagnosis as most tests rely on SNP genotyping panels developed from population based genetic analyses. Analysis of family genetic inheritance using whole genomes, whole exomes or a panel of genes has been shown to be effective in identifying disease-causing mutations. In this review, we discuss next generation and third generation sequencing platforms, bioinformatic tools and genetic resources commonly used to analyze family based genomic data with a focus on identifying inherited or novel disease-causing mutations. Additionally, we also highlight the analytical, ethical and regulatory challenges associated with analyzing personal genomes which constitute the data used for family genetic inheritance.

Ecological Divergence with Gene Flow in a Thermophilic Cyanobacterium.

How ecological diversity is maintained and distributed within populations is a longstanding question in microbial ecology. In the thermophilic cyanobacterium Synechococcus B', high observed levels of recombination are predicted to maintain ecological variation despite the simultaneous action of diverse selective pressures on different regions of the genome. To investigate ecological diversity in these bacteria, we directly isolated laboratory strains of Synechococcus B' from samples collected along the thermal gradients of two geothermal environments in Yellowstone National Park. Extensive recombination was evident for a multi-locus sequence data set, and, consequently, our sample did not exhibit the sequence clustering expected for distinct ecotypes evolving by periodic clonal selection. Evidence for local selective sweeps at specific loci suggests that sweeps may be common but that recombination is effective for maintaining diversity of unlinked genomic regions. Thermal performance for strain growth was positively associated with the temperature of the environment, indicating that Synechococcus B' populations consist of locally adapted ecological specialists that occupy specific thermal niches. Because this ecological differentiation is observed despite the absence of dispersal barriers among sites, we conclude that these bacteria may freely exchange much of the genome but that barriers to gene flow exist for loci under direct temperature selection.

High-resolution mapping reveals hundreds of genetic incompatibilities in hybridizing fish species.

Hybridization is increasingly being recognized as a common process in both animal and plant species. Negative epistatic interactions between genes from different parental genomes decrease the fitness of hybrids and can limit gene flow between species. However, little is known about the number and genome-wide distribution of genetic incompatibilities separating species. To detect interacting genes, we perform a high-resolution genome scan for linkage disequilibrium between unlinked genomic regions in naturally occurring hybrid populations of swordtail fish. We estimate that hundreds of pairs of genomic regions contribute to reproductive isolation between these species, despite them being recently diverged. Many of these incompatibilities are likely the result of natural or sexual selection on hybrids, since intrinsic isolation is known to be weak. Patterns of genomic divergence at these regions imply that genetic incompatibilities play a significant role in limiting gene flow even in young species.

Regression-based approach for testing the association between multi-region haplotype configuration and complex trait

BackgroundIt is quite common that the genetic architecture of complex traits involves many genes and their interactions. Therefore, dealing with multiple unlinked genomic regions simultaneously is desirable.ResultsIn this paper we develop a regression-based approach to assess the interactions of haplotypes that belong to different unlinked regions, and we use score statistics to test the null hypothesis of non-genetic association. Additionally, multiple marker combinations at each unlinked region are considered. The multiple tests are settled via the minP approach. The P value of the "best" multi-region multi-marker configuration is corrected via Monte-Carlo simulations. Through simulation studies, we assess the performance of the proposed approach and demonstrate its validity and power in testing for haplotype interaction association.ConclusionOur simulations showed that, for binary trait without covariates, our proposed methods prove to be equal and even more powerful than htr and hapcc which are part of the FAMHAP program. Additionally, our model can be applied to a wider variety of traits and allow adjustment for other covariates. To test the validity, our methods are applied to analyze the association between four unlinked candidate genes and pig meat quality.

Analysis of Allele and Haplotype Diversity Across 25 Genomic Regions in Three Eastern European Populations

Objective: Individual population history is the main reason for the variability of linkage disequilibrium (LD) patterns and haplotype frequencies among populations. Such diversity may influence the transferability of tag SNPs from one population to another. Our goal was to compare patterns of pairwise LD and allele and haplotype frequencies in Estonian and Russian populations, to estimate the genetic variation between populations and assess the potential transferability of tag SNPs. Methods: 452 SNPs from 25 unlinked genomic regions on 12 chromosomes were genotyped in 140 Estonians and 207 Russians from Northern and Western regions of the European area of Russia. Results: The allele frequency distributions were highly consistent between populations (R > 0.90 for all pairwise comparisons). The overall frequency variation among populations was low (F_ST = 0.0054). The number of SNPs with high-range F_ST values (0.02–0.09) was most prominent for the MC5R genomic region. Haplotype heterogeneity among populations was low (F_ST values ranged within 0.00–0.01, with the exception of haploblocks in the ADIPOR2 and MC5R regions). The interpopulation proximity was also evaluated using haplotype diversity. Conclusion: Our data showed a high concordance between the populations studied, which may be considered as the result of their historical formation on a cognate genetic basis.

Equilibrium Processes Cannot Explain High Levels of Short- and Medium-Range Linkage Disequilibrium in the Domesticated Grass Sorghum bicolor

Patterns of linkage disequilibrium (LD) are of interest because they provide evidence of both equilibrium (e.g., mating system or long-term population structure) and nonequilibrium (e.g., demographic or selective) processes, as well as because of their importance in strategies for identifying the genetic basis of complex phenotypes. We report patterns of short and medium range (up to 100 kb) LD in six unlinked genomic regions in the partially selfing domesticated grass, Sorghum bicolor. The extent of allelic associations in S. bicolor, as assessed by pairwise measures of LD, is higher than in maize but lower than in Arabidopsis, in qualitative agreement with expectations based on mating system. Quantitative analyses of the population recombination parameter, rho, however, based on empirical estimates of rates of recombination, mutation, and self-pollination, show that LD is more extensive than expected under a neutral equilibrium model. The disparity between rho and the population mutation parameter, , is similar to that observed in other species whose population history appears to be complex. From a practical standpoint, these results suggest that S. bicolor is well suited for association studies using reasonable numbers of markers, since LD typically extends at least several kilobases but has largely decayed by 15 kb.

Haplotype interaction analysis of unlinked regions

Genetically complex diseases are caused by interacting environmental factors and genes. As a consequence, statistical methods that consider multiple unlinked genomic regions simultaneously are desirable. Such consideration, however, may lead to a vast number of different high-dimensional tests whose appropriate analysis pose a problem. Here, we present a method to analyze case-control studies with multiple SNP data without phase information that considers gene-gene interaction effects while correcting appropriately for multiple testing. In particular, we allow for interactions of haplotypes that belong to different unlinked regions, as haplotype analysis often proves to be more powerful than single marker analysis. In addition, we consider different marker combinations at each unlinked region. The multiple testing issue is settled via the minP approach; the P value of the "best" marker/region configuration is corrected via Monte-Carlo simulations. Thus, we do not explicitly test for a specific pre-defined interaction model, but test for the global hypothesis that none of the considered haplotype interactions shows association with the disease. We carry out a simulation study for case-control data that confirms the validity of our approach. When simulating two-locus disease models, our test proves to be more powerful than association methods that analyze each linked region separately. In addition, when one of the tested regions is not involved in the etiology of the disease, only a small amount of power is lost with interaction analysis as compared to analysis without interaction. We successfully applied our method to a real case-control data set with markers from two genes controlling a common pathway. While classical analysis failed to reach significance, we obtained a significant result even after correction for multiple testing with our proposed haplotype interaction analysis. The method described here has been implemented in FAMHAP.

Inheritance and RAPD tagging of multiple genes for resistance to net blotch in barley

A doubled haploid barley (Hordeum vulgare L.) population that was created from a cross between cultivars 'Léger' and 'CI 9831' was characterized by RAPD (random amplified polymorphic DNA) markers for resistance to isolate WRS857 of Pyrenophora teres Drechs. f. sp. maculata Smedeg., the causal agent of the spot form of net blotch. Resistance, which initially appeared to be conferred by a single gene from the approximate 1:1 (resistant : susceptible) segregation ratio of the doubled-haploid (DH) progeny, was found to be associated with three different genomic regions by RAPD analysis. Of 500 RAPD random primers that were screened against the parents, 195 revealed polymorphic bands, seven showed an association to the resistance in bulks, and these seven markers were mapped to three unlinked genomic regions. Two of these regions, one of which was mapped to chromosome 2, have major resistance genes. The third region has some homology to the chromosome 2 region. This study demonstrates the simultaneous location of markers for more than one gene governing a trait by using RAPD and bulked segregant analysis (BSA).Key words: net blotch, RAPD markers, bulked segregant analysis, barley, doubled haploids.

Inheritance and RAPD tagging of multiple genes for resistance to net blotch in barley

A doubled haploid barley (Hordeum vulgare L.) population that was created from a cross between cultivars 'Léger' and 'CI 9831' was characterized by RAPD (random amplified polymorphic DNA) markers for resistance to isolate WRS857 of Pyrenophora teres Drechs. f. sp. maculata Smedeg., the causal agent of the spot form of net blotch. Resistance, which initially appeared to be conferred by a single gene from the approximate 1:1 (resistant : susceptible) segregation ratio of the doubled-haploid (DH) progeny, was found to be associated with three different genomic regions by RAPD analysis. Of 500 RAPD random primers that were screened against the parents, 195 revealed polymorphic bands, seven showed an association to the resistance in bulks, and these seven markers were mapped to three unlinked genomic regions. Two of these regions, one of which was mapped to chromosome 2, have major resistance genes. The third region has some homology to the chromosome 2 region. This study demonstrates the simultaneous location of markers for more than one gene governing a trait by using RAPD and bulked segregant analysis (BSA).

Divergence Time and Population Size in the Lineage Leading to Modern Humans

We have developed maximum likelihood (ML) methods for comparisons of nucleotide sequences from unlinked genomic regions. In the case of a single species, the ML method primarily estimates the effective population size ( N e) under both constant size and abrupt expansion conditions. In the case of two or three species, the ML method simultaneously estimates the species divergence time and the effective size of ancestral populations. This allows us to trace the evolutionary history of the human population over the past several million years (my). Available sequences at human autosomal loci indicate N e = 10,000 in the Late Pleistocene, a figure concordant with the results obtained from mitochondrial DNA sequence and allele-frequency data analysis, and there is no indication of population expansion. The ML analysis of two species shows that humans diverged from chimpanzees 4.6 my ago and that the human and chimpanzee clade diverged from the gorilla 7.2 my ago. Furthermore, the effective population size of humans more than 4.6 my ago is nearly 10 times larger than N e of modern humans. The effective population size in the human lineage does not seem to have remained constant over the past several million pears. The ML model for three species predicts slightly different, but consistent results to those obtained by the two-species analysis.

Evidence for orthologous seed weight genes in cowpea and mung bean based on RFLP mapping.

A well saturated genomic map is a necessity for a breeding program based on marker assisted selection. To this end, we are developing genomic maps for cowpea (Vigna unguiculata 2N = 22) and mung bean (Vigna radiata 2N = 22) based on restriction fragment length polymorphism (RFLP) markers. Using these maps, we have located major quantitative trait loci (QTLs) for seed weight in both species. Two unlinked genomic regions in cowpea contained QTLs accounting for 52.7% of the variation for seed weight. In mung bean there were four unlinked genomic regions accounting for 49.7% of the variation for seed weight. In both cowpea and mung bean the genomic region with the greatest effect on seed weight spanned the same RFLP markers in the same linkage order. This suggests that the QTLs in this genomic region have remained conserved through evolution. This inference is supported by the observation that a significant interaction (i.e., epistasis) was detected between the QTL(s) in the conserved region and an unlinked RFLP marker locus in both species.

Genetic analysis of tolerance to low-phosphorus stress in maize using restriction fragment length polymorphisms

An understanding of the genetic nature underlying tolerance to low-phosphorus (low-P) stress could aid in the efficient development of tolerant plant strains. The objective of this study was to identify the number of loci in a maize (Zea mays L.) population segregating for tolerance to low-P stress, their approximate location, and the magnitude of their effect.Seventy-seven restriction fragment length polymorphisms (RFLPs) were identified and scored in a maize F2 population derived from a cross between line NY821 and line H99. The F2 individuals were self-pollinated to produce F3 families. Ninety F3 families were grown in a sand-alumina system, which simulated diffusion-limited, low-P soil conditions. The F3 families were evaluated for vegetative growth in a controlled-environment experiment. To identify quantitative trait loci (QTLs) underlying tolerance to low-P stress, the mean phenotypic performances of the F3 families were contrasted based on genotypic classification at each of 77 RFLP marker loci.Six RFLP marker loci were significantly associated with performance under low-P stress (P<0.01). One marker locus accounted for 25% of the total phenotypic variation. Additive gene action was predominant for all of the QTLs identified. Significant marker loci were located on four separate chromosomes representing five unlinked genomic regions. Two marker loci were associated with an additive by additive epistatic interaction. A multiple regression model including three marker loci and the significant epistatic interaction accounted for 46% of the total phenotypic variation. Heterozygosity per se was not predictive of phenotypic performance.