Signatures of genetic variation in human microRNAs point to processes of positive selection and population-specific disease risks.

The occurrence of natural variation in human microRNAs has been the focus of numerous studies during the last 20 years. Most of them have been focused on the role of specific mutations in disease, while a minor proportion seek to analyse microRNA diversity in the genomes of human populations. We analyse the latest human microRNA annotations in the light of the most updated catalogue of genetic variation provided by the 1000 Genomes Project. By means of the in silico analysis of microRNA genetic variation we show that the level of evolutionary constraint of these sequences is governed by the interplay of different factors, like their evolutionary age or genomic location. The role of mutations in the shaping of microRNA-driven regulatory interactions is emphasized with the acknowledgement that, while the whole microRNA sequence is highly conserved, the seed region shows a pattern of higher genetic diversity that appears to be caused by the dramatic frequency shifts of a fraction of human microRNAs. We highlight the participation of these microRNAs in population-specific processes by identifying that not only the seed, but also the loop, are particularly differentiated regions among human populations. The quantitative computational comparison of signatures of population differentiation showed that candidate microRNAs with the largest differences are enriched in variants implicated in gene expression levels (eQTLs), selective sweeps and pathological processes. We explore the implication of these evolutionary-driven microRNAs and their SNPs in human diseases, such as different types of cancer, and discuss their role in population-specific disease risk.

Y-chromosome diversity in Catalan surname samples: insights into surname origin and frequency.

The biological behavior of the Y chromosome, which is paternally inherited, implies that males sharing the same surname may also share a similar Y chromosome. However, socio-cultural factors, such as polyphyletism, non-paternity, adoption, or matrilineal surname transmission, may prevent the joint transmission of the surname and the Y chromosome. By genotyping 17 Y-STRs and 68 SNPs in ~2500 male samples that each carried one of the 50 selected Catalan surnames, we could determine sets of descendants of a common ancestor, the population of origin of the common ancestor, and the date when such a common ancestor lived. Haplotype diversity was positively correlated with surname frequency, that is, rarer surnames showed the strongest signals of coancestry. Introgression rates of Y chromosomes into a surname by non-paternity, adoption, and transmission of the maternal surname were estimated at 1.5-2.6% per generation, with some local variation. Average ages for the founders of the surnames were estimated at ~500 years, suggesting a delay between the origin of surnames (twelfth and thirteenth centuries) and the systematization of their paternal transmission. We have found that, in general, a foreign etymology for a surname does not often result in a non-indigenous origin of surname founders; however, bearers of some surnames with an Arabic etymology show an excess of North African haplotypes. Finally, we estimate that surname prediction from a Y-chromosome haplotype, which may have interesting forensic applications, has a ~60% sensitivity but a 17% false discovery rate.

Identification of risk loci for Crohn's disease phenotypes using a genome-wide association study.

BACKGROUND & AIMS: Crohn's disease is a highly heterogeneous inflammatory bowel disease comprising multiple clinical phenotypes. Genome-wide association studies (GWASs) have associated a large number of loci with disease risk but have not associated any specific genetic variants with clinical phenotypes. We performed a GWAS of clinical phenotypes in Crohn's disease. METHODS: We genotyped 576,818 single-nucleotide polymorphisms in a well-characterized cohort of 1090 Crohn's disease patients of European ancestry. We assessed their association with 17 phenotypes of Crohn's disease (based on disease location, disease behavior, disease course, age at onset, and extraintestinal manifestations). A total of 57 markers with strong associations to Crohn's disease phenotypes (P < 2 × 10(-4)) were subsequently analyzed in an independent replication cohort of 1296 patients of European ancestry. RESULTS: We replicated the association of 4 loci with different Crohn's disease phenotypes. Variants in MAGI1, CLCA2, 2q24.1, and LY75 loci were associated with a complicated stricturing disease course (Pcombined = 2.01 × 10(-8)), disease location (Pcombined = 1.3 × 10(-6)), mild disease course (Pcombined = 5.94 × 10(-7)), and erythema nodosum (Pcombined = 2.27 × 10(-6)), respectively. CONCLUSIONS: In a GWAS, we associated 4 loci with clinical phenotypes of Crohn's disease. These findings indicate a genetic basis for the clinical heterogeneity observed for this inflammatory bowel disease.

A genome-wide association study on a southern European population identifies a new Crohn's disease susceptibility locus at RBX1-EP300.

