The X chromosome in population genetics

Investigation of 12 X-STR loci in Mongolian and Eastern Han populations of China with comparison to other populations

Due to the unique inheritance pattern, X-chromosomal short tandem repeats (X-STRs) have several advantages in complex kinship cases, such as deficiency cases or grandparent-grandchild and half-sisters testing. In our study, 541 unrelated individuals gathered from Mongolian and Eastern Chinese Han populations were successfully genotyped using the Investigator Argus X-12 kit. We calculated allele/haplotype frequencies and other forensic parameters of the two populations and further explored their genetic distance with already published Chinese populations and six global populations. Our results showed that the 12 X-STR markers were highly informative in the two populations when compared with nine other Chinese populations: significant differences were found at several loci. Geographically neighboring populations or different ethnic groups within the same area appeared to have closer evolutionary relationships. We also analyzed population genetic structure by performing clustering with the STRUCTURE program and Principal Coordinate Analysis (PCoA), and we found that the Chinese and other populations enrolled in this study could be distinguished. Furthermore, Mongolian males were distinguishable from the other studied males by a moderate genetic distance. Our study also expanded the X-STR database, which could facilitate the appropriate application of the 12 X-STR markers in the forensic field in China.

Trauma-Induced Acute X Chromosome Skewing in White Blood Cells Represents an Immuno-Modulatory Mechanism Unique to Females and a Likely Contributor to Sex-Based Outcome Differences

Sex-related outcome disparities following severe trauma have been demonstrated in human and animal studies; however, sex hormone status could not fully account for the differences. This study tested whether X-linked cellular mosaicism, which is unique to females, could represent a genetically based mechanism contributing to sex-related immuno-modulation following trauma. Serial blood samples collected for routine laboratory tests were analyzed for ChrX inactivation (XCI) ratios in white blood cells. Thirty-nine severely injured (mean ISS 19) female trauma patients on mixed racial and ethnic background were tested for initial (baseline) and trauma-induced changes in XCI ratios and their associations with severity of injury and clinical outcome. At admission, two-thirds of the patients showed XCI-ratio values between one and three, about a third presented skewed XCI ratios (3-7 range) and three patients displayed extremely skewed XCI ratios (8-30 range). Serial blood samples during the clinical course showed additional changes in XCI ratios ranging between 20% and 900% over initial. Increasing XCI ratios during the injury course correlated with the severity of trauma, subsequent need for ventilator support and pneumonia. In contrast, initial XCI ratios did not show correlations with injury severity or clinical complications. Initial XCI ratios showed a positive correlation with age but older patients retained the ability to mount trauma-induced secondary XCI changes. These data show that trauma results in X-linked cell selection in females, which is likely to be driven by polymorphic differences between the parental ChrXs. X-linked white blood cell skewing correlates with injury severity and a complicated postinjury clinical course. Female X-linked cellular mosaicism and its capacity to change dynamically during the injury course compared with the lack of this machinery in males may represent a novel immuno-modulatory mechanism contributing to sex-based outcome differences after injury and infection.

Dinucleotide (CA)n tandem repeats on the human X-chromosome and the history of the Mediterranean populations

BACKGROUND

Tandem repeats (STRs) are genomic markers of particular interest in forensic and population genetics. Most of the population data currently available correspond to the variation of STRs of forensic panels, which barely include dinucleotide tandem repeats.

AIMS

The aim of the study is to test the usefulness of a battery of dinucleotide STRs on the X chromosome for population and forensic studies.

SUBJECTS AND METHODS

A total of 672 individuals from 12 Mediterranean populations and two external references were analysed for 15 X-STR following the instructions of the commercial company and using control DNA from the CEPH centre whose sequences are published in GenBank. Genotypic results were analysed using standard population genetics methods including estimates of linkage disequilibrium, population structure and gene flow. Common forensic efficiency parameters were calculated.

RESULTS

The analysed X-STRs show high values of genetic diversity, comparable to other STRs of more common use. No significant associations between markers were found. A slight population structure was detected between the two shores of the Mediterranean. The X-STRs studied here present a similar degree of variability to that of other X-STRs used in forensics.

CONCLUSION

Tandem-repeated dinucleotides are a good tool for evidencing population differences here. Forensic parameters indicate that the dinucleotide X-STRs are suitable for forensic use.

Genomic diversity and population structure of three autochthonous Greek sheep breeds assessed with genome-wide DNA arrays

In the present study, genome-wide genotyping was applied to characterize the genetic diversity and population structure of three autochthonous Greek breeds: Boutsko, Karagouniko and Chios. Dairy sheep are among the most significant livestock species in Greece numbering approximately 9 million animals which are characterized by large phenotypic variation and reared under various farming systems. A total of 96 animals were genotyped with the Illumina's OvineSNP50K microarray beadchip, to study the population structure of the breeds and develop a specialized panel of single-nucleotide polymorphisms (SNPs), which could distinguish one breed from the others. Quality control on the dataset resulted in 46,125 SNPs, which were used to evaluate the genetic structure of the breeds. Population structure was assessed through principal component analysis (PCA) and admixture analysis, whereas inbreeding was estimated based on runs of homozygosity (ROHs) coefficients, genomic relationship matrix inbreeding coefficients (F_GRM) and patterns of linkage disequilibrium (LD). Associations between SNPs and breeds were analyzed with different inheritance models, to identify SNPs that distinguish among the breeds. Results showed high levels of genetic heterogeneity in the three breeds. Genetic distances among breeds were modest, despite their different ancestries. Chios and Karagouniko breeds were more genetically related to each other compared to Boutsko. Analysis revealed 3802 candidate SNPs that can be used to identify two-breed crosses and purebred animals. The present study provides, for the first time, data on the genetic background of three Greek indigenous dairy sheep breeds as well as a specialized marker panel that can be applied for traceability purposes as well as targeted genetic improvement schemes and conservation programs.

Landscape of genome-wide age-related DNA methylation in breast tissue

