Genome scans and elusive candidate genes: detecting the variation that matters for speciation

Next‐generation sequencing is providing us with vast amounts of genetic data, yet we are currently struggling in our attempts to make sense of them. In particular, it has proven difficult to link phenotypic divergence and speciation to genome level divergence. In the current issue of Molecular Ecology, Ruegg et al. ( 2014 ) present new empirical results from two closely related bird taxa. They use a promising approach combining genome scan and candidate gene analysis. Their results suggest that we may have been looking in vain for candidate speciation genes by focusing only on genes found within genomic islands of divergence. This is because genes important in divergence and speciation may not be detected by genome scans and because features of the genomic architecture per se may have a large effect on the pattern of genome divergence.     

Genome scans of variation and adaptive change: extended analysis of a candidate locus close to the phantom gene region in Drosophila melanogaster

Nucleotide variation in populations originating from the recent range expansion of a species should reflect their adaptation to new habitats as well as their demographic history. A survey of nucleotide variation at 109 noncoding X-chromosome fragments in a European population of Drosophila melanogaster allowed identifying some candidates to have been recently affected by positive selection. Adaptive changes leave a spatial differential footprint that can be used to discriminate among candidates by extending their study to neighboring regions. Here, we surveyed variation at an approximately 190-kb region spanning a locus exhibiting a significantly skewed frequency spectrum. A stretch of approximately 12 kb with reduced variation was detected within a continuously sequenced region that included the focal fragment. Moreover, the regions flanking this stretch exhibited an excess of high-frequency derived variants. Application of maximum likelihood ratio and goodness-of-fit tests suggested that the pattern of variation detected at the studied region (at cytological bands 17C-17D) might have been shaped by a recent selective change, most probably at or around the phantom gene that encodes CYP306A1, a cytochrome P450 enzyme in the ecdysteroidogenic pathway.     

Genome scans reveal candidate regions involved in the adaptation to host plant in the pea aphid complex

A major goal in evolutionary biology is to uncover the genetic basis of adaptation. Divergent selection exerted on ecological traits may result in adaptive population differentiation and reproductive isolation and affect differentially the level of genetic divergence along the genome. Genome‐wide scan of large sets of individuals from multiple populations is a powerful approach to identify loci or genomic regions under ecologically divergent selection. Here, we focused on the pea aphid, a species complex of divergent host races, to explore the organization of the genomic divergence associated with host plant adaptation and ecological speciation. We analysed 390 microsatellite markers located at variable distances from predicted genes in replicate samples of sympatric populations of the pea aphid collected on alfalfa, red clover and pea, which correspond to three common host‐adapted races reported in this species complex. Using a method that accounts for the hierarchical structure of our data set, we found a set of 11 outlier loci that show higher genetic differentiation between host races than expected under the null hypothesis of neutral evolution. Two of the outliers are close to olfactory receptor genes and three other nearby genes encoding salivary proteins. The remaining outliers are located in regions with genes of unknown functions, or which functions are unlikely to be involved in interactions with the host plant. This study reveals genetic signatures of divergent selection across the genome and provides an inventory of candidate genes responsible for plant specialization in the pea aphid, thereby setting the stage for future functional studies.     

Identification of candidate disease genes in patients with common variable immunodeficiency

Background: Common variable immunodeficiency (CVID), the most prevalent form of primary immunodeficiency (PID), is characterized by hypogammaglobulinemia and recurrent infections. Understanding protein-protein interaction (PPI) networks of CVID genes and identifying candidate CVID genes are critical steps in facilitating the early diagnosis of CVID. Here, the aim was to investigate PPI networks of CVID genes and identify candidate CVID genes using computation techniques. Methods: Network density and biological distance were used to study PPI data for CVID and PID genes obtained from the STRING database. Gene expression data of patients with CVID were obtained from the Gene Expression Omnibus, and then Pearson’s correlation coefficient, a PPI database, and Kyoto Encyclopedia of Genes and Genomes were used to identify candidate CVID genes. We then evaluated our predictions and identified differentially expressed CVID genes. Results: The majority of CVID genes are characterized by a high network density and small biological distance, whereas most PID genes are characterized by a low network density and large biological distance, indicating that CVID genes are more functionally similar to each other and closely interact with one other compared with PID genes. Subsequently, we identified 172 CVID candidate genes that have similar biological functions to known CVID genes, and eight genes were recently reported as CVID-related genes. MYC, a candidate gene, was down-regulated in CVID duodenal biopsies, but up-regulated in blood samples compared with levels in healthy controls. Conclusion: Our findings will aid in a better understanding of the complex of CVID genes, possibly further facilitating the early diagnosis of CVID.     