OBJECTIVE: Genome-wide association studies (GWAS) have identified multiple risk loci for Crohn's disease (CD). However, the cumulative risk exerted by these loci is low, and the likelihood that additional, as-yet undiscovered loci contribute to the risk of CD is very high. We performed a GWAS on a southern European population to identify new CD risk loci. DESIGN: We genotyped 620 901 genome markers on 1341 CD patients and 1518 controls from Spain. The top association signals representing new candidate risk loci were subsequently analysed in an independent replication cohort of 1365 CD patients and 1396 controls. RESULTS: We identified a genome-wide significant association on chromosome 22q13.2 in the intergenic region between the RBX1 and EP300 genes (single nucleotide polymorphism rs4820425, OR 1.27, 95% CI 1.17 to 1.38, p=3.42E-8). We also found suggestive evidence for the association of the IFNGR2 (21q22.11), FOXP2 (7q31), MACROD2 (20p12.1) and AIF1 (6p21.3) loci with CD risk. CONCLUSIONS: In this GWAS performed on a southern European cohort, we have identified a new risk locus for CD between RBX1 and EP300. This study demonstrates that using populations of different ancestry is a useful strategy to identify new risk loci for CD.

Nucleotide substitution rates for the full set of mitochondrial protein-coding genes in Coleoptera.

The ages of cladogenetic events in Coleoptera are frequently estimated with mitochondrial protein-coding genes (MPCGs) and the "standard" mitochondrial nucleotide substitution rate for arthropods. This rate has been used for different mitochondrial gene combinations and time scales despite it was estimated on short mitochondrial sequences from few comparisons of close related species. These shortcomings may cause greater impact at deep phylogenetic levels as errors in rates and ages increase with branch lengths. We use the full set of MPCGs of 15 species of beetles (two of them newly sequenced here) to estimate the nucleotide evolutionary rates in a reconstructed phylogeny among suborders, paying special attention to the effect of data partitioning and model choices on these estimations. The optimal strategy for nucleotide data, as measured with Bayes factors, was partitioning by codon position. This retrieved Adephaga as a sister group to Myxophaga with strong support (expected-likelihood weights test 0.94-1) and both sisters to Polyphaga, in contradiction with the most currently accepted views. The hypothesis of Archostemata being sister to the remaining Coleoptera, which is in agreement with morphology, was increasingly supported when third codon sites were recoded or completely removed, sequences were analyzed as AA, and heterogeneous models were implemented but the support levels remained low. Nucleotide substitution rates were strongly affected by the choice of data partitioning (codon position versus individual genes), with up to sixfold levels of variation, whereas differences in the molecular clock algorithm produced changes of only about 20%. The global mitochondrial protein coding rate using codon partitioning and an estimated age of 250 million years (MY) for the origin of the Coleoptera was 1.34% per branch per MY, which closely matches the 'standard' clock of 1.15% per MY. The estimation of the rates on alternative topologies gave similar results. Using local molecular clocks, the evolutionary rate in the Polyphaga and Archostemata was estimated to be nearly twice as fast as in the Adephaga and Myxophaga (1.03% versus 0.53% per MY). Rates across individual genes varied from 0.55% to 8.61% per MY. Our results suggest that cox1 might not be an optimal gene for implementing molecular clocks in deep phylogenies for beetles because it shows relatively slow rates at first and second codon positions but very fast rates at third ones. In contrast, nad5, nad4 and nad2 perform better, as they exhibit more homogeneous rates among codon positions.

The Mediterranean paradox for susceptibility factors in coronary heart disease extends to genetics.

The incidence of coronary heart disease (CHD) shows a North to South gradient in Europe. We tested whether that gradient could be accounted for by the distribution of putative susceptibility genotypes. We correlated the published frequencies of susceptibility genotypes for the genes most often associated with CHD (ACE, AGT, APOE, F2, F5, MTHFR, PON1, and SERPINE1) with the incidence of the disease, controlling for the effects of smoking, systolic pressure, total cholesterol, and body-mass index. In three polymorphisms a negative correlation between the incidence of CHD and the frequency of a suceptibility genotype was observed. For ACE this correlation was significantly negative even when discounting classical susceptibility factors. This suggests that some alleles described as susceptibility factors cannot account for disease incidence at the population level. A genetic component must be added to the "Mediterranean paradox": genetic variants deemed to be risk factors for CHD show a geographical pattern uncorrelated with the disease incidence. This pattern can be understood from the history of populations which has shaped the genetic diversity of the European populations in North-South clines, similar to what is observed for CHD incidence, which will tend to create spurious correlations with polymorphisms related, or not related, to the disease.