Despite known age-related DNA methylation (aDNAm) changes in breast tumors, little is known about aDNAm in normal breast tissues. Breast tissues from a cross-sectional study of 121 cancer-free women, were assayed for genome-wide DNA methylation. mRNA expression was assayed by microarray technology. Analysis of covariance was used to identify aDNAm’s. Altered methylation was correlated with expression of the corresponding gene and with DNA methyltransferase protein DNMT3A, assayed by immunohistochemistry. Publically-available TCGA-BRCA data were used for replication. 1,214 aDNAm’s were identified; 97% with increased methylation, and all on autosomes. Sites with increased methylation were predominantly in CpG lslands and non-enhancers. aDNAm’s with decreased methylation were generally located in intergenic regions, non-CpG Islands, and enhancers. Of the aDNAm’s identified, 650 are known to be involved in cancer, including ESR1 and beta-estradiol responsive genes. Expression of DNMT3A was positively associated with age. Two aDNAm’s showed borderline significant associations with DNMT3A expression; KRR1 (OR 6.57, 95% CI: 2.51–17.23) and DHRS12 (OR 6.08, 95% CI: 2.33–15.86). A subset of aDNAm’s co-localized within vulnerable regions for somatic mutations in cancers including breast cancer. Expression of C19orf48 was inversely and significantly correlated with its methylation level. In the TCGA dataset, 84% and 64% of the previously identified aDNAm’s were correlated with age in both normal-adjacent and tumor breast tissues, with differential associations by histological subtype. Given the similarity of findings in the breast tissues of healthy women and breast tumors, aDNAm’s may be one pathway for increased breast cancer risk with age.

Inherent X-Linked Genetic Variability and Cellular Mosaicism Unique to Females Contribute to Sex-Related Differences in the Innate Immune Response

Females have a longer lifespan and better general health than males. Considerable number of studies also demonstrated that, after trauma and sepsis, females present better outcomes as compared to males indicating sex-related differences in the innate immune response. The current notion is that differences in the immuno-modulatory effects of sex hormones are the underlying causative mechanism. However, the field remains controversial and the exclusive role of sex hormones has been challenged. Here, we propose that polymorphic X-linked immune competent genes, which are abundant in the population are important players in sex-based immuno-modulation and play a key role in causing sex-related outcome differences following trauma or sepsis. We describe the differences in X chromosome (ChrX) regulation between males and females and its consequences in the context of common X-linked polymorphisms at the individual as well as population level. We also discuss the potential pathophysiological and immune-modulatory aspects of ChrX cellular mosaicism, which is unique to females and how this may contribute to sex-biased immune-modulation. The potential confounding effects of ChrX skewing of cell progenitors at the bone marrow is also presented together with aspects of acute trauma-induced de novo ChrX skewing at the periphery. In support of the hypothesis, novel observations indicating ChrX skewing in a female trauma cohort as well as case studies depicting the temporal relationship between trauma-induced cellular skewing and the clinical course are also described. Finally, we list and discuss a selected set of polymorphic X-linked genes, which are frequent in the population and have key regulatory or metabolic functions in the innate immune response and, therefore, are primary candidates for mediating sex-biased immune responses. We conclude that sex-related differences in a variety of disease processes including the innate inflammatory response to injury and infection may be related to the abundance of X-linked polymorphic immune-competent genes, differences in ChrX regulation, and inheritance patterns between the sexes and the presence of X-linked cellular mosaicism, which is unique to females.

Genome Architecture and Evolution of a Unichromosomal Asexual Nematode

Asexual reproduction in animals, though rare, is the main or exclusive mode of reproduction in some long-lived lineages. The longevity of asexual clades may be correlated with the maintenance of heterozygosity by mechanisms that rearrange genomes and reduce recombination. Asexual species thus provide an opportunity to gain insight into the relationship between molecular changes, genome architecture, and cellular processes. Here we report the genome sequence of the parthenogenetic nematode Diploscapter pachys with only one chromosome pair. We show that this unichromosomal architecture is shared by a long-lived clade of asexual nematodes closely related to the genetic model organism Caenorhabditis elegans. Analysis of the genome assembly reveals that the unitary chromosome arose through fusion of six ancestral chromosomes, with extensive rearrangement among neighboring regions. Typical nematode telomeres and telomeric protection-encoding genes are lacking. Most regions show significant heterozygosity; homozygosity is largely concentrated to one region and attributed to gene conversion. Cell-biological and molecular evidence is consistent with the absence of key features of meiosis I, including synapsis and recombination. We propose that D. pachys preserves heterozygosity and produces diploid embryos without fertilization through a truncated meiosis. As a prelude to functional studies, we demonstrate that D. pachys is amenable to experimental manipulation by RNA interference.

Genetic characterization of Guinea-Bissau using a 12 X-chromosomal STR system: Inferences from a multiethnic population

A male West African sample from Guinea-Bissau (West-African coast) was genetically analyzed using 12 X chromosomal short tandem repeats that are grouped into four haplotype groups. Linkage disequilibrium was tested (p≤0.0008) and association was detected for the majority of markers in three out of the four studied haplotype clusters. The sample of 332 unrelated individuals analyzed in this study belonged to several recognized ethnic groups (n=18) which were used to evaluate the genetic variation of Guinea-Bissau's population. Pairwise genetic distances (F_ST) did not reveal significant differences among the majority of groups. An additional 110 samples from other countries also belonging to West Africa were as well compared with the sample of Guinea-Bissau. No significant differences were found between these two groups of West African individuals, supporting the genetic homogeneity of this region on the X chromosome level. The generation of over 100 DNA West African sequences provided new insights into the repeat sequence structure of some of the present X-STRs. Parameters for forensic evaluation were also calculated for each X-STR, supporting the potential application of these markers in typical kinship scenarios. Also, the high power of discrimination values for samples of female and male origin observed in this study, confirms the usefulness of the present X-STRs in identification analysis.

Deconstructing isolation-by-distance: The genomic consequences of limited dispersal

Geographically limited dispersal can shape genetic population structure and result in a correlation between genetic and geographic distance, commonly called isolation-by-distance. Despite the prevalence of isolation-by-distance in nature, to date few studies have empirically demonstrated the processes that generate this pattern, largely because few populations have direct measures of individual dispersal and pedigree information. Intensive, long-term demographic studies and exhaustive genomic surveys in the Florida Scrub-Jay (Aphelocoma coerulescens) provide an excellent opportunity to investigate the influence of dispersal on genetic structure. Here, we used a panel of genome-wide SNPs and extensive pedigree information to explore the role of limited dispersal in shaping patterns of isolation-by-distance in both sexes, and at an exceedingly fine spatial scale (within ~10 km). Isolation-by-distance patterns were stronger in male-male and male-female comparisons than in female-female comparisons, consistent with observed differences in dispersal propensity between the sexes. Using the pedigree, we demonstrated how various genealogical relationships contribute to fine-scale isolation-by-distance. Simulations using field-observed distributions of male and female natal dispersal distances showed good agreement with the distribution of geographic distances between breeding individuals of different pedigree relationship classes. Furthermore, we built coalescent simulations parameterized by the observed dispersal curve, population density, and immigration rate, and showed how incorporating these extensions to Malécot’s theory of isolation-by-distance allows us to accurately reconstruct observed sex-specific isolation-by-distance patterns in autosomal and Z-linked SNPs. Therefore, patterns of fine-scale isolation-by-distance in the Florida Scrub-Jay can be well understood as a result of limited dispersal over contemporary timescales.

