Chromosomal location effects on gene sequence evolution in mammals.

BACKGROUND: Nucleotide substitution rates and G + C content vary considerably among mammalian genes. It has been proposed that the mammalian genome comprises a mosaic of regions - termed isochores - with differing G + C content. The regional variation in gene G + C content might therefore be a reflection of the isochore structure of chromosomes, but the factors influencing the variation of nucleotide substitution rate are still open to question. RESULTS: To examine whether nucleotide substitution rates and gene G + C content are influenced by the chromosomal location of genes, we compared human and murid (mouse or rat) orthologues known to belong to one of the chromosomal (autosomal) segments conserved between these species. Multiple members of gene families were excluded from the dataset. Sets of neighbouring genes were defined as those lying within 1 centiMorgan (cM) of each other on the mouse genetic map. For both synonymous substitution rates and G + C content at silent sites, neighbouring genes were found to be significantly more similar to each other than sets of genes randomly drawn from the dataset. Moreover, we demonstrated that the regional similarities in G + C content (isochores) and synonymous substitution rate were independent of each other. CONCLUSIONS: Our results provide the first substantial statistical evidence for the existence of a regional variation in the synonymous substitution rate within the mammalian genome, indicating that different chromosomal regions evolve at different rates. This regional phenomenon which shapes gene evolution could reflect the existence of 'evolutionary rate units' along the chromosome.     

Clustering of Genes Coding for DNA Binding Proteins in a Regionof Atypical Evolution of the Human Genome

Comparison of the human and mouse genomes has revealed that significant variations in evolutionary rates exist among genomic regions and that a large part of this variation is interchromosomal. We confirm in this work, using a large collection of introns, that human chromosome 19 is the one that shows the highest divergence with respect to mouse. To search for other differences among chromosomes, we examine the distribution of gene functions in human and mouse chromosomes using the Gene Ontology definitions. We found by correspondence analysis that among the strongest clusterings of gene functions in human chromosomes is a group of genes coding for DNA binding proteins in chromosome 19. Interestingly, chromosome 19 also has a very high GC content, a feature that has been proposed to promote an opening of the chromatin, thereby facilitating binding of proteins to the DNA helix. In the mouse genome, however, a similar aggregation of genes coding for DNA binding proteins and high GC content cannot be found. This suggests that the distribution of genes coding for DNA binding proteins and the variations of the chromatin accessibility to these proteins are different in the human and mouse genomes. It is likely that the overall high synonymous and intron rates in chromosome 19 are a by-product of the high GC content of this chromosome.Department of Physiology and Molecular Biodiversity, Institut de Biologia Molecular de Barcelona, CSIC, Jordi Girona 18, 08034 Barcelona, Spain     

The Effects of Mutation and Natural Selection on Codon Bias in the Genes of Drosophila

Codon bias varies widely among the loci of Drosophila melanogaster, and some of this diversity has been explained by variation in the strength of natural selection. A study of correlations between intron and coding region base composition shows that variation in mutation pattern also contributes to codon bias variation. This finding is corroborated by an analysis of variance (ANOVA), which shows a tendency for introns from the same gene to be similar in base composition. The strength of base composition correlations between introns and codon third positions is greater for genes with low codon bias than for genes with high codon bias. This pattern can be explained by an overwhelming effect of natural selection, relative to mutation, in highly biased loci. In particular, this correlation is absent when examining fourfold degenerate sites of highly biased genes. In general, it appears that selection acts more strongly in choosing among fourfold degenerate codons than among twofold degenerate codons. Although the results indicate regional variation in mutational bias, no evidence is found for large scale regions of compositional homogeneity.     

Long-term balancing selection at the west nile virus resistance gene, Oas1b, maintains transspecific polymorphisms in the house mouse

