Václav Janoušek

Genome-wide architecture of reproductive isolation in a naturally occurring hybrid zone between Mus musculus musculus and M. m. domesticus

Studies of a hybrid zone between two house mouse subspecies (Mus musculus musculus and M. m. domesticus) along with studies using laboratory crosses reveal a large role for the X chromosome and multiple autosomal regions in reproductive isolation as a consequence of disrupted epistasis in hybrids. One limitation of previous work has been that most of the identified genomic regions have been large. The goal here is to detect and characterize precise genomic regions underlying reproductive isolation. We surveyed 1401 markers evenly spaced across the genome in 679 mice collected from two different transects. Comparisons between transects provide a means for identifying common patterns that likely reflect intrinsic incompatibilities. We used a genomic cline approach to identify patterns that correspond to epistasis. From both transects, we identified contiguous regions on the X chromosome in which markers were inferred to be involved in epistatic interactions. We then searched for autosomal regions showing the same patterns and found they constitute about 5% of autosomal markers. We discovered substantial overlap between these candidate regions underlying reproductive isolation and QTL for hybrid sterility identified in laboratory crosses. Analysis of gene content in these regions suggests a key role for several mechanisms, including the regulation of transcription, sexual conflict and sexual selection operating at both the postmating prezygotic and postzygotic stages of reproductive isolation. Taken together, these results indicate that speciation in two recently diverged (c. 0.5 Ma) house mouse subspecies is complex, involving many genes dispersed throughout the genome and associated with distinct functions.

Measures of linkage disequilibrium among neighbouring SNPs indicate asymmetries across the house mouse hybrid zone

Theory predicts that naturally occurring hybrid zones between genetically distinct taxa can move over space and time as a result of selection and/or demographic processes, with certain types of hybrid zones being more or less likely to move. Determining whether a hybrid zone is stationary or moving has important implications for understanding evolutionary processes affecting interactions in hybrid populations. However, direct observations of hybrid zone movement are difficult to make unless the zone is moving rapidly. Here, evidence for movement in the house mouse Mus musculus domesticus × Mus musculus musculus hybrid zone is provided using measures of LD and haplotype structure among neighbouring SNP markers from across the genome. Local populations of mice across two transects in Germany and the Czech Republic were sampled, and a total of 1301 mice were genotyped at 1401 markers from the nuclear genome. Empirical measures of LD provide evidence for extinction and (re)colonization in single populations and, together with simulations, suggest hybrid zone movement because of either geography‐dependent asymmetrical dispersal or selection favouring one subspecies over the other.

Functional Organization of the Genome May Shape the Species Boundary in the House Mouse

Genomic features such as rate of recombination and differentiation have been suggested to play a role in species divergence. However, the relationship of these phenomena to functional organization of the genome in the context of reproductive isolation remains unexplored. Here, we examine genomic characteristics of the species boundaries between two house mouse subspecies (Mus musculus musculus/M. m. domesticus). These taxa form a narrow semipermeable zone of secondary contact across Central Europe. Due to the incomplete nature of reproductive isolation, gene flow in the zone varies across the genome. We present an analysis of genomic differentiation, rate of recombination, and functional composition of genes relative to varying amounts of introgression. We assessed introgression using 1,316 autosomal single nucleotide polymorphism markers, previously genotyped in hybrid populations from three transects. We found a significant relationship between amounts of introgression and both genomic differentiation and rate of recombination with genomic regions of reduced introgression associated with higher genomic differentiation and lower rates of recombination, and the opposite for genomic regions of extensive introgression. We also found a striking functional polarization of genes based on where they are expressed in the cell. Regions of elevated introgression exhibit a disproportionate number of genes involved in signal transduction functioning at the cell periphery, among which olfactory receptor genes were found to be the most prominent group. Conversely, genes expressed intracellularly and involved in DNA binding were the most prevalent in regions of reduced introgression. We hypothesize that functional organization of the genome is an important driver of species divergence.

Experimental hepatic porphyria induced by oxadiazon in male mice and rats

Abstract Male mice and rats were fed a diet containing 1000 mg/kg of oxadiazon for 6 days. Fecal porphyrin was elevated in both mice and rats (from 30.2 to 1673.1 nmol/g dry wt and from 98.3 to 547.6 nmol/g dry wt, respectively). Porphyrin content of mouse bile was elevated from 5.7 to 439.0 nmol/g wet wt. Livers of treated animals contained mainly uroporphyrin, while protoporphyrin constituted the major porphyrin fraction in the feces of treated animals and in mouse bile. Urine of oxadiazon-treated rats contained coproporphyrin and uroporphyrin; total urine porphyrin content was 0.4 nmol/ml for control and 15.2 nmol/ml for treated rats. Porphobilinogen concentration was elevated in liver, kidneys, and blood plasma of treated animals. The activity of 5-aminolevulinic acid synthase was increased in the liver, but not in the kidneys of treated animals. Hepatic microsomal cytochrome P450 levels were decreased to about 50% of control values. The pattern of porphyrin accumulation and excretion in oxadiazon-induced experimental porphyria resembles the porphyrin pattern observed in human variegate porphyria.