Dispersal is a fundamental component of the life history of most organisms and therefore influences many biological processes. Dispersal is particularly important in creating genetic structure on the landscape. We often observe a pattern of decreased genetic relatedness between individuals as geographic distances increases, or isolation-by-distance. This pattern is particularly pronounced in organisms with extremely short dispersal distances. Despite the ubiquity of isolation-by-distance patterns in nature, there are few examples that explicitly demonstrate how limited dispersal influences spatial genetic structure. Here we investigate the processes that result in spatial genetic structure using the Florida Scrub-Jay, a bird with extremely limited dispersal behavior and extensive genome-wide data. We take advantage of the long-term monitoring of a contiguous population of Florida Scrub-Jays, which has resulted in a detailed pedigree and measurements of dispersal for hundreds of individuals. We show how limited dispersal results in close genealogical relatives living closer together geographically, which generates a strong pattern of isolation-by-distance at an extremely small spatial scale (<10 km) in just a few generations. Given the detailed dispersal, pedigree, and genomic data, we can achieve a fairly complete understanding of how dispersal shapes patterns of genetic diversity over short spatial scales.

DNA Commission of the International Society for Forensic Genetics (ISFG): Guidelines on the use of X-STRs in kinship analysis

Forensic genetic laboratories perform an increasing amount of genetic analyses of the X chromosome, in particular to solve complex cases of kinship analysis. For some biological relationships X-chromosomal markers can be more informative than autosomal markers, and there are a large number of markers, methods and databases that have been described for forensic use. Due to their particular mode of inheritance, and their physical location on a single chromosome, some specific considerations are required when estimating the weight of evidence for X-chromosomal marker DNA data. The DNA Commission of the International Society for Forensic Genetics (ISFG) hereby presents guidelines and recommendations for the use of X-chromosomal markers in kinship analysis with a special focus on the biostatistical evaluation. Linkage and linkage disequilibrium (association of alleles) are of special importance for such evaluations and these concepts and the implications for likelihood calculations are described in more detail. Furthermore it is important to use appropriate computer software that accounts for linkage and linkage disequilibrium among loci, as well as for mutations. Even though some software exist, there is still a need for further improvement of dedicated software.

Signatures of human European Palaeolithic expansion shown by resequencing of non-recombining X-chromosome segments

Human genetic diversity in Europe has been extensively studied using uniparentally inherited sequences (mitochondrial DNA (mtDNA) and the Y chromosome), which reveal very different patterns indicating sex-specific demographic histories. The X chromosome, haploid in males and inherited twice as often from mothers as from fathers, could provide insights into past female behaviours, but has not been extensively investigated. Here, we use HapMap single-nucleotide polymorphism data to identify genome-wide segments of the X chromosome in which recombination is historically absent and mutations are likely to be the only source of genetic variation, referring to these as phylogeographically informative haplotypes on autosomes and X chromosome (PHAXs). Three such sequences on the X chromosome spanning a total of ~49 kb were resequenced in 240 males from Europe, the Middle East and Africa at an average coverage of 181 ×. These PHAXs were confirmed to be essentially non-recombining across European samples. All three loci show highly homogeneous patterns across Europe and are highly differentiated from the African sample. Star-like structures of European-specific haplotypes in median-joining networks indicate past population expansions. Bayesian skyline plots and time-to-most-recent-common-ancestor estimates suggest expansions pre-dating the Neolithic transition, a finding that is more compatible with data on mtDNA than the Y chromosome, and with the female bias of X-chromosomal inheritance. This study demonstrates the potential of the use of X-chromosomal haplotype blocks, and the utility of the accurate ascertainment of rare variants for inferring human demographic history.

Differences in the effective population sizes of males and females do not require differences in their distribution of offspring number

Difference in male and female effective population sizes has, at times, been attributed to both sexes having unequal variance in their number of offspring. Such difference is paralleled by the relative effective sizes of autosomes, sex chromosomes, and mitochondrial DNA. I develop a simple framework to calculate the inbreeding effective population sizes for loci with different modes of inheritance. In this framework, I separate the effects due to mating strategy and those due to genetic transmission. I then show that, in addition to differences in the variance in offspring number, skew in the male/female effective sizes can also be caused by family composition. This approach can be used to illustrate the effect of induced behaviors on the relative male and female effective population sizes. In particular, I show the impact of the one-child policy formerly implemented in the People's Republic of China on the relative male and female effective population sizes. Furthermore, I argue that, under some strong constraints on family structure, the concepts of male and female effective population sizes are invalid.

Deep sequencing of 10,000 human genomes

We report on the sequencing of 10,545 human genomes at 30×-40× coverage with an emphasis on quality metrics and novel variant and sequence discovery. We find that 84% of an individual human genome can be sequenced confidently. This high-confidence region includes 91.5% of exon sequence and 95.2% of known pathogenic variant positions. We present the distribution of over 150 million single-nucleotide variants in the coding and noncoding genome. Each newly sequenced genome contributes an average of 8,579 novel variants. In addition, each genome carries on average 0.7 Mb of sequence that is not found in the main build of the hg38 reference genome. The density of this catalog of variation allowed us to construct high-resolution profiles that define genomic sites that are highly intolerant of genetic variation. These results indicate that the data generated by deep genome sequencing is of the quality necessary for clinical use.

A model for the clustered distribution of SNPs in the human genome

Motivated by a non-random but clustered distribution of SNPs, we introduce a phenomenological model to account for the clustering properties of SNPs in the human genome. The phenomenological model is based on a preferential mutation to the closer proximity of existing SNPs. With the Hapmap SNP data, we empirically demonstrate that the preferential model is better for illustrating the clustered distribution of SNPs than the random model. Moreover, the model is applicable not only to autosomes but also to the X chromosome, although the X chromosome has different characteristics from autosomes. The analysis of the estimated parameters in the model can explain the pronounced population structure and the low genetic diversity of the X chromosome. In addition, correlation between the parameters reveals the population-wise difference of the mutation probability. These results support the mutational non-independence hypothesis against random mutation.

Impact of the X Chromosome and sex on regulatory variation