Oligoadenylate synthetases (OASs) are interferon-inducible enzymes that participate in the first line of defense against a wide range of viral infection. Recent studies have determined that Oas1b, a member of the OAS gene family in the house mouse (Mus musculus), provides specific protection against flavivirus infection (e.g., West Nile virus, dengue fever virus, and yellow fever virus). We characterized the nucleotide sequence variation in coding and noncoding regions of the Oas1b gene for a large number of wild-derived strains of M. musculus and related species. Our sequence analyses determined that this gene is one of the most polymorphic genes ever described in any mammal. The level of variation in noncoding regions of Oas1b is an order of magnitude higher than the level reported for other regions of the mouse genome and is significantly different from the level of intraspecific variation expected under neutrality. Furthermore, a phylogenetic analysis of intronic sequences demonstrated that Oas1b alleles are ancient and that their divergence predates several speciation events, resulting in transspecific polymorphisms. The amino acid sequence of Oas1b is also extremely variable, with 1 out of 7 amino acid positions being polymorphic within M. musculus. Oas1b alleles are comparatively more divergent at synonymous positions than most autosomal genes and the ratio of nonsynonymous to synonymous substitution is remarkably high, suggesting that positive selection has been acting on Oas1b. The ancestry of Oas1b polymorphisms and the high level of amino acid polymorphisms strongly suggest that the allelic variation at Oas1b has been maintained in mouse populations by long-term balancing selection.     

Genetic dissection of ethanol tolerance in the budding yeast Saccharomyces cerevisiae

Uncovering genetic control of variation in ethanol tolerance in natural populations of yeast Saccharomyces cerevisiae is essential for understanding the evolution of fermentation, the dominant lifestyle of the species, and for improving efficiency of selection for strains with high ethanol tolerance, a character of great economic value for the brewing and biofuel industries. To date, as many as 251 genes have been predicted to be involved in influencing this character. Candidacy of these genes was determined from a tested phenotypic effect following gene knockout, from an induced change in gene function under an ethanol stress condition, or by mutagenesis. This article represents the first genomics approach for dissecting genetic variation in ethanol tolerance between two yeast strains with a highly divergent trait phenotype. We developed a simple but reliable experimental protocol for scoring the phenotype and a set of STR/SNP markers evenly covering the whole genome. We created a mapping population comprising 319 segregants from crossing the parental strains. On the basis of the data sets, we find that the tolerance trait has a high heritability and that additive genetic variance dominates genetic variation of the trait. Segregation at five QTL detected has explained approximately 50% of phenotypic variation; in particular, the major QTL mapped on yeast chromosome 9 has accounted for a quarter of the phenotypic variation. We integrated the QTL analysis with the predicted candidacy of ethanol resistance genes and found that only a few of these candidates fall in the QTL regions.     

Comparative genomic sequencing reveals a strikingly similar architecture of a conserved syntenic region on human chromosome 11p15.3 (including gene ST5) and mouse chromosome 7

Comparative genomics is a superior way to identify phylogenetically conserved features like genes or regions involved in gene regulation. The comparison of extended orthologous chromosomal regions should also reveal other characteristic traits essential for chromosome or gene function. In the present study we have sequenced and compared a region of conserved synteny from human chromosome 11p15.3 and mouse chromosome 7. In human, this region is known to contain several genes involved in the development of various disorders like Beckwith-Wiedemann overgrowth syndrome and other tumor diseases. Furthermore, in the neighboring chromosome region 11p15.5 extensive imprinting of genes has been reported which might extend to region 11p15.3. The analysis of approximately 730 kb in human and 620 kb in mouse led to the identification of eleven genes. All putative genes found in the mouse DNA were also present in the same order and orientation in the human chromosome. However, in the human DNA one putative gene of unknown function could be identified which is not present in the orthologous position of the mouse chromosome. The sequence similarity between human and mouse is higher in transcribed and exon regions than in non-transcribed segments. Dot plot analysis, however, reveals a surprisingly well-conserved sequence similarity over the entire analyzed region. In particular, the positions of CpG islands, short regions of very high GC content in the 5' region of putative genes, are similar in human and mouse. With respect to base composition, two distinct segments of significantly different GC content exist as well in human as in the mouse. With a GC content of 45% the one segment would correspond to "isochore H1" and the other segment (39% GC in human, 40% GC in mouse) to "isochore L1/L2". The gene density (one gene per 66 kb) is slightly higher than the average calculated for the complete human genome (one gene per 90 kb). The comparison of the number and distribution of repetitive elements shows that the proportion of human DNA made up by interspersed repeats (43.8%) is significantly higher than in the corresponding mouse DNA (30.1%). This partly explains why the human DNA is longer between the landmark genes used to define the orthologous positions in human and mouse.     

An Integrative Genomic Analysis of the Superior Fecundity Phenotype in QSi5 Mice