Genome scans and gene expression microarrays converge to identify gene regulatory loci relevant in schizophrenia

Multiple linkage regions have been reported in schizophrenia, and some appear to harbor susceptibility genes that are differentially expressed in postmortem brain tissue derived from unrelated individuals. We combined traditional genome-wide linkage analysis in a multiplex family with lymphocytic genome-wide expression analysis. A genome scan suggested linkage to a chromosome 4q marker (D4S1530, LOD 2.17, θ=0) using a dominant model. Haplotype analysis using flanking microsatellite markers delineated a 14 Mb region that cosegregated with all those affected. Subsequent genome-wide scan with SNP genotypes supported the evidence of linkage to 4q33–35.1 (LOD=2.39) using a dominant model. Genome-wide microarray analysis of five affected and five unaffected family members identified two differentially expressed genes within the haplotype AGA and GALNT7 (aspartylglucosaminidase and UDP-N-acetyl-alpha-D-galactosamine: polypeptide N-acetylgalactosaminyltransferase 7) with nominal significance; however, these genes did not remain significant following analysis of covariance. We carried out genome-wide linkage analyses between the quantitative expression phenotype and genetic markers. AGA expression levels showed suggestive linkage to multiple markers in the haplotype (maximum LOD=2.37) but to no other genomic region. GALNT7 expression levels showed linkage to regulatory loci at 4q28.1 (maximum LOD=3.15) and in the haplotype region at 4q33–35.1 (maximum LOD=2.37). ADH1B (alcohol dehydrogenase IB) was linked to loci at 4q21–q23 (maximum LOD=3.08) and haplotype region at 4q33–35.1 (maximum LOD=2.27). Seven differentially expressed genes were validated with RT-PCR. Three genes in the 4q33–35.1 haplotype region were also differentially expressed in schizophrenia in postmortem dorsolateral prefrontal cortex: AGA, HMGB2, and SCRG1. These results indicate that combining differential gene expression with linkage analysis may help in identifying candidate genes and potential regulatory sites. Moreover, they also replicate recent findings of complex trans- and cis- regulation of genes.     

Infant and adult human B cell responses to rotavirus share common immunodominant variable gene repertoires

Ab repertoires exhibit marked restrictions during fetal life characterized by biases of variable gene usage and lack of junctional diversity. We tested the hypothesis that Ab repertoire restriction contributes to the observed poor quality of specific Ab responses made by infants to viral infections. We analyzed the molecular determinants of B cell responses in humans to two Ags of rotavirus (RV), a common and clinically important infection of human infants. We sequenced Ab H and L chain V region genes (V(H) and V(L)) of clones expanded from single B cells responding to RV virus protein 6 or virus protein 7. We found that adults exhibited a distinct bias in use of gene segments in the V(H)1 and V(H)4 families, for example, V(H)1-46, V(H)4-31, and V(H)4-61. This gene segment bias differed markedly from the V(H)3 dominant bias seen in randomly selected adult B cells. Recombinant Abs incorporating any of those three immunodominant V(H) segments bound to RV-infected cells and also to purified RV particles. The RV-specific B cell repertoires of infants aged 2-11 mo and those of adults were highly related when compared by V(H), D, J(H), V(L), and J(L) segment selection, extent of junctional diversity, and mean H chain complementarity determining region 3 length. These data suggest that residual fetal bias of the B cell repertoire is not a limiting determinant of the quality of Ab responses to viruses of infants beyond the neonatal period.     

Development of candidate gene markers associated to common bacterial blight resistance in common bean