The X Chromosome, with its unique mode of inheritance, contributes to differences between the sexes at a molecular level, including sex-specific gene expression and sex-specific impact of genetic variation. Improving our understanding of these differences offers to elucidate the molecular mechanisms underlying sex-specific traits and diseases. However, to date, most studies have either ignored the X Chromosome or had insufficient power to test for the sex-specific impact of genetic variation. By analyzing whole blood transcriptomes of 922 individuals, we have conducted the first large-scale, genome-wide analysis of the impact of both sex and genetic variation on patterns of gene expression, including comparison between the X Chromosome and autosomes. We identified a depletion of expression quantitative trait loci (eQTL) on the X Chromosome, especially among genes under high selective constraint. In contrast, we discovered an enrichment of sex-specific regulatory variants on the X Chromosome. To resolve the molecular mechanisms underlying such effects, we generated chromatin accessibility data through ATAC-sequencing to connect sex-specific chromatin accessibility to sex-specific patterns of expression and regulatory variation. As sex-specific regulatory variants discovered in our study can inform sex differences in heritable disease prevalence, we integrated our data with genome-wide association study data for multiple immune traits identifying several traits with significant sex biases in genetic susceptibilities. Together, our study provides genome-wide insight into how genetic variation, the X Chromosome, and sex shape human gene regulation and disease.

An expanded sequence context model broadly explains variability in polymorphism levels across the human genome

The rate of single nucleotide polymorphism varies substantially across the human genome and fundamentally influences evolution and incidence of genetic disease. Previous studies have only considered the immediate flanking nucleotides around a polymorphic site –the site’s trinucleotide sequence context– to study polymorph levels across the genome. Moreover, the impact of larger sequence contexts has not been fully clarified, even though context substantially influences rates of polymorphism. Using a new statistical framework and data from the 1000 Genomes Project, we demonstrate that a heptanucleotide context explains >81% of variability in substitution probabilities, revealing new mutation-promoting motifs at ApT dinucleotide, CAAT, and TACG sequences. Our approach also identifies previously undocumented variability in C-to-T substitutions at CpG sites, which is not immediately explained by differential methylation intensity. Using our model, we present informative substitution intolerance scores for genes and a new intolerance score for amino acids, and we demonstrate clinical use of the model in neuropsychiatric diseases.

Complexity and diversity of F8 genetic variations in the 1000 genomes

BACKGROUND

Hemophilia A (HA) is an X-linked bleeding disorder caused by deleterious mutations in the coagulation factor VIII gene (F8). To date, F8 mutations have been documented predominantly in European subjects and in American subjects of European descent. Information on F8 variants in individuals of more diverse ethnic backgrounds is limited.

OBJECTIVES

To discover novel and rare F8 variants, and to characterize F8 variants in diverse population backgrounds.

PATIENTS/METHODS

We analyzed 2535 subjects, including 26 different ethnicities, whose data were available from the 1000 Genomes Project (1000G) phase 3 dataset, for F8 variants and their potential functional impact.

RESULTS

We identified 3030 single nucleotide variants, 31 short deletions and insertions (Indels) and a large, 497 kb, deletion. Among all variants, 86.4% were rare variants and 55.6% were novel. Eighteen variants previously associated with HA were found in our study. Most of these 'HA variants' were ethnic-specific with low allele frequency; however, one variant (p.M2257V) was present in 27% of African subjects. The p.E132D, p.T281A, p.A303V and p.D422H 'HA variants' were identified only in males. Twelve novel missense variants were predicted to be deleterious. The large deletion was discovered in eight female subjects without affecting F8 transcription and the transcription of genes on the X chromosome.

CONCLUSION

Characterizing F8 in the 1000G project highlighted the complexity of F8 variants and the importance of interrogating genetic variants on multiple ethnic backgrounds for associations with bleeding and thrombosis. The haplotype analysis and the orientation of duplicons that flank the large deletion suggested that the deletion was recurrent and originated by homologous recombination.

Genetic sharing and heritability of paediatric age of onset autoimmune diseases

Autoimmune diseases (AIDs) are polygenic diseases affecting 7-10% of the population in the Western Hemisphere with few effective therapies. Here, we quantify the heritability of paediatric AIDs (pAIDs), including JIA, SLE, CEL, T1D, UC, CD, PS, SPA and CVID, attributable to common genomic variations (SNP-h(2)). SNP-h(2) estimates are most significant for T1D (0.863±s.e. 0.07) and JIA (0.727±s.e. 0.037), more modest for UC (0.386±s.e. 0.04) and CD (0.454±0.025), largely consistent with population estimates and are generally greater than that previously reported by adult GWAS. On pairwise analysis, we observed that the diseases UC-CD (0.69±s.e. 0.07) and JIA-CVID (0.343±s.e. 0.13) are the most strongly correlated. Variations across the MHC strongly contribute to SNP-h(2) in T1D and JIA, but does not significantly contribute to the pairwise rG. Together, our results partition contributions of shared versus disease-specific genomic variations to pAID heritability, identifying pAIDs with unexpected risk sharing, while recapitulating known associations between autoimmune diseases previously reported in adult cohorts.

Chromosomal Rearrangements as Barriers to Genetic Homogenization between Archaic and Modern Humans

Chromosomal rearrangements, which shuffle DNA throughout the genome, are an important source of divergence across taxa. Using a paired-end read approach with Illumina sequence data for archaic humans, I identify changes in genome structure that occurred recently in human evolution. Hundreds of rearrangements indicate genomic trafficking between the sex chromosomes and autosomes, raising the possibility of sex-specific changes. Additionally, genes adjacent to genome structure changes in Neanderthals are associated with testis-specific expression, consistent with evolutionary theory that new genes commonly form with expression in the testes. I identify one case of new-gene creation through transposition from the Y chromosome to chromosome 10 that combines the 5'-end of the testis-specific gene Fank1 with previously untranscribed sequence. This new transcript experienced copy number expansion in archaic genomes, indicating rapid genomic change. Among rearrangements identified in Neanderthals, 13% are transposition of selfish genetic elements, whereas 32% appear to be ectopic exchange between repeats. In Denisovan, the pattern is similar but numbers are significantly higher with 18% of rearrangements reflecting transposition and 40% ectopic exchange between distantly related repeats. There is an excess of divergent rearrangements relative to polymorphism in Denisovan, which might result from nonuniform rates of mutation, possibly reflecting a burst of transposable element activity in the lineage that led to Denisovan. Finally, loci containing genome structure changes show diminished rates of introgression from Neanderthals into modern humans, consistent with the hypothesis that rearrangements serve as barriers to gene flow during hybridization. Together, these results suggest that this previously unidentified source of genomic variation has important biological consequences in human evolution.

Association of 370-371insACA, 494T>C, and 1423C>T haplotype in ubiquitin-specific protease 26 gene and male infertility: a meta-analysis

