Environmental selection on transcriptome-derived SNPs in a high gene flow marine fish, the Atlantic herring (Clupea harengus)

High gene flow is considered the norm for most marine organisms and is expected to limit their ability to adapt to local environments. Few studies have directly compared the patterns of differentiation at neutral and selected gene loci in marine organisms. We analysed a transcriptome-derived panel of 281 SNPs in Atlantic herring (Clupea harengus), a highly migratory small pelagic fish, for elucidating neutral and selected genetic variation among populations and to identify candidate genes for environmental adaptation. We analysed 607 individuals from 18 spawning locations in the northeast Atlantic, including two temperature clines (5-12 °C) and two salinity clines (5-35‰). By combining genome scan and landscape genetic analyses, four genetically distinct groups of herring were identified: Baltic Sea, Baltic-North Sea transition area, North Sea/British Isles and North Atlantic; notably, samples exhibited divergent clustering patterns for neutral and selected loci. We found statistically strong evidence for divergent selection at 16 outlier loci on a global scale, and significant correlations with temperature and salinity at nine loci. On regional scales, we identified two outlier loci with parallel patterns across temperature clines and five loci associated with temperature in the North Sea/North Atlantic. Likewise, we found seven replicated outliers, of which five were significantly associated with low salinity across both salinity clines. Our results reveal a complex pattern of varying spatial genetic variation among outlier loci, likely reflecting adaptations to local environments. In addition to disclosing the fine scale of local adaptation in a highly vagile species, our data emphasize the need to preserve functionally important biodiversity.

Is there evidence of selection in the dopamine receptor D4 gene in Australian invasive starling populations?

 Although population genetic theory is largely based on the premise that loci under study are selectively neutral, it has been acknowledged that the study of DNA sequence data under the influence of selection can be useful. In some circumstances, these loci show increased population differentiation and gene diversity. Highly polymorphic loci may be especially useful when studying populations having low levels of diversity overall, such as is often the case with threatened or newly established invasive populations. Using common starlings Sturnus vulgaris sampled from invasive Australian populations, we investigated sequence data of the dopamine receptor D4 gene (DRD4), a locus suspected to be under selection for novelty-seeking behaviour in a range of taxa including humans and passerine birds. We hypothesised that such behaviour may be advantageous when species encounter novel environments, such as during invasion. In addition to analyses to detect the presence of selection, we also estimated population differentiation and gene diversity using DRD4 data and compared these estimates to those from microsatellite and mitochondrial DNA sequence data, using the same individuals. We found little evidence for selection on DRD4 in starlings. However, we did find elevated levels of within-population gene diversity when compared to microsatellites and mitochondrial DNA sequence, as well as a greater degree of population differentiation. We suggest that sequence data from putatively nonneutral loci are a useful addition to studies of invasive populations, where low genetic variability is expected.

Selection and validation of reference genes for gene expression studies by reverse transcription quantitative PCR in Xanthomonas citri subsp citri during infection of Citrus sinensis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Positive selection results in frequent reversible amino acid replacements in the G protein gene of human respiratory syncytial virus

Human respiratory syncytial virus (HRSV) is the major cause of lower respiratory tract infections in children under 5 years of age and the elderly, causing annual disease outbreaks during the fall and winter. Multiple lineages of the HRSVA and HRSVB serotypes co-circulate within a single outbreak and display a strongly temporal pattern of genetic variation, with a replacement of dominant genotypes occurring during consecutive years. In the present study we utilized phylogenetic methods to detect and map sites subject to adaptive evolution in the G protein of HRSVA and HRSVB. A total of 29 and 23 amino acid sites were found to be putatively positively selected in HRSVA and HRSVB, respectively. Several of these sites defined genotypes and lineages within genotypes in both groups, and correlated well with epitopes previously described in group A. Remarkably, 18 of these positively selected tended to revert in time to a previous codon state, producing a flipflop phylogenetic pattern. Such frequent evolutionary reversals in HRSV are indicative of a combination of frequent positive selection, reflecting the changing immune status of the human population, and a limited repertoire of functionally viable amino acids at specific amino acid sites.