Common bacterial blight (CBB), caused by Xanthomonas axonopodis pv. phaseoli (Xap), is a major yield-limiting factor of common bean (Phaseolus vulgaris L.) production around the world. Two major CBB-resistant quantitative trait loci (QTL), linked to the sequence characterized amplified region markers BC420 and SU91, are located at chromosomes 6 and 8, respectively. Using map-based cloning approach, four bacterial artificial chromosome (BAC) clones from the BC420-QTL locus and one BAC clone containing SU91 were sequenced by Roche 454 technique and subsequently assembled using merged assemblies from three different programs. Based on the quality of the assembly, only the sequences of BAC 32H6 and 4K7 were used for candidate gene marker (CGM) development and candidate gene (CG) selection. For the BC420-QTL locus, 21 novel genes were predicted in silico by FGENESH using Medicago gene model, whereas 16 genes were identified in the SU91-QTL locus. For each putative gene, one or more primer pairs were designed and tested in the contrasting near isogenic lines. Overall, six and nine polymorphic markers were found in the SU91- and BC420-QTL loci, respectively. Afterwards, association mapping was conducted in a breeding population of 395 dry bean lines to discover marker-trait associations. Two CGMs per each locus showed better association with CBB resistance than the BC420 and SU91 markers, which include BC420-CG10B and BC420-CG14 for BC420_QTL locus, and SU91-CG10 and SU91-CG11 for SU91_QTL locus. The strong associations between CBB resistance and the CGs 10 and 14 from BC420_QTL locus and the CGs 10 and 11 from SU91_QTL locus indicate that the genes 10 and 14 from the BC420 locus are potential CGs underlying the BC420_QTL locus, whereas the genes 10 and 11 from the SU91 locus are potential CGs underlying the SU91_QTL locus. The superiority of SU91-CG11 was further validated in a recombinant inbred line population Sanilac?×?OAC 09-3. Thus, co-dominant CGMs, BC420-CG14 and SU91-CG11, are recommended to replace BC420 and SU91 for marker-assisted selection of common bean with resistance to CBB.     

Genome scans on experimentally evolved populations reveal candidate regions for adaptation to plant resistance in the potato cyst nematode Globodera pallida

Improving resistance durability involves to be able to predict the adaptation speed of pathogen populations. Identifying the genetic bases of pathogen adaptation to plant resistances is a useful step to better understand and anticipate this phenomenon. Globodera pallida is a major pest of potato crop for which a resistance QTL, GpaV_vrn, has been identified in Solanum vernei. However, its durability is threatened as G. pallida populations are able to adapt to the resistance in few generations. The aim of this study was to investigate the genomic regions involved in the resistance breakdown by coupling experimental evolution and high‐density genome scan. We performed a whole‐genome resequencing of pools of individuals (Pool‐Seq) belonging to G. pallida lineages derived from two independent populations having experimentally evolved on susceptible and resistant potato cultivars. About 1.6 million SNPs were used to perform the genome scan using a recent model testing for adaptive differentiation and association to population‐specific covariables. We identified 275 outliers and 31 of them, which also showed a significant reduction in diversity in adapted lineages, were investigated for their genic environment. Some candidate genomic regions contained genes putatively encoding effectors and were enriched in SPRYSECs, known in cyst nematodes to be involved in pathogenicity and in (a)virulence. Validated candidate SNPs will provide a useful molecular tool to follow frequencies of virulence alleles in natural G. pallida populations and define efficient strategies of use of potato resistances maximizing their durability.     

Coexistence of trichome variation in a natural plant population: a combined study using ecological and candidate gene approaches

The coexistence of distinct phenotypes within populations has long been investigated in evolutionary ecology. Recent studies have identified the genetic basis of distinct phenotypes, but it is poorly understood how the variation in candidate loci is maintained in natural environments. In this study, we examined fitness consequences and genetic basis of variation in trichome production in a natural population of Arabidopsis halleri subsp. gemmifera. Half of the individuals in the study population produced trichomes while the other half were glabrous, and the leaf beetle Phaedon brassicae imposed intensive damage to both phenotypes. The fitness of hairy and glabrous plants showed no significant differences in the field during two years. A similar result was obtained when sibling hairy and glabrous plants were transplanted at the same field site, whereas a fitness cost of trichome production was detected under a weak herbivory condition. Thus, equivalent fitness of hairy and glabrous plants under natural herbivory allows their coexistence in the contemporary population. The pattern of polymorphism of the candidate trichome gene GLABROUS1 (GL1) showed no evidence of long-term maintenance of trichome variation within the population. Although balancing selection under fluctuating biotic environments is often proposed to explain the maintenance of defense variation, the lack of clear evidence of balancing selection in the study population suggests that other factors such as gene flow and neutral process may have played relatively large roles in shaping trichome variation at least for the single population level.     

Patterns of genetic diversity in the Andean gene pool of common bean reveal a candidate domestication gene