The role of retrotransposons in gene family expansions: insights from the mouse Abp gene family

BackgroundRetrotransposons have been suggested to provide a substrate for non-allelic homologous recombination (NAHR) and thereby promote gene family expansion. Their precise role, however, is controversial. Here we ask whether retrotransposons contributed to the recent expansions of the Androgen-binding protein (Abp) gene families that occurred independently in the mouse and rat genomes.ResultsUsing dot plot analysis, we found that the most recent duplication in the Abp region of the mouse genome is flanked by L1Md_T elements. Analysis of the sequence of these elements revealed breakpoints that are the relicts of the recombination that caused the duplication, confirming that the duplication arose as a result of NAHR using L1 elements as substrates. L1 and ERVII retrotransposons are considerably denser in the Abp regions than in one Mb flanking regions, while other repeat types are depleted in the Abp regions compared to flanking regions. L1 retrotransposons preferentially accumulated in the Abp gene regions after lineage separation and roughly followed the pattern of Abp gene expansion. By contrast, the proportion of shared vs. lineage-specific ERVII repeats in the Abp region resembles the rest of the genome.ConclusionsWe confirmed the role of L1 repeats in Abp gene duplication with the identification of recombinant L1Md_T elements at the edges of the most recent mouse Abp gene duplication. High densities of L1 and ERVII repeats were found in the Abp gene region with abrupt transitions at the region boundaries, suggesting that their higher densities are tightly associated with Abp gene duplication. We observed that the major accumulation of L1 elements occurred after the split of the mouse and rat lineages and that there is a striking overlap between the timing of L1 accumulation and expansion of the Abp gene family in the mouse genome. Establishing a link between the accumulation of L1 elements and the expansion of the Abp gene family and identification of an NAHR-related breakpoint in the most recent duplication are the main contributions of our study.

Genomic islands of differentiation in two songbird species reveal candidate genes for hybrid female sterility

Hybrid sterility is a common first step in the evolution of postzygotic reproductive isolation. According to Haldanes Rule, it affects predominantly the heterogametic sex. While the genetic basis of hybrid male sterility in organisms with heterogametic males has been studied for decades, the genetic basis of hybrid female sterility in organisms with heterogametic females has received much less attention. We investigated the genetic basis of reproductive isolation in two closely related avian species, the common nightingale (Luscinia megarhynchos) and the thrush nightingale (L. luscinia), that hybridize in a secondary contact zone and produce viable hybrid progeny. In accordance with Haldanes Rule, hybrid females are sterile, while hybrid males are fertile, allowing gene flow to occur between the species. Using transcriptomic data from multiple individuals of both nightingale species, we identified genomic islands of high differentiation (FST) and of high divergence (Dxy), and we analysed gene content and patterns of molecular evolution within these islands. Interestingly, we found that these islands were enriched for genes related to female meiosis and metabolism. The islands of high differentiation and divergence were also characterized by higher levels of linkage disequilibrium than the rest of the genome in both species indicating that they might be situated in genomic regions of low recombination. This study provides one of the first insights into genetic basis of hybrid female sterility in organisms with heterogametic females.

The Role of Retrotransposons in Gene Family Expansions in the Human and Mouse Genomes

Abstract Retrotransposons comprise a large portion of mammalian genomes. They contribute to structural changes and more importantly to gene regulation. The expansion and diversification of gene families have been implicated as sources of evolutionary novelties. Given the roles retrotransposons play in genomes, their contribution to the evolution of gene families warrants further exploration. In this study, we found a significant association between two major retrotransposon classes, LINEs and LTRs, and lineage-specific gene family expansions in both the human and mouse genomes. The distribution and diversity differ between LINEs and LTRs, suggesting that each has a distinct involvement in gene family expansion. LTRs are associated with open chromatin sites surrounding the gene families, supporting their involvement in gene regulation, whereas LINEs may play a structural role promoting gene duplication. Our findings also suggest that gene family expansions, especially in the mouse genome, undergo two phases. The first phase is characterized by elevated deposition of LTRs and their utilization in reshaping gene regulatory networks. The second phase is characterized by rapid gene family expansion due to continuous accumulation of LINEs and it appears that, in some instances at least, this could become a runaway process. We provide an example in which this has happened and we present a simulation supporting the possibility of the runaway process. Altogether we provide evidence of the contribution of retrotransposons to the expansion and evolution of gene families. Our findings emphasize the putative importance of these elements in diversification and adaptation in the human and mouse lineages.

