Immunoglobulin genes, reproductive isolation and vertebrate speciation

Blood transcriptome analysis in a buck-ewe hybrid points towards an nuclear factor-kappa B lymphoproliferative autoimmune disorder

Mammal hybridization is a speciation mechanism and an evolutionary driver. Goat-sheep, especially buck-ewe hybrids, are very rare with only one case reported in 2016, which is the subject of the work presented here. Blood transcriptome analysis revealed that the hybrid largely deviated from imprinting schemes previously described in sheep and other mammals. Furthermore, transcriptome regulation seems to differ from the parent transcriptomes, which is most likely a product of partially incompatible imprinting mechanisms from two closely related species. To gain a deeper understanding of hybridization in mammals we re-analyzed the RNA sequencing data of the buck-ewe hybrid and its parents. We found parent-of-origin-specific expression of genes that functionally clustered, which we explain with the Dobzhansky-Muller incompatibility (DMI) model. According to the DMI model, proteins which interact have a high probability of being barrier loci and hence are prone to monoallelic expression. We discovered enrichment of genes uniquely expressed by the buck-ewe hybrid, which implicate that it suffered from an NF-κB lymphoproliferative autoimmune disorder. Similar findings were reported in the F1 generation of hybrid mice. We propose that hybridization of two related species may lead to an autoimmune phenotype, due to immunoglobulin incompatibilities and incomplete silencing of barrier loci.

Adaptive immune receptor genotyping using the corecount program

We present a new Rep-Seq analysis tool called corecount, for analyzing genotypic variation in immunoglobulin (IG) and T cell receptor (TCR) genes. corecount is highly efficient at identifying V alleles, including those that are infrequently used in expressed repertoires and those that contain 3' end variation that are otherwise refractory to reliable identification during germline inference from expressed libraries. Furthermore, corecount facilitates accurate D and J gene genotyping. The output is highly reproducible and facilitates the comparison of genotypes from multiple individuals, such as those from clinical cohorts. Here, we applied corecount to the genotypic analysis of IgM libraries from 16 individuals. To demonstrate the accuracy of corecount, we Sanger sequenced all the heavy chain IG alleles (65 IGHV, 27 IGHD and 7 IGHJ) from one individual from whom we also produced two independent IgM Rep-seq datasets. Genomic analysis revealed that 5 known IGHV and 2 IGHJ sequences are truncated in current reference databases. This dataset of genomically validated alleles and IgM libraries from the same individual provides a useful resource for benchmarking other bioinformatic programs that involve V, D and J assignments and germline inference, and may facilitate the development of AIRR-Seq analysis tools that can take benefit from the availability of more comprehensive reference databases.

Analysis of the Whole-Genome Sequences from an Equus Parent-Offspring Trio Provides Insight into the Genomic Incompatibilities in the Hybrid Mule

Interspecific hybridization often shows negative effects on hybrids. However, only a few multicellular species, limited to a handful of plants and animals, have shown partial genetic mechanisms by which hybridization leads to low fitness in hybrids. Here, to explore the outcome of combining the two genomes of a horse and donkey, we analyzed the whole-genome sequences from an Equus parent-offspring trio using Illumina platforms. We generated 41.39× and 46.21× coverage sequences for the horse and mule, respectively. For the donkey, a 40.38× coverage sequence was generated and stored in our laboratory. Approximately 24.86 million alleles were discovered that varied from the reference genome. Single nucleotide polymorphisms were used as polymorphic markers for assigning alleles to their parental genomic inheritance. We identified 25,703 Mendelian inheritance error single nucleotide polymorphisms in the mule genome that were not inherited from the parents through Mendelian inheritance. A total of 555 de novo single nucleotide polymorphisms were also identified. The rate of de novo single nucleotide polymorphisms was 2.21 × 10^-7 in the mule from the Equus parent-offspring trio. This rate is obviously higher than the natural mutation rate for Equus, which is also consistent with the previous hypothesis that interracial crosses may have a high mutation rate. The genes associated with these single nucleotide polymorphisms are mainly involved in immune processes, DNA repair, and cancer processes. The results of the analysis of three genomes from an Equus parent-offspring trio improved our knowledge of the consequences of the integration of parental genomes in mules.