Whether the 370-371insACA, 494T>C, and 1423C>T haplotype in ubiquitin-specific protease 26 (USP26) gene is associated with male infertility is controversial. To clarify this issue, we conducted a meta-analysis based on the most recent studies. Eligible studies were screened by using PubMed and Embase. Pooled odd ratio (OR) with 95% confidence interval (CI) was calculated with fixed effect models. Ten studies with 1603 patients and 2505 controls were included. Overall, the results indicated that there was an association between the haplotype and male infertile risk (OR = 1.74, 95% CI: 1.09–2.77). The OR calculated based on the five studies in Asia and three in Europe was 1.96 (95% CI: 1.05–3.67) and 1.54 (95% CI: 0.75–3.16) respectively, however, the OR was 0.86 (95% CI: 0.05–15.29) based on the two investigations in America. In addition, the data from the patients with azoospermia (AZO) showed an increased pooled OR of 2.35 (95% CI: 1.22–4.50). In contrast, the studies with oligoasthenoteratozoospermia (OAT) exhibited that the pooled OR was 0.97 (95% CI: 0.43–2.16). Our analyses indicate that there is an association of alteration in USP26 with male infertility, especially in AZO and Asian population.

A Genome-Wide Survey of Date Palm Cultivars Supports Two Major Subpopulations in Phoenix dactylifera

The date palm (Phoenix dactylifera L.) is one of the oldest cultivated trees and is intimately tied to the history of human civilization. There are hundreds of commercial cultivars with distinct fruit shapes, colors, and sizes growing mainly in arid lands from the west of North Africa to India. The origin of date palm domestication is still uncertain, and few studies have attempted to document genetic diversity across multiple regions. We conducted genotyping-by-sequencing on 70 female cultivar samples from across the date palm–growing regions, including four Phoenix species as the outgroup. Here, for the first time, we generate genome-wide genotyping data for 13,000–65,000 SNPs in a diverse set of date palm fruit and leaf samples. Our analysis provides the first genome-wide evidence confirming recent findings that the date palm cultivars segregate into two main regions of shared genetic background from North Africa and the Arabian Gulf. We identify genomic regions with high densities of geographically segregating SNPs and also observe higher levels of allele fixation on the recently described X-chromosome than on the autosomes. Our results fit a model with two centers of earliest cultivation including date palms autochthonous to North Africa. These results adjust our understanding of human agriculture history and will provide the foundation for more directed functional studies and a better understanding of genetic diversity in date palm.

Genetic ancestry of a Moroccan population as inferred from autosomal STRs

Detecting population substructure and ancestry is a critical issue for both association studies of health behaviors and forensic genetics. Determining aspects of a population's genetic history as potential sources of substructure can aid in design of future genetic studies. Within this context, fifteen autosomal short tandem repeat (STR), were used to examine population genetic structure and hypotheses of the origin of the modern Moroccan population from individuals belonging to three different ethnical groups from Morocco (Arab, Berber and Sahrawi), by comparing their autosomal STR variation with that of neighboring and non-neighboring populations in North Africa, Europe and Middle East as well as proposed ancestral populations in Morocco (Berber). We report on the results that the gradient of North African ancestry accounts for previous observations of low levels of sharing with Near East and a substantially increased gene flow especially from Morocco and Spain.

Linkage disequilibrium patterns and persistence of phase in purebred and crossbred pig (Sus scrofa) populations

Background

Genomic selection and genomic wide association studies are widely used methods that aim to exploit the linkage disequilibrium (LD) between markers and quantitative trait loci (QTL). Securing a sufficiently large set of genotypes and phenotypes can be a limiting factor that may be overcome by combining data from multiple breeds or using crossbred information. However, the estimated effect of a marker in one breed or a crossbred can only be useful for the selection of animals in another breed if there is a correspondence of the phase between the marker and the QTL across breeds. Using data of five pure pig (Sus scrofa) lines (SL1, SL2, SL3, DL1, DL2), one F₁ cross (DLF1) and two commercial finishing crosses (TER1 and TER2), the objectives of this study were: (i) to compare the equality of LD decay curves of different pig populations; and (ii) to evaluate the persistence of the LD phase across lines or final crosses.

Results

Almost all of the lines presented different extents of LD, except for the SL2 and DL3, both of which exhibited the same extent of LD. Similar levels of LD over large distances were found in crossbred and pure lines. The crossbred animals (DLF1, TER1 and TER2) presented a high persistence of phase with their parental lines, suggesting that the available porcine single nucleotide polymorphism (SNP) chip should be dense enough to include markers that have the same LD phase with QTL across crossbred and parental pure lines. The persistence of phase across pure lines varied considerably between the different line comparisons; however, correlations were above 0.8 for all line comparisons when marker distances were smaller than 50 kb.

Conclusions

This study showed that crossbred populations could be very useful as a reference for the selection of pure lines by means of the available SNP chip panel. Here, we also pinpoint pure lines that could be combined in a multiline training population. However, if multiline reference populations are used for genomic selection, the required density of SNP panels should be higher compared with a single breed reference population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-014-0126-3) contains supplementary material, which is available to authorized users.

Multilocus nuclear DNA markers reveal population structure and demography of Anopheles minimus

Utilization of multiple putatively neutral DNA markers for inferring evolutionary history of species population is considered to be the most robust approach. Molecular population genetic studies have been conducted in many species of Anopheles genus, but studies based on single nucleotide polymorphism (SNP) data are still very scarce. Anopheles minimus is one of the principal malaria vectors of Southeast (SE) Asia including the Northeastern (NE) India. Although population genetic studies with mitochondrial genetic variation data have been utilized to infer phylogeography of the SE Asian populations of this species, limited information on the population structure and demography of Indian An. minimus is available. We herewith have developed multilocus nuclear genetic approach with SNP markers located in X chromosome of An. minimus in eight Indian and two SE Asian population samples (121 individual mosquitoes in total) to infer population history and test several hypotheses on the phylogeography of this species. While the Thai population sample of An. minimus presented the highest nucleotide diversity, majority of the Indian samples were also fairly diverse. In general, An. minimus populations were moderately substructured in the distribution range covering SE Asia and NE India, largely falling under three distinct genetic clusters. Moreover, demographic expansion events could be detected in the majority of the presently studied populations of An. minimus. Additional DNA sequencing of the mitochondrial COII region in a subset of the samples (40 individual mosquitoes) corroborated the existing hypothesis of Indian An. minimus falling under the earlier reported mitochondrial lineage B.