Repeat associated mechanisms of genome evolution and function revealed by the Mus caroli and Mus pahari genomes

Understanding the mechanisms driving lineage-specific evolution in both primates and rodents has been hindered by the lack of sister clades with a similar phylogenetic structure having high-quality genome assemblies. Here, we have created chromosome-level assemblies of the Mus caroli and Mus pahari genomes. Together with the Mus musculus and Rattus norvegicus genomes, this set of rodent genomes is similar in divergence times to the Hominidae (human-chimpanzee-gorilla-orangutan). By comparing the evolutionary dynamics between the Muridae and Hominidae, we identified punctate events of chromosome reshuffling that shaped the ancestral karyotype of Mus musculus and Mus caroli between 3 and 6 million yr ago, but that are absent in the Hominidae. Hominidae show between four- and sevenfold lower rates of nucleotide change and feature turnover in both neutral and functional sequences, suggesting an underlying coherence to the Muridae acceleration. Our system of matched, high-quality genome assemblies revealed how specific classes of repeats can play lineage-specific roles in related species. Recent LINE activity has remodeled protein-coding loci to a greater extent across the Muridae than the Hominidae, with functional consequences at the species level such as reproductive isolation. Furthermore, we charted a Muridae-specific retrotransposon expansion at unprecedented resolution, revealing how a single nucleotide mutation transformed a specific SINE element into an active CTCF binding site carrier specifically in Mus caroli, which resulted in thousands of novel, species-specific CTCF binding sites. Our results show that the comparison of matched phylogenetic sets of genomes will be an increasingly powerful strategy for understanding mammalian biology.

Alu insertion polymorphisms in the African Sahel and the origin of Fulani pastoralists

Abstract Background: The origin of Western African pastoralism, represented today by the Fulani nomads, has been a highly debated issue for the past decades, and has not yet been conclusively resolved. Aim: This study focused on Alu polymorphisms in sedentary and nomadic populations across the African Sahel to investigate patterns of diversity that can complement the existing results and contribute to resolving issues concerning the origin of West African pastoralism. Subjects and methods: A new dataset of 21 Alu biallelic markers covering a substantial part of the African Sahel has been analysed jointly with several published North African populations. Results: Interestingly, with regard to Alu variation, the relationship of Fulani pastoralists to North Africans is not as evident as was earlier revealed by studies of uniparental loci such as mtDNA and NRY. Alu insertions point rather to an affinity of Fulani pastoralists to Eastern Africans also leading a pastoral lifestyle. Conclusions: It is suggested that contemporary Fulani pastoralists might be descendants of an ancestral Eastern African population that, while crossing the Sahara in the Holocene, admixed slightly with a population of Eurasian (as evidenced by uniparental polymorphisms) ancestry. It seems that, in the Fulani pastoralists, Alu elements reflect more ancient genetic relationships than do uniparental genetic systems.

The roles of LINEs, LTRs and SINEs in lineage-specific gene family expansions in the human and mouse genomes

We explored genome-wide patterns of RT content surrounding lineage-specific gene family expansions in the human and mouse genomes. Our results suggest that the size of a gene family is an important predictor of the RT distribution in close proximity to the family members. The distribution differs considerably between the three most common RT classes (LINEs, LTRs and SINEs). LINEs and LTRs tend to be more abundant around genes of multi-copy gene families, whereas SINEs tend to be depleted around such genes. Detailed analysis of the distribution and diversity of LINEs and LTRs with respect to gene family size suggests that each has a distinct involvement in gene family expansion. LTRs are associated with open chromatin sites surrounding the gene families, supporting their involvement in gene regulation, whereas LINEs may play a structural role, promoting gene duplication. This suggests that gene family expansions, especially in the mouse genome, might undergo two phases, the first is characterized by elevated deposition of LTRs and their utilization in reshaping gene regulatory networks. The second phase is characterized by rapid gene family expansion due to continuous accumulation of LINEs and it appears that, in some instances at least, this could become a runaway process. We provide an example in which this has happened and we present a simulation supporting the possibility of the runaway process. Our observations also suggest that specific differences exist in this gene family expansion process between human and mouse genomes.