Most studies on the genetic diversity of common bean (Phaseolus vulgaris L.) have focussed on accessions from the Mesoamerican gene pool compared to the Andean gene pool. A deeper knowledge of the genetic structure of Argentinian germplasm would enable researchers to determine how the Andean domestication event affected patterns of genetic diversity in domesticated beans and to identify candidates for genes targeted by selection during the evolution of the cultivated common bean. A collection of 116 wild and domesticated accessions representing the diversity of the Andean bean in Argentina was genotyped by means of 114 simple sequence repeat (SSR) markers. Forty-seven Mesoamerican bean accessions and 16 Andean bean accessions representing the diversity of Andean landraces and wild accessions were also included. Using the Bayesian algorithm implemented in the software STRUCTURE we identified five major groups that correspond to Mesoamerican and Argentinian wild accessions and landraces and a group that corresponds to accessions from different Andean and Mesoamerican countries. The neighbour-joining algorithm and principal coordinate clustering analysis confirmed the genetic relationships among accessions observed with the STRUCTURE analysis. Argentinian accessions showed a substantial genetic variation with a considerable number of unique haplotypes and private alleles, suggesting that they may have played an important role in the evolution of the species. The results of statistical analyses aimed at identifying genomic regions with consistent patterns of variation were significant for 35 loci (~20 % of the SSRs used in the Argentinian accessions). One of these loci mapped in or near the genomic region of the glutamate decarboxylase gene. Our data characterize the population structure of the Argentinian germplasm. This information on its diversity will be very valuable for use in introgressing Argentinian genes into commercial varieties because the majority of present-day common bean varieties are of Andean origin.     

Natural products drug discovery: accelerating the clinical candidate development using reverse pharmacology approaches

The pharmaceutical industry is facing serious challenges as the drug discovery process is becoming extremely expensive, riskier and critically inefficient. A significant shift from single to multi targeted drugs especially for polygenic syndromes is being witnessed. Strategic options based on natural product drug discovery, ethnopharmacology and traditional medicines are re-emerging to offer good base as an attractive discovery engine. Approaches based on reverse pharmacology may offer efficient development platforms for herbal formulations. Relevant case studies from India and other countries where such approaches have expedited the drug discovery and development process by reducing time and economizing investments with better safety are discussed.     

Abnormalities of lymphokine gene expression in patients with common variable immunodeficiency

Common variable immunodeficiency (CVI) is a syndrome characterized by hypogammaglobulinemia, recurrent bacterial infections, and increased occurrence of both autoimmune disease and malignancy. In our study we examine the expression of lymphokine genes in mitogen-activated T cells from four patients with CVI. T cells from patients with CVI did not differ significantly from normals in total T cell number, CD4/CD8 ratio, CD45R expression, or proliferation in response to PHA. However, T cells from this group of patients did exhibit significant abnormalities of mitogen-induced lymphokine gene expression. T cells from patients exhibited significantly decreased expression of IL-2, IL-4, IL-5, and IFN-gamma when compared to normal controls. In contrast to these abnormal findings, mitogen-activated T cells from patients with CVI expressed normal amounts of IL-2R alpha and c-myc suggesting that these patients have a selective abnormality of T cell activation. Furthermore, it is likely that the deficient production of IFN-gamma by patient T cells is partially due to the abnormality of IL-2 production as the levels of IFN-gamma mRNA detected during the initial IL-2-independent phase of T cell activation were normal and the addition of exogenous rIL-2 was able to normalize IFN-gamma production by PHA-stimulated patient cells. Finally, supernatants from PHA-activated cultures of patients PBMC were deficient in their ability to support Ig secretion by Staphylococcus A Cowan's-activated normal B cells suggesting that these T cell abnormalities may contribute to the pathogenesis of this syndrome.     

Genome scans for the contemporary response to selection in quantitative traits

Genome scans have been an important approach for discovering historical signatures of selection in both model and nonmodel species. An exciting new experimental design for genome scans is to measure the change in allele frequency before and after contemporary selection within a generation, from a single population. The most widely‐used methods, however, have two major limitations: they are based on testing one locus at a time, and they only have power to uncover loci that have evolved under relatively strong selection. On the other hand, complex quantitative traits are common in nature and are caused by several loci of small effect. Selection on a quantitative trait at the phenotypic level is predicted to be accompanied by subtle allele frequency changes in many loci that covary (a polygenic soft sweep), rather than a large, single‐effect allele (a selective sweep). In this issue of Molecular Ecology, Bourret et al. (2014) measure the contemporary response to natural selection across the genome in multiple cohorts of Atlantic salmon during their first year at sea. They introduce a multilocus framework based on groups of markers that covary in their genotypic distribution. While the traditional, single‐locus approach did not find evidence for repeated patterns of selection, the multivariate approach found that a group of covarying SNPs was selected for in different cohorts at one site. Their multilocus framework has potential to be a more fruitful approach for uncovering the genomic basis of adaptation in quantitative traits, although caution should be applied as the framework has yet to be validated with simulated data.     