Molecular population genetics of the NADPH cytochrome P450 reductase (CPR) gene in Anopheles minimus

Development of insecticide resistance (IR) in mosquito vectors is a primary huddle to malaria control program. Since IR has genetic basis, and genes constantly evolve with response to environment for adaptation to organisms, it is important to know evolutionary pattern of genes conferring IR in malaria vectors. The mosquito Anopheles minimus is a major malaria vector of the Southeast (SE) Asia and India and is susceptible to all insecticides, and thus of interest to know if natural selection has shaped variations in the gene conferring IR. If not, the DNA fragment of such a gene could be used to infer population structure and demography of this species of malaria vector. We have therefore sequenced a ~569 bp DNA segment of the NADPH cytochrome P450 reductase (CPR) gene (widely known to confer IR) in 123 individuals of An. minimus collected in 10 different locations (eight Indian, one Thai and one Vietnamese). Two Indian population samples were completely mono-morphic in the CPR gene. In general, low genetic diversity was found with no evidence of natural selection in this gene. The data were therefore analyzed to infer population structure and demography of this species. The 10 populations could be genetically differentiated into four different groups; the samples from Thailand and Vietnam contained high nucleotide diversity. All the 10 populations conform to demographic equilibrium model with signature of past population expansion in four populations. The results in general indicate that the An. minimus mosquitoes sampled in the two SE Asian localities contain several genetic characteristics of being parts of the ancestral population.

The peopling of Greenland: further insights from the analysis of genetic diversity using autosomal and X-chromosomal markers

The peopling of Greenland has a complex history shaped by population migrations, isolation and genetic drift. The Greenlanders present a genetic heritage with components of European and Inuit groups; previous studies using uniparentally inherited markers in Greenlanders have reported evidence of a sex-biased, admixed genetic background. This work further explores the genetics of the Greenlanders by analysing autosomal and X-chromosomal data to obtain deeper insights into the factors that shaped the genetic diversity in Greenlanders. Fourteen Greenlandic subsamples from multiple geographical settlements were compared to assess the level of genetic substructure in the Greenlandic population. The results showed low levels of genetic diversity in all sets of the genetic markers studied, together with an increased number of X-chromosomal loci in linkage disequilibrium in relation to the Danish population. In the broader context of worldwide populations, Greenlanders are remarkably different from most populations, but they are genetically closer to some Inuit groups from Alaska. Admixture analyses identified an Inuit component in the Greenlandic population of approximately 80%. The sub-populations of Ammassalik and Nanortalik are the least diverse, presenting the lowest levels of European admixture. Isolation-by-distance analyses showed that only 16% of the genetic substructure of Greenlanders is most likely to be explained by geographic barriers. We suggest that genetic drift and a differentiated settlement history around the island explain most of the genetic substructure of the population in Greenland.

X-chromosome STR markers data in a Cabo Verde immigrant population of Lisboa

Population genetic data of 12 X chromosomal short tandem repeats markers (DXS10074, DXS10079, DXS10101, DXS10103, DXS10134, DXS10135, DXS10146, DXS10148, DXS7132, DXS7423, DXS8378 and HPRTB) were analysed in 54 females and 95 males of an immigrant population from Cabo Verde living in Lisboa. The obtained results for forensic statistical parameters such as observed heterozigosity, polymorphism information content, power of discrimination and mean exclusion chance, based on single allele frequencies, reveal that this multiplex system is highly informative and can represent an important tool for genetic identification purposes in the immigrant population of Cabo Verde. Since the studied short tandem repeats genetic markers are distributed on four linkage groups, that can provide independent genotype information, we studied those groups as haploytes. The forensic efficiency parameters for the linked groups were all higher than 0.97, with linkage group I being the most polymorphic and linkage group III the less informative.

Genomic divergence of zebu and taurine cattle identified through high-density SNP genotyping

BACKGROUND

Natural selection has molded evolution across all taxa. At an arguable date of around 330,000 years ago there were already at least two different types of cattle that became ancestors of nearly all modern cattle, the Bos taurus taurus more adapted to temperate climates and the tropically adapted Bos taurus indicus. After domestication, human selection exponentially intensified these differences. To better understand the genetic differences between these subspecies and detect genomic regions potentially under divergent selection, animals from the International Bovine HapMap Experiment were genotyped for over 770,000 SNP across the genome and compared using smoothed F(ST). The taurine sample was represented by ten breeds and the contrasting zebu cohort by three breeds.

RESULTS

Each cattle group evidenced similar numbers of polymorphic markers well distributed across the genome. Principal components analyses and unsupervised clustering confirmed the well-characterized main division of domestic cattle. The top 1% smoothed F(ST), potentially associated to positive selection, contained 48 genomic regions across 17 chromosomes. Nearly half of the top F(ST) signals (n = 22) were previously detected using a lower density SNP assay. Amongst the strongest signals were the BTA7:~50 Mb and BTA14:~25 Mb; both regions harboring candidate genes and different patterns of linkage disequilibrium that potentially represent intrinsic differences between cattle types. The bottom 1% of the smoothed F(ST) values, potentially associated to balancing selection, included 24 regions across 13 chromosomes. These regions often overlap with copy number variants, including the highly variable region at BTA23:~24 Mb that harbors a large number of MHC genes. Under these regions, 318 unique Ensembl genes are annotated with a significant overrepresentation of immune related pathways.

CONCLUSIONS

Genomic regions that are potentially linked to purifying or balancing selection processes in domestic cattle were identified. These regions are of particular interest to understand the natural and human selective pressures to which these subspecies were exposed to and how the genetic background of these populations evolved in response to environmental challenges and human manipulation.

Exome sequencing identifies potential risk variants for Mendelian disorders at high prevalence in Qatar

Exome sequencing of families of related individuals has been highly successful in identifying genetic polymorphisms responsible for Mendelian disorders. Here, we demonstrate the value of the reverse approach, where we use exome sequencing of a sample of unrelated individuals to analyze allele frequencies of known causal mutations for Mendelian diseases. We sequenced the exomes of 100 individuals representing the three major genetic subgroups of the Qatari population (Q1 Bedouin, Q2 Persian-South Asian, Q3 African) and identified 37 variants in 33 genes with effects on 36 clinically significant Mendelian diseases. These include variants not present in 1000 Genomes and variants at high frequency when compared with 1000 Genomes populations. Several of these Mendelian variants were only segregating in one Qatari subpopulation, where the observed subpopulation specificity trends were confirmed in an independent population of 386 Qataris. Premarital genetic screening in Qatar tests for only four out of the 37, such that this study provides a set of Mendelian disease variants with potential impact on the epidemiological profile of the population that could be incorporated into the testing program if further experimental and clinical characterization confirms high penetrance.