Laboratory inbred mouse models are a valuable resource to identify quantitative trait loci (QTL) for complex reproductive performance traits. Advances in mouse genomics and high density single nucleotide polymorphism mapping has enabled genome-wide association studies to identify genes linked with specific phenotypes. Gene expression profiles of reproductive tissues also provide potentially useful information for identifying genes that play an important role. We have developed a highly fecund inbred strain, QSi5, with accompanying genotyping for comparative analysis of reproductive performance. Here we analyzed the QSi5 phenotype using a comparative analysis with fecundity data derived from 22 inbred strains of mice from the Mouse Phenome Project, and integration with published expression data from mouse ovary development. Using a haplotype association approach, 400 fecundity-associated regions (FDR < 0.05) with 499 underlying genes were identified. The most significant associations were located on Chromosomes 14, 8, and 6, and the genes underlying these regions were extracted. When these genes were analyzed for expression in an ovarian development profile (GSE6916) several distinctive co-expression patterns across each developmental stage were identified. The genetic analysis also refined 21 fecundity associated intervals on Chromosomes 1, 6, 9, 13, and 17 that overlapped with previously reported reproductive performance QTL. The combined use of phenotypic and in silico data with an integrative genomic analysis provides a powerful tool for elucidating the molecular mechanisms underlying fecundity.     

The molecular genetics of male infertility

Spermatogenesis is an elaborate process involving both cell division and differentiation, and cell-cell interactions. Defects in any of these processes can result in infertility, and in some cases these can be genetic in cause. Mapping experiments have defined at least three regions of the human Y chromosome that are required for normal spermatogenesis. Two of these contain the genes encoding the RNA binding proteins RBM and DAZ, suggesting that the control of RNA metabolism is likely to be an important control point for human spermatogenesis. A similar analysis in mice has shown that at least two regions of the mouse Y chromosome are essential for spermatogenesis. Both genetic and reverse genetic approaches have been used to identify mouse autosomal genes required for spermatogenesis. These studies have shown that genes in a number of different pathways are essential for normal spermatogenesis, and also provide putative models of human infertility.     

Global analysis of regional differences in craniometric diversity and population substructure.

Estimates of genetic diversity in major geographic regions are frequently made by pooling all individuals into regional aggregates. This method can potentially bias results if there are differences in population substructure within regions, since increased variation among local populations could inflate regional diversity. A preferred method of estimating regional diversity is to compute the mean diversity within local populations. Both methods are applied to a global sample of craniometric data consisting of 57 measurements taken on 1734 crania from 18 local populations in six geographic regions: sub-Saharan Africa, Europe, East Asia, Australasia, Polynesia, and the Americas. Each region is represented by three local populations. Both methods for estimating regional diversity show sub-Saharan Africa to have the highest levels of phenotypic variation, consistent with many genetic studies. Polynesia and the Americas both show high levels of regional diversity when regional aggregates are used, but the lowest mean local population diversity. Regional estimates of F(ST) made using quantitative genetic methods show that both Polynesia and the Americas also have the highest levels of differentiation among local populations, which inflates regional diversity. Regional differences in F(ST) are directly related to the geographic dispersion of samples within each region; higher F(ST) values occur when the local populations are geographically dispersed. These results show that geographic sampling can affect results, and suggest caution in making inferences regarding regional diversity when population substructure is ignored.     

Genetic regulation of bone mineral density in mice

Peak bone mass is a major determinant of risk of osteoporotic fracture. Family and twin studies have found a strong genetic component to the determination of bone mineral density (BMD). However, BMD is a complex trait whose expression is confounded by environmental influences and polygenic inheritance. The number, locations and effects of the individual genes contributing to natural variation in this trait are all unknown. The extreme difficulty of dissecting out environmental factors from genetic ones in humans has motivated the investigation of animal models. Genetically distinct animal strains raised under strict environmental control are critical tools for defining genetic regulation. The availability of inbred strains, combined with its relative fecundity, has established the mouse as the best model system for the study of mammalian genetics and physiology. Importantly, genes identified in murine analyses can usually be readily mapped to particular human chromosomal regions because of the high degree of synteny that exists between the mouse and human genomes. We employed quantitative trait locus (QTL) analysis to examine peak BMD in 24 recombinant inbred (RI) mouse strains, derived from a cross between C57BL/6 (B6) and DBA/2 (D2) progenitors (BXD RI). The distribution of BMD values among these strains clearly indicated the presence of strong genetic influences, with an estimated narrow sense heritability of 35%. The differences in peak whole body BMD in the BXD strains were integrated with a large database of genetic markers previously defined in the RI BXD strains to generate chromosome map sites for QTL locations. This QTL analysis provisionally identified a number of chromosomal sites linked to BMD. In the second phase of our BMD QTL mapping efforts, we used three independent mouse populations (all derived from B6 and D2 progenitor strains) to confirm and narrow the genetic locations of 4 QTLs (on chromosomes 1, 2, 4, and 11) that strongly influence the acquisition of peak BMD in mice. Using a novel, fine-mapping approach (recombinant inbred segregation testing), we have succeeded in narrowing two of the BMD-related chromosomal regions and in the process eliminated a number of candidate genes. The homologous regions in the human genome for each of these murine QTLs have been identified in recent human genetic studies. In light of this, we believe that findings in mice should aid in the identification of specific candidate genes for study in humans.     