Integrating candidate gene and quantitative genetic approaches to understand variation in timing of breeding in wild tit populations

Two commonly used techniques for estimating the effect of genes on traits in wild populations are the candidate gene approach and quantitative genetic analyses. However, whether these two approaches measure the same underlying processes remains unresolved. Here, we use these two methods to test whether they are alternative or complementary approaches to understanding genetic variation in the timing of reproduction - a key trait involved in adaptation to climate change - in wild tit populations. Our analyses of the candidate gene Clock show weak correlates with timing variables in blue tits, but no association in great tits, confirming earlier results. Quantitative genetic analyses revealed very low levels of both direct (female) and indirect (male) additive genetic variation in timing traits for both species, in contrast to previous studies on these traits, and much lower than generally assumed. Hence, neither method suggests strong genetic effects on the timing of breeding in birds, and further work should seek to assess the generality of these conclusions. We discuss how differences in the genetic control of traits, species life-history and confounding environmental variables may determine how useful integrating these two techniques is to understand the phenotypic variation in wild populations.     

Genome wide association mapping and candidate gene analysis for pod shatter resistance in Brassica juncea and its progenitor species

Molecular Biology Reports - We investigated phenotypic variations for pod shattering, pod length and number of seeds per pod in large germplasm collections of Brassica juncea (2n?=?36;...     

Combining gene expression and genetic analyses to identify candidate genes involved in cold responses in pea

Cold stress affects plant growth and development. In order to better understand the responses to cold (chilling or freezing tolerance), we used two contrasted pea lines. Following a chilling period, the Champagne line becomes tolerant to frost whereas the Terese line remains sensitive.     

Genome size variation in the common frog Rana temporaria

Genome size variation in the common frog (Rana temporaria) was investigated with flow cytometry in three latitudinally separated populations in Sweden to see whether it could provide a useful tool for sex-identification in this species. Depending on the sex and population, per cell DNA content (2C value) varied from 8.823 to 11.266 pg with a mean (+/- SE) 2C value of 9.961+/-0.083 pg. Analysis of variance revealed significant differences in genome size among populations and between sexes. Females had ca 3% larger genomes (x=10.133+/-0.068 pg) than males (x=9.832+/-0.068 pg) in all of the populations (sex x population interaction: P>0.10). Individuals from the southern-most population had significantly (x=9.330+/-0.081 pg) smaller genomes than those from the more northern populations (x=10.032+/-0.085 and x=10.584+/-0.085 pg, respectively). These results are in line with the interpretation that males in the common frog are the heterogametic sex, and that there exists large (up to 12%) geographic variation in genome size in this species. However, the sex differences in the genome size are too small to be useful in individual sex identification.     

Genome scans of facial features in East Africans and cross-population comparisons reveal novel associations

Facial morphology is highly variable, both within and among human populations, and a sizable portion of this variation is attributable to genetics. Previous genome scans have revealed more than 100 genetic loci associated with different aspects of normal-range facial variation. Most of these loci have been detected in Europeans, with few studies focusing on other ancestral groups. Consequently, the degree to which facial traits share a common genetic basis across diverse sets of humans remains largely unknown. We therefore investigated the genetic basis of facial morphology in an East African cohort. We applied an open-ended data-driven phenotyping approach to a sample of 2,595 3D facial images collected on Tanzanian children. This approach segments the face into hierarchically arranged, multivariate features that capture the shape variation after adjusting for age, sex, height, weight, facial size and population stratification. Genome scans of these multivariate shape phenotypes revealed significant (p < 2.5 × 10⁻⁸) signals at 20 loci, which were enriched for active chromatin elements in human cranial neural crest cells and embryonic craniofacial tissue, consistent with an early developmental origin of the facial variation. Two of these associations were in highly conserved regions showing craniofacial-specific enhancer activity during embryological development (5q31.1 and 12q21.31). Six of the 20 loci surpassed a stricter threshold accounting for multiple phenotypes with study-wide significance (p < 6.25 × 10⁻¹⁰). Cross-population comparisons indicated 10 association signals were shared with Europeans (seven sharing the same associated SNP), and facilitated fine-mapping of causal variants at previously reported loci. Taken together, these results may point to both shared and population-specific components to the genetic architecture of facial variation.