Cellular mosaicism for X-linked polymorphisms and IRAK1 expression presents a distinct phenotype and improves survival following sepsis

ChrX cellular mosaicism for X-linked genetic polymorphisms in females versus the single ChrX representation in males denotes a genetic difference, which may contribute to gender bias in the inflammatory response. This hypothesis was tested in female F1 offspring of consomic mice (BL6J-ChrX(A/J)/NaJ) that were homokaryotic or mosaic for the active BL6 and AJ ChrXs or for IRAK1 deficiency linked to the BL6 ChrX. Sepsis was initiated by CLP. IRAK1-deficient and IRAK1-mosaic mice showed similar protection from sepsis-induced mortality and reduced IL-6 and IL-10 release compared with WT. BM cellularity and blood B cell counts were increased in naive IRAK1-mosaic mice compared with WT-mosaic or IRAK1-deficient animals. Sepsis-induced BM cell depletion was greater in IRAK1-mosaic mice compared with WT-mosaic or IRAK1-deficient subjects, whereas splenic B and T cell depletion was less in IRAK1-mosaic and IRAK1-deficient than WT-mosaic mice. Skewing toward AJ or BL6-ChrX-expressing cells was assessed by testing allele-specific expression of strain-variant Xkrx and BTK genes. In naive IRAK1-mosaic mice, BM and blood cells with the active BL6-ChrX, were greater than cells expressing the AJ-ChrX (cell ratio 2.5 in IRAK1-mosaic; 1.5 in WT-mosaic mice). Sepsis decreased cell ratios more in IRAK1-mosaic than in WT-mosaic mice. The study reveals functional variability in cellular mosaicism for IRAK1 expression and natural X-linked polymorphisms during sepsis. Mosaicism for IRAK1 expression is accompanied by skewing toward deficient immune cell populations, producing a phenotype that is preconditioned for improved sepsis outcome similar to that observed in IRAK1 deficiency.

Sex-linked and autosomal microsatellites provide new insights into island populations of the tammar wallaby

The emerging availability of microsatellite markers from mammalian sex chromosomes provides opportunities to investigate both male- and female-mediated gene flow in wild populations, identifying patterns not apparent from the analysis of autosomal markers alone. Tammar wallabies (Macropus eugenii), once spread over the southern mainland, have been isolated on several islands off the Western Australian and South Australian coastlines for between 10,000 and 13,000 years. Here, we combine analyses of autosomal, Y-linked and X-linked microsatellite loci to investigate genetic variation in populations of this species on two islands (Kangaroo Island, South Australia and Garden Island, Western Australia). All measures of diversity were higher for the larger Kangaroo Island population, in which genetic variation was lowest at Y-linked markers and highest at autosomal markers (θ=3.291, 1.208 and 0.627 for autosomal, X-linked and Y-linked data, respectively). Greater relatedness among females than males provides evidence for male-biased dispersal in this population, while sex-linked markers identified genetic lineages not apparent from autosomal data alone. Overall genetic diversity in the Garden Island population was low, especially on the Y chromosome where most males shared a common haplotype, and we observed high levels of inbreeding and relatedness among individuals. Our findings highlight the utility of this approach for management actions, such as the selection of animals for translocation or captive breeding, and the ecological insights that may be gained by combining analyses of microsatellite markers on sex chromosomes with those derived from autosomes.

Amish revisited: next-generation sequencing studies of psychiatric disorders among the Plain people

The rapid development of next-generation sequencing (NGS) technology has led to renewed interest in the potential contribution of rarer forms of genetic variation to complex non-mendelian phenotypes such as psychiatric illnesses. Although challenging, family-based studies offer some advantages, especially in communities with large families and a limited number of founders. Here we revisit family-based studies of mental illnesses in traditional Amish and Mennonite communities--known collectively as the Plain people. We discuss the new opportunities for NGS in these populations, with particular emphasis on investigating psychiatric disorders. We also address some of the challenges facing NGS-based studies of complex phenotypes in founder populations.

Genomic conflict drives patterns of X-linked population structure in Drosophila neotestacea

Intragenomic conflict has the potential to cause widespread changes in patterns of genetic diversity and genome evolution. In this study, we investigate the consequences of sex-ratio (SR) drive on the population genetic patterns of the X-chromosome in Drosophila neotestacea. An SR X-chromosome prevents the maturation of Y-bearing sperm during male spermatogenesis and thus is transmitted to ~100% of the offspring, nearly all of which are daughters. Selection on the rest of the genome to suppress SR can be strong, and the resulting conflict over the offspring sex ratio can result in the accumulation of multiple loci on the X-chromosome that are necessary for the expression of drive. We surveyed variation at 12 random X-linked microsatellites across 16 populations of D. neotestacea that range in SR frequency from 0% to 30%. First, every locus was differentiated between SR and wild-type chromosomes, and this drives genetic structure at the X-chromosome. Once the association with SR is accounted for, the patterns of differentiation among populations are similar to the autosomes. Second, within wild-type chromosomes, the relative heterozygosity is reduced in populations with an increased prevalence of drive, and the heterozygosity of SR chromosomes is higher than expected based on its prevalence. The combination of the relatively high prevalence of SR drive and the structuring of polymorphism between the SR and wild-type chromosomes suggests that genetic conflict because of SR drive has had significant consequences on the patterns of X-linked polymorphism and thus also probably affects the tempo of X-chromosome evolution in D. neotestacea.

Uncovering networks from genome-wide association studies via circular genomic permutation

Genome-wide association studies (GWAS) aim to detect single nucleotide polymorphisms (SNP) associated with trait variation. However, due to the large number of tests, standard analysis techniques impose highly stringent significance thresholds, leaving potentially associated SNPs undetected, and much of the trait genetic variation unexplained. Pathway- and network-based methodologies applied to GWAS aim to detect associations missed by standard single-marker approaches. The complex and non-random architecture of the genome makes it a challenge to derive an appropriate testing framework for such methodologies. We developed a rapid and simple permutation approach that uses GWAS SNP association results to establish the significance of pathway associations while accounting for the linkage disequilibrium structure of SNPs and the clustering of functionally related elements in the genome. All SNPs used in the GWAS are placed in a “circular genome” according to their location. Then the complete set of SNP association P values are permuted by rotation with respect to the genomic locations of the SNPs. Once these “simulated” P values are assigned, the joint gene P values are calculated using Fisher’s combination test, and the association of pathways is tested using the hypergeometric test. The circular genomic permutation approach was applied to a human genome-wide association dataset. The data consists of 719 individuals from the ORCADES study genotyped for ∼300,000 SNPs and measured for 51 traits ranging from physical to biochemical measurements. KEGG pathways (n = 225) were used as the sets of pathways to be tested. Our results demonstrate that the circular genomic permutations provide robust association P values. The non-permuted hypergeometric analysis generates ∼1400 pathway-trait combination results with an association P value more significant than P ≤ 0.05, whereas applying circular genomic permutation reduces the number of significant results to a more credible 40% of that value. The circular permutation software (“genomicper”) is available as an R package at http://cran.r-project.org/.