Selection of endogenous genes for gene expression studies in Eucalyptus under biotic (Puccinia psidii) and abiotic (acibenzolar-S-methyl) stresses using RT-qPCR

Background: Rust caused by Puccinia psidii Winter has been limiting for the establishment of new Eucalyptus plantations, as well as for resprouting of susceptible genetic materials. Identifying host genes involved in defense responses is important to elucidate resistance mechanisms. Reverse transcription-quantitative PCR is the most common method of mRNA quantitation for gene expression analysis. This method generally employs a reference gene as an internal control to normalize results. A good endogenous control transcript shows minimal variation due to experimental conditions. Findings. We analyzed the expression of 13 genes to identify transcripts with minimal variation in leaves of 60-day-old clonal seedlings of two Eucalyptus clones (rust-resistant and susceptible) subjected to biotic (P. psidii) and abiotic (acibenzolar-S-methyl, ASM) stresses. Conclusions. For tissue samples of clones that did not receive any stimulus, a combination of the eEF2 and EglDH genes was the best control for normalization. When pathogen-inoculated and uninoculated plant samples were compared, eEF2 and UBQ together were more appropriate as normalizers. In ASM-treated and untreated leaves of both clones, transcripts of the CYP and elF4B genes combined were the ones with minimal variation. Finally, when comparing expression in both clones for ASM-treated leaves, P. psidii-inoculated leaves, ASM-treated plus P. psidii-inoculated leaves, and their respective controls, the genes with the most stable expression were EgIDH and UBQ. The chitinase gene, which is highly expressed in studies on plant resistance to phytopathogens, was used to confirm variation in gene expression due to the treatments. © 2010 Laia et al; licensee BioMed Central Ltd.

Divergent evolution and purifying selection of the H (FUT1) gene in New World monkeys (Primates, Platyrrhini)

In the present study, the coding region of the H gene was sequenced and analyzed in fourteen genera of New World primates (Alouatta, Aotus, Ateles, Brachyteles, Cacajao, Callicebus, Callithrix, Cebus, Chiropotes, Lagothrix, Leontopithecus, Pithecia, Saguinus, and Saimiri), in order to investigate the evolution of the gene. The analyses revealed that this coding region contains 1,101 nucleotides, with the exception of Brachyteles, the callitrichines (Callithrix, Leontopithecus, and Saguinus) and one species of Callicebus (moloch), in which one codon was deleted. In the primates studied, the high GC content (63%), the nonrandom distribution of codons and the low evolution rate of the gene (0.513 substitutions/site/MA in the order Primates) suggest the action of a purifying type of selective pressure, confirmed by the Z-test. Our analyses did not identify mutations equivalent to those responsible for the H-deficient phenotypes found in humans, nor any other alteration that might explain the lack of expression of the gene in the erythrocytes of Neotropical monkeys. The phylogenetic trees obtained for the H gene and the distance matrix data suggest the occurrence of divergent evolution in the primates.

Multivariate approach in the selection of superior soybean progeny which carry the RR gene

Efficiency in the use of genetic variability, whether existing or created, increases when properly explored and analysed. Incorporation of biotechnology into breeding programs has been the general practice. The challenge for the researcher is the constant development of new and improved cultivars. The aim of this experiment was to select progenies with superior characteristics, whether or not carriers of the RR gene, derived from bi-parental crosses in the soybean, with the help of multivariate techniques. The experiment was carried out in a family-type experimental design, including controls, during the agricultural year 2010/2011 and 2011/2012 in Jaboticabal in the Brazilian State of São Paulo. From the F3 generation, phenotypically superior plants were selected, which were evaluated for the following traits: number of days to flowering; number of days to maturity; height of first pod insertion; plant height at maturity; lodging; agronomic value; number of branches; number of pods per plant; 100-seed weight; number of seeds per plant; grain yield per plant. Given the results, it appears possible to select superior progeny by principal component analysis. Cluster analysis using the K-means method links progeny according to the most important characteristics in each group and identifies, by the Ward method and by means of a dendrogram, the structure of similarity and divergence between selected progeny. Both methods are effective in aiding progeny selection.