Heterogeneity in Rates of Recombination across the Mouse Genome

If loci are randomly distributed on a physical map, the density of markers on a genetic map will be inversely proportional to recombination rate. First proposed by MARY LYON, we have used this idea to estimate recombination rates from the Drosophila melanogaster linkage map. These results were compared with results of two other studies that estimated regional recombination rates in D. melanogaster using both physical and genetic maps. The three methods were largely concordant in identifying large-scale genomic patterns of recombination. The marker density method was then applied to the Mus musculus microsatellite linkage map. The distribution of microsatellites provided evidence for heterogeneity in recombination rates. Centromeric regions for several mouse chromosomes had significantly greater numbers of markers than expected, suggesting that recombination rates were lower in these regions. In contrast, most telomeric regions contained significantly fewer markers than expected. This indicates that recombination rates are elevated at the telomeres of many mouse chromosomes and is consistent with a comparison of the genetic and cytogenetic maps in these regions. The density of markers on a genetic map may provide a generally useful way to estimate regional recombination rates in species for which genetic, but not physical, maps are available.     

An essential gene mutagenesis screen across the highly conserved piebald deletion region of mouse chromosome 14

The piebald deletion complex is a set of overlapping chromosomal deficiencies on distal mouse chromosome 14. We surveyed the functional genetic content of the piebald deletion region in an essential gene mutagenesis screen of 952 genomes to recover seven lethal mutants. The ENU‐induced mutations were mapped to define genetic intervals using the piebald deletion panel. Lethal mutations included loci required for establishment of the left‐right embryonic axis and a loss‐of‐function allele of Phr1 resulting in respiratory distress at birth. A functional map of the piebald region integrates experimental genetic data from the deletion panel, mutagenesis screen, and the targeted disruption of specific genes. A comparison of several genomic intervals targeted in regional mutagenesis screens suggests that the piebald region is characterized by a low gene density and high essential gene density with a distinct genomic content and organization that supports complex regulatory interactions and promotes evolutionary stability. genesis 47:392–403, 2009. © 2009 Wiley‐Liss, Inc.     

Chromosomal assignments of the genes for neuroendocrine convertase PC1 (NEC1) to human 5q15-21, neuroendocrine convertase PC2 (NEC2) to human 20p11.1-11.2, and furin (mouse 7[D1-E2] region).

The chromosomal localization of the genes coding for the pro-protein and pro-hormone convertases PC1, PC2, and Furin has been achieved by in situ hybridization. The genes for PC1 and PC2 were located on human chromosomes 5q15-21 and 20p11.1-11.2, respectively. The gene for Furin was assigned to the mouse chromosome 7D1-7E2 region. These data complete the chromosomal localization of these three convertases in both human and mouse. The results confirm the regional correspondence of the human chromosomes 15 and mouse chromosomes 7, as well as between human chromosome 20 and mouse chromosome 2. Furthermore, the identification of the NEC1 locus on human chromosome 5 and mouse chromosome 13 suggests a conservation of synthenic regions between these regions of the human and mouse genomes.     

Characterization by saturation studies of the in vivo binding of [3H]flunitrazepam in mouse brain