Comparative analysis of cancer genes in the human and chimpanzee genomes.

BACKGROUND: Cancer is a major medical problem in modern societies. However, the incidence of this disease in non-human primates is very low. To study whether genetic differences between human and chimpanzee could contribute to their distinct cancer susceptibility, we have examined in the chimpanzee genome the orthologous genes of a set of 333 human cancer genes. RESULTS: This analysis has revealed that all examined human cancer genes are present in chimpanzee, contain intact open reading frames and show a high degree of conservation between both species. However, detailed analysis of this set of genes has shown some differences in genes of special relevance for human cancer. Thus, the chimpanzee gene encoding p53 contains a Pro residue at codon 72, while this codon is polymorphic in humans and can code for Arg or Pro, generating isoforms with different ability to induce apoptosis or interact with p73. Moreover, sequencing of the BRCA1 gene has shown an 8 Kb deletion in the chimpanzee sequence that prematurely truncates the co-regulated NBR2 gene. CONCLUSION: These data suggest that small differences in cancer genes, as those found in tumor suppressor genes, might influence the differences in cancer susceptibility between human and chimpanzee. Nevertheless, further analysis will be required to determine the exact contribution of the genetic changes identified in this study to the different cancer incidence in non-human primates.

Radiation and phylogeography in the Japanese macaque, Macaca fuscata.

The Japanese macaque (Macaca fuscata) presumably differentiated from eastern rhesus macaque (Macaca mulatta) populations during the Pleistocene and the two species are closely related. In order to analyse speciation and subspeciation events in the Japanese macaque and to describe historical and current relationships among their populations, we sequenced and analysed a fragment of 392bp of mitochondrial DNA (mtDNA) control region in 50 individuals belonging to six populations of Japanese macaque and compared these sequences with 89 eastern rhesus macaque control region sequences from GenBank/EMBL database. There were high genetic similarities between both species and only two positions were fixed within each species, which supports the inclusion of the Japanese macaque in a single species with eastern populations of rhesus macaques. Japanese macaque ancestors colonised Japan after the separation of the two species, estimated at between 0.31 and 0.88 million years ago (Mya). The star-like phylogeny, multimodal mismatch distribution, and lack of correlation between geographic and genetic distances are in accordance with a rapid dispersion of macaques throughout the archipelago after the arrival into Japan. The species shows low genetic variation within populations and high levels of genetic differentiation among populations with no mtDNA haplotype shared across populations. Genetic distances between Yakushima macaques (Macaca fuscata yakui) and any other population of Macaca fuscata fuscata subspecies are comparable to the distances between populations of Honshu, Awajishima, and Kyushu, not supporting the classification of Yakushima macaques as a different subspecies.

Spatial patterns of cystic fibrosis mutation spectra in European populations.

Cystic fibrosis (CF) is the most frequent severe recessive disorder in European populations. We have analyzed its mutation frequency spectrum in 94 European, North African and SW Asian populations taken from the literature. Most major mutations as well as the incidence of CF mutations showed clinals patterns as demonstrated by autocorrelogram analysis. More importantly, measures of mutation diversity did also show clinal patterns, with mutation spectra being more diverse in southern than in northern Europe. This increased diversity would imply roughly a three-fold long-term effective population size in southern than in northern Europe. Distances were computed among populations based on their CF mutation frequencies and compared with distances based on other genic regions. CF-based distances correlated with mtDNA but not with Y-chromosome-based distances, which may be a consequence of the relatively homogeneous CF mutation frequencies in European populations.

Can a place of origin of the main cystic fibrosis mutations be identified?

The genetic background of the mutations that most often cause cystic fibrosis (CF) is different from that of non-CF chromosomes in populations of European origin. It is not known whether these haplotype backgrounds could be found at high frequencies in populations in which CF is, at present, not common; such populations would be candidates for the place of origin of CF mutations. An analysis of haplotypes of CF transmembrane conductance regulator, together with their variation in specific CF chromosomes, in a worldwide survey of normal chromosomes shows (1) a very low frequency or absence of the most common CF haplotypes in all populations analyzed and (2) a strong genetic variability and divergence, among various populations, of the chromosomes that carry disease-causing mutations. The depth of the gene genealogy associated with disease-causing mutations may be greater than that of the evolutionary process that gave rise to present-day human populations. The concept of "population of origin" lacks either spatial or temporal meaning for mutations that are likely to have been present in Europeans before the ethnogenesis of present populations; subsequent population processes may have erased the traces of their geographic origin.