Development of multilocus putatively neutral DNA markers in the X-chromosome for population genetic studies in humans

BACKGROUND

It has now been well documented that the type (coding, non-coding) and location (nuclear, mitochondrial etc.) of genetic markers heavily influence evolutionary inferences; realistic assumptions can be drawn if multiple putatively neutral DNA fragments spread across the genome are used.

AIM

To infer human population history, Single Nucleotide Polymorphisms (SNPs), located in the non-coding regions of different genes in the X-chromosome have been developed as 'putatively neutral markers'.

SUBJECTS AND METHODS

A population sample consisting of 16 male individuals from the western part of India was utilized for sequencing eight DNA fragments located in introns of three genes (Duchenne muscular dystrophy, Factor IX and Pyruvate dehydrogenase E1 sub-unit) on the human X-chromosome. PCR amplification and DNA sequencing confirmed the polymorphic status of all the fragments.

RESULTS

Twenty nine SNPs were found to be segregating in the Western Indian population samples. Using these SNPs the nucleotide diversity and demographic parameters of the Western Indian population were estimated. Several tests of neutrality ascertained that all eight fragments evolve putatively neutrally. Further, linkage disequilibrium analyses confirmed this fact.

CONCLUSION

All eight DNA fragments seem to bear the characteristics to be considered as 'putatively neutral genetic markers' and thus, could be utilized for inference of human population and demographic histories.

What Can Domesticated Genes Tell Us about the Intron Gain in Mammals?

Domesticated genes, originating from retroelements or from DNA-transposons, constitute an ideal system for testing the hypothesis on the absence of intron gain in mammals. Since single-copy domesticated genes originated from the intronless multicopy transposable elements, the ancestral intron state for domesticated genes is zero. A phylogenomic approach has been used to analyse all domesticated genes in mammals and chordates that originated from the coding parts of transposable elements. A significant amount of intron gain was found only in domesticated genes of placental mammals, where more than 70 cases were identified. De novo gained introns show clear positional bias, since they are distributed mainly in 5′ UTR and coding regions, while 3′ UTR introns are very rare. In the coding regions of some domesticated genes up to 8 de novo gained introns have been found. Surprisingly, the majority of intron gains have occurred in the ancestor of placental mammals. Domesticated genes could constitute an excellent system on which to analyse the mechanisms of intron gain. This paper summarizes the current understanding of intron gain in mammals.

Worldwide genetic relationships of pigs as inferred from X chromosome SNPs

The phylogeography of the porcine X chromosome has not been studied despite the unique characteristics of this chromosome. Here, we genotyped 59 single nucleotide polymorphisms (SNPs) in 312 pigs from around the world, representing 39 domestic breeds and wild boars in 30 countries. Overall, widespread commercial breeds showed the highest heterozygosity values, followed by African and American populations. Structuring, as inferred from FST and analysis of molecular variance, was consistently larger in the non-pseudoautosomal (NPAR) than in the pseudoautosomal regions (PAR). Our results show that genetic relationships between populations can vary widely between the NPAR and the PAR, underscoring the fact that their genetic trajectories can be quite different. NPAR showed an increased commercial-like genetic component relative to the PAR, probably because human selection processes to obtain individuals with high productive parameters were mediated by introgressing boars rather than sows.

Refining the genetic portrait of Portuguese Roma through X-chromosomal markers

Due to differences in transmission between X-chromosomal and autosomal DNA, the comparison of data derived from both markers allows deeper insight into the forces that shape the patterns of genetic diversity in populations. In this study, we applied this comparative approach to a sample of Portuguese Roma (Gypsies) by analyzing 43 X-chromosomal markers and 53 autosomal markers. Portuguese individuals of non-Gypsy ancestry were also studied. Compared with the host population, reduced levels of diversity on the X chromosome and autosomes were detected in Gypsies; this result was in line with known patterns of genetic diversity typical of Roma groups. As a consequence of the complex demographic past of the Roma, during which admixture and genetic drift played major roles, the amount of linkage disequilibrium (LD) on the X chromosome in Gypsies was considerably higher than that observed in non-Gypsies. When the pattern of differentiation on the X chromosome was compared with that of autosomes, there was evidence for asymmetries in female and male effective population sizes during the admixture between Roma and non-Roma. This result supplements previous data provided by mtDNA and the Y chromosome, underlining the importance of using combined information from the X chromosome and autosomes to dissect patterns of genetic diversity. Following the out-of-India dispersion, the Roma acquired a complex genetic pattern that was influenced by drift and introgression with surrounding populations, with important contributions from both males and females. We provide evidence that a sex-biased admixture with Europeans is probably associated with the founding of the Portuguese Gypsies.

Model-based comparisons of phylogeographic scenarios resolve the intraspecific divergence of cactophilic Drosophila mojavensis

The cactophilic fly Drosophila mojavensis exhibits considerable intraspecific genetic structure across allopatric geographic regions and shows associations with different host cactus species across its range. The divergence between these populations has been studied for more than 60years, yet their exact historical relationships have not been resolved. We analysed sequence data from 15 intronic X-linked loci across populations from Baja California, mainland Sonora-Arizona and Mojave Desert regions under an isolation-with-migration model to assess multiple scenarios of divergence. We also compared the results with a pre-existing sequence data set of eight autosomal loci. We derived a population tree with Baja California placed at its base and link their isolation to Pleistocene climatic oscillations. Our estimates suggest the Baja California population diverged from an ancestral Mojave Desert/mainland Sonora-Arizona group around 230,000-270,000years ago, while the split between the Mojave Desert and mainland Sonora-Arizona populations occurred one glacial cycle later, 117,000-135,000years ago. Although we found these three populations to be effectively allopatric, model ranking could not rule out the possibility of a low level of gene flow between two of them. Finally, the Mojave Desert population showed a small effective population size, consistent with a historical population bottleneck. We show that model-based inference from multiple loci can provide accurate information on the historical relationships of closely related groups allowing us to set into historical context a classic system of incipient ecological speciation.