The in vivo binding of [3H]flunitrazepam [( 3H]Fln) was characterized in seven regions of the mouse brain. The binding showed saturability and linear Scatchard plots. Hill coefficients were close to unity. Data fitting to a hyperbola by least squares yielded consistent Kd values for all regions studied (0.36-0.6 pmol/mg protein). Bmax values ranged from 0.14 to 0.89 pmol/mg protein, a sixfold regional variation. The order of binding is as follows: cortex greater than hippocampus greater than midbrain = thalamus/hypothalamus greater than striatum much greater than cerebellum greater than brainstem, consistent with that obtained by in vitro binding. The in vivo receptor density and affinity are apparently lower in comparison with in vitro parameters. This is consistent with the observation that the Kd increases and Bmax decreases in vitro when the incubation temperature is increased from 0 degrees C. Non-specific binding has been estimated by displacement of in vivo binding by unlabelled ligand in vitro as well as by pretreatment with unlabelled ligand. The two alternative methods were compared and evaluated. It is concluded that the displacement method provides more reliable estimates of the nonspecific binding. Diazepam-sensitive mice did not differ from the control mice in the in vivo [3H]Fln binding. However, mice pretreated with diazepam 1 or 2 days before have binding reduced by 70 or 30%, respectively. The reduced binding may be explained by receptor occupancy by residual oxazepam. However, the low concentration of the residual oxazepam is an unlikely cause of the phenomenon of "acute tolerance" observed in these mice.     

Accounting for regional niche variation in habitat suitability models

Ecological niche modeling has become an increasingly important tool to address issues in many fields of basic and applied ecology. The ecological niche space occupied across the geographic range, particularly for wide-ranging species, may vary for a variety of evolutionary and non-evolutionary reasons. However, ecological niche models are often applied over large geographic areas without regard for the potential effects of regional variation in adaptation, environmental conditions and their interactions, and species responses, thus significantly reducing their accuracy and utility. We develop regionally partitioned ecological niche models, using GARP, for the wide-ranging North American tree Gleditsia triacanthos (Fabaceae) . Models were constructed based on known tree occurrences at peripheral and range-centre locations, as well as across the geographic range as a whole.
Our results suggest that the niche space occupied by G. triacanthos varies regionally and that between some regions in particular there may be a complete absence of niche overlap. In particular, while there is some overlap between the niche space occupied by trees in the western and central regions of the range, there appears to be virtually no overlap in the niche space occupied by trees in the south of the range and that occupied by western and central trees. This lack of overlap appears to be driven primarily by regional differences in abiotic conditions, rather than regional adaptation per se. The results of our study have several important implications for the future development of habitat suitability models over large geographic areas. Spatial partitioning of data is clearly necessary to improve predictions of models where regional niche variation occurs. For wide-ranging species in particular, regional differences in ecological characteristics may cause apparent niche variation.     

Variation in recombination rate may bias human genetic disease mapping studies

The availability of the human genome sequence and variability information (as from the International HapMap project) will enhance our ability to map genetic disorders and choose targets for therapeutic intervention. However, several factors, such as regional variation in recombination rate, can bias conclusions from genetic mapping studies. Here, we examine the impact of regional variation in recombination rate across the human genome. Through computer simulations and literature surveys, we conclude that genetic disorders have been mapped to regions of low recombination more often than expected if such diseases were randomly distributed across the genome. This concentration in low recombination regions may be an artifact, and disorders appearing to be caused by a few genes of large effect may be polygenic. Future genetic mapping studies should be conscious of this potential complication by noting the regional recombination rate of regions implicated in diseases.     

On numerical uncertainty in evaluation of pest population size

Obtaining reliable estimates of pest insect species abundance is an essential part of ecological monitoring programs. It is often the case that data available for obtaining such estimates are sparse which in turn makes achieving an accurate evaluation difficult. This is especially true for strongly heterogeneous pest population density distributions. In our paper we discuss the accuracy of a mean density estimate when a certain class of high aggregation density distributions is considered and a standard statistical method is employed to handle sparse sampled data. It will be shown in the paper that conventional conclusions about the accuracy of the pest population size evaluation do not work when the data are sparse and a new approach is required. Namely, if the number of traps is small, an estimate of the mean density becomes a random variable with an error of high magnitude and we have to compute the probability of an accurate estimate rather than computing the estimate itself. We have obtained a probability of an accurate estimate based on the assumption that only one trap falls within a sub-domain where the pest population density is different from zero. The probability has been calculated for the one-dimensional and the two-dimensional problem.