Worldwide genetic variation at the 3 ' untranslated region of the HLA-G gene: balancing selection influencing genetic diversity

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Isolated autosomal dominant growth hormone deficiency: stimulating mutant GH-1 gene expression drives GH-1 splice-site selection, cell proliferation, and apoptosis

The majority of mutations that cause isolated GH deficiency type II (IGHD II) affect splicing of GH-1 transcripts and produce a dominant-negative GH isoform lacking exon 3 resulting in a 17.5-kDa isoform, which further leads to disruption of the GH secretory pathway. A clinical variability in the severity of the IGHD II phenotype depending on the GH-1 gene alteration has been reported, and in vitro and transgenic animal data suggest that the onset and severity of the phenotype relates to the proportion of 17.5-kDa produced. The removal of GH in IGHD creates a positive feedback loop driving more GH expression, which may itself increase 17.5-kDa isoform productions from alternate splice sites in the mutated GH-1 allele. In this study, we aimed to test this idea by comparing the impact of stimulated expression by glucocorticoids on the production of different GH isoforms from wild-type (wt) and mutant GH-1 genes, relying on the glucocorticoid regulatory element within intron 1 in the GH-1 gene. AtT-20 cells were transfected with wt-GH or mutated GH-1 variants (5'IVS-3 + 2-bp T->C; 5'IVS-3 + 6 bp T->C; ISEm1: IVS-3 + 28 G->A) known to cause clinical IGHD II of varying severity. Cells were stimulated with 1 and 10 mum dexamethasone (DEX) for 24 h, after which the relative amounts of GH-1 splice variants were determined by semiquantitative and quantitative (TaqMan) RT-PCR. In the absence of DEX, only around 1% wt-GH-1 transcripts were the 17.5-kDa isoform, whereas the three mutant GH-1 variants produced 29, 39, and 78% of the 17.5-kDa isoform. DEX stimulated total GH-1 gene transcription from all constructs. Notably, however, DEX increased the amount of 17.5-kDa GH isoform relative to the 22- and 20-kDa isoforms produced from the mutated GH-1 variants, but not from wt-GH-1. This DEX-induced enhancement of 17.5-kDa GH isoform production, up to 100% in the most severe case, was completely blocked by the addition of RU486. In other studies, we measured cell proliferation rates, annexin V staining, and DNA fragmentation in cells transfected with the same GH-1 constructs. The results showed that that the 5'IVS-3 + 2-bp GH-1 gene mutation had a more severe impact on those measures than the splice site mutations within 5'IVS-3 + 6 bp or ISE +28, in line with the clinical severity observed with these mutations. Our findings that the proportion of 17.5-kDa produced from mutant GH-1 alleles increases with increased drive for gene expression may help to explain the variable onset progression, and severity observed in IGHD II.

PHYLOGENETIC ANALYSES, GENE DOSAGE EFFECT AND SELECTION IN MEMBERS OF THE LATERAL BOUNDARY DOMAIN (LBD) GENES TRANSCRIPTION FACTOR FAMILY

Wood formation is an economically and environmentally important process and has played a significant role in the evolution of terrestrial plants. Despite its significance, the molecular underpinnings of the process are still poorly understood. We have previously shown that four Lateral Boundary Domain (LBD) transcription factors have important roles in the regulation of wood formation with two (LBD1 and LBD4) involved in secondary phloem and ray cell development and two (LBD15 and LBD18) in secondary xylem formation. Here, we used comparative phylogenetic analyses to test potential roles of the four LBD genes in the evolution of woodiness. We studied the copy number and variation in DNA and amino acid sequences of the four LBDs in a wide range of woody and herbaceous plant taxa with fully sequenced and annotated genomes. LBD1 showed the highest gene copy number across the studied species, and LBD1 gene copy number was strongly and significantly correlated with the level of ray seriation. The lianas, cucumber and grape, with multiseriate ray cells showed the highest gene copy number (12 and 11, respectively). Because lianas’ growth habit requires significant twisting and bending, the less lignified ray parenchyma cells likely facilitate stem flexibility and maintenance of xylem conductivity. We further demonstrate conservation of amino acids in the LBD18 protein sequences that are specific to woody taxa. Neutrality tests showed evidence for strong purifying selection on these gene regions across various orders, indicating adaptive convergent evolution of LBD18. Structural modeling demonstrates that the conserved amino acids have a significant impact on the tertiary protein structure and thus are likely of significant functional importance.

Gene expression changes in the coccolithophore Emiliania huxleyi after 500 generations of selection to ocean acidification

Coccolithophores are unicellular marine algae that produce biogenic calcite scales and substantially contribute to marine primary production and carbon export to the deep ocean. Ongoing ocean acidification particularly impairs calcifying organisms, mostly resulting in decreased growth and calcification. Recent studies revealed that the immediate physiological response in the coccolithophore Emiliania huxleyi to ocean acidification may be partially compensated by evolutionary adaptation, yet the underlying molecular mechanisms are currently unknown. Here, we report on the expression levels of 10 candidate genes putatively relevant to pH regulation, carbon transport, calcification and photosynthesis in E. huxleyi populations short-term exposed to ocean acidification conditions after acclimation (physiological response) and after 500 generations of high CO2 adaptation (adaptive response). The physiological response revealed downregulation of candidate genes, well reflecting the concomitant decrease of growth and calcification. In the adaptive response, putative pH regulation and carbon transport genes were up-regulated, matching partial restoration of growth and calcification in high CO2-adapted populations. Adaptation to ocean acidification in E. huxleyi likely involved improved cellular pH regulation, presumably indirectly affecting calcification. Adaptive evolution may thus have the potential to partially restore cellular pH regulatory capacity and thereby mitigate adverse effects of ocean acidification.

Possible ecological risks of transgenic organism release when transgenes affect mating success: Sexual selection and the Trojan gene hypothesis

Widespread interest in producing transgenic organisms is balanced by concern over ecological hazards, such as species extinction if such organisms were to be released into nature. An ecological risk associated with the introduction of a transgenic organism is that the transgene, though rare, can spread in a natural population. An increase in transgene frequency is often assumed to be unlikely because transgenic organisms typically have some viability disadvantage. Reduced viability is assumed to be common because transgenic individuals are best viewed as macromutants that lack any history of selection that could reduce negative fitness effects. However, these arguments ignore the potential advantageous effects of transgenes on some aspect of fitness such as mating success. Here, we examine the risk to a natural population after release of a few transgenic individuals when the transgene trait simultaneously increases transgenic male mating success and lowers the viability of transgenic offspring. We obtained relevant life history data by using the small cyprinodont fish, Japanese medaka (Oryzias latipes) as a model. Our deterministic equations predict that a transgene introduced into a natural population by a small number of transgenic fish will spread as a result of enhanced mating advantage, but the reduced viability of offspring will cause eventual local extinction of both populations. Such risks should be evaluated with each new transgenic animal before release.

Purifying selection and birth-and-death evolution in the ubiquitin gene family

Ubiquitin is a highly conserved protein that is encoded by a multigene family. It is generally believed that this gene family is subject to concerted evolution, which homogenizes the member genes of the family. However, protein homogeneity can be attained also by strong purifying selection. We therefore studied the proportion (pS) of synonymous nucleotide differences between members of the ubiquitin gene family from 28 species of fungi, plants, and animals. The results have shown that pS is generally very high and is often close to the saturation level, although the protein sequence is virtually identical for all ubiquitins from fungi, plants, and animals. A small proportion of species showed a low level of pS values, but these values appeared to be caused by recent gene duplication. It was also found that the number of repeat copies of the gene family varies considerably with species, and some species harbor pseudogenes. These observations suggest that the members of this gene family evolve almost independently by silent nucleotide substitution and are subjected to birth-and-death evolution at the DNA level.

B lymphocyte involvement in ankylosing spondylitis: the heavy chain variable segment gene repertoire of B lymphocytes from germinal center-like foci in the synovial membrane indicates antigen selection

The synovial membrane (SM) of affected joints in ankylosing spondylitis (AS) is infiltrated by germinal center-like aggregates (foci) of lymphocytes similar to rheumatoid arthritis (RA). We characterized the rearranged heavy chain variable segment (VH) genes in the SM for gene usage and the mutational pattern to elucidate the B lymphocyte involvement in AS.