The genomic structure of isolation across breed, country and strain for important South African and Australian sheep populations

Background South Africa and Australia shares multiple important sheep breeds. For some of these breeds, genomic breeding values are provided to breeders in Australia, but not yet in South Africa. Combining genomic resources could facilitate development for across country selection, but the influence of population structures could be important to the compatability of genomic data from varying origins. The genetic structure within and across breeds, countries and strains was evaluated in this study by population genomic parameters derived from SNP-marker data. Populations were first analysed by breed and country of origin and then by subpopulations of South African and Australian Merinos. Results Mean estimated relatedness according to the genomic relationship matrix varied by breed (-0.11 to 0.16) and bloodline (-0.08 to 0.06) groups and depended on co-ancestry as well as recent genetic links. Measures of divergence across bloodlines (FST: 0.04–0.12) were sometimes more distant than across some breeds (FST: 0.05–0.24), but the divergence of common breeds from their across-country equivalents was weak (FST: 0.01–0.04). According to mean relatedness, FST, PCA and Admixture, the Australian Ultrafine line was better connected to the SA Cradock Fine Wool flock than with other AUS bloodlines. Levels of linkage disequilibrium (LD) between adjacent markers was generally low, but also varied across breeds (r2: 0.14–0.22) as well as bloodlines (r2: 0.15–0.19). Patterns of LD decay was also unique to breeds, but bloodlines differed only at the absolute level. Estimates of effective population size (Ne) showed genetic diversity to be high for the majority of breeds (Ne: 128–418) but also for bloodlines (Ne: 137–369). Conclusions This study reinforced the genetic complexity and diversity of important sheep breeds, especially the Merino breed. The results also showed that implications of isolation can be highly variable and extended beyond breed structures. However, knowledge of useful links across these population substructures allows for a fine-tuned approach in the combination of genomic resources. Isolation across country rarely proved restricting compared to other structures considered. Consequently, research into the accuracy of across-country genomic prediction is recommended.

Coronavirus Genome Sequence Similarity and Protein Sequence Classification

The world currently is going through a serious pandemic due to the coronavirus disease (COVID-19). In this study, we investigate the gene structure similarity of coronavirus genomes isolated from COVID-19 patients, Severe Acute Respiratory Syndrome (SARS) patients and bats genes. We also explore the extent of similarity between their genome structures to find if the new coronavirus is similar to either of the other genome structures. Our experimental results show that there is 82.42% similarity between the CoV-2 genome structure and the bat genome structure. Moreover, we have used a bidirectional Gated Recurrent Unit (GRU) model as the deep learning technique and an improved variant of Recurrent Neural networks (i.e., Bidirectional Long Short Term Memory model) to classify the protein families of these genomes to isolate the prominent protein family accession. The accuracy of Gated Recurrent Unit (GRU) is 98% for labeled protein sequences against the protein families. By comparing the performance of the Gated Recurrent Unit (GRU) model with the Bidirectional Long Short Term Memory (Bi-LSTM) model results, we found that the GRU model is 1.6% more accurate than the Bi-LSTM model for our multiclass protein classification problem. Our experimental results would be further support medical research purposes in targeting the protein family similarity to better understand the coronavirus genomic structure.

The Organization of the Pig T-cell Receptor Gamma (TRG) Locus Provides Insights Into the Evolutionary Patterns of the TRG Genes Across Cetartiodactyla

The domestic pig (Sus scrofa) is a species representative of the Suina, one of the four suborders within Cetartiodactyla. In this paper, we reported our analysis of the pig TRG locus in comparison with the loci of species representative of the Ruminantia, Tylopoda and Cetacea suborders. The pig TRG genomic structure reiterates the peculiarity of the organization of Cetartiodactyla loci in TRGC “cassettes”, each containing the basic V-J-J-C unit. Eighteen genes arranged in four TRGC cassettes, form the pig TRG locus. All the functional TRG genes were expressed, and the TRGV genes preferentially rearrange with the TRGJ genes within their own cassette, which correlates the diversity of the gamma-chain repertoire with the number of cassettes. Among them, the TRGC5, located at the 5’ end of the locus, is the only cassette that retains a marked homology with the corresponding TRGC cassettes of all the analyzed species. The preservation of the TRGC5 cassette for such a long evolutionary time presumes a highly specialized function of its genes, which could be essential for the survival of species. Therefore, the maintenance of this cassette in pigs confirms that it is the most evolutionarily ancient within Cetartiodactyla, and it has undergone a process of duplication to give rise to the other TRGC cassettes in the different artiodactyl species in a lineage-specific manner.

Genomics Reveal Admixture and Unexpected Patterns of Diversity in a Parapatric Pair of Butterflies

We studied the evolutionary relationship of two widely distributed parapatric butterfly species, Melitaea athalia and Melitaea celadussa, using the ddRAD sequencing approach, as well as genital morphology and mtDNA data. M. athalia was retrieved as paraphyletic with respect to M. celadussa. Several cases of mito-nuclear discordance and morpho-genetic mismatch were found in the contact zone. A strongly diverged and marginally sympatric clade of M. athalia from the Balkans was revealed. An in-depth analysis of genomic structure detected high levels of admixture between M. athalia and M. celadussa at the contact zone, though not reaching the Balkan clade. The demographic modelling of populations supported the intermediate genetic make-up of European M. athalia populations with regards to M. celadussa and the Balkan clade. However, the dissimilarity matrix of genotype data (PCoA) suggested the Balkan lineage having a genetic component that is unrelated to the athalia-celadussa group. Although narrowly sympatric, almost no signs of gene flow were found between the main M. athalia group and the Balkan clade. We propose two possible scenarios on the historical evolution of our model taxa and the role of the last glacial maximum in shaping their current distribution. Finally, we discuss the complexities regarding the taxonomic delimitation of parapatric taxa.

Population genomic structure in Goodman's mouse lemur reveals long-standing separation of Madagascar's Central Highlands and eastern rainforests

Madagascar’s Central Highlands are largely composed of grasslands, interspersed with patches of forest. The pre-human extent of these grasslands is a topic of vigorous debate, with conventional wisdom holding that they are anthropogenic in nature and emerging evidence supporting that grasslands were a component of the pre-human Central Highlands vegetation. Here, we shed light on the temporal dynamics of Madagascar’s vegetative composition by conducting a population genomic investigation of Goodman’s mouse lemur (Microcebus lehilahytsara; Cheirogaleidae). These small-bodied primates occur both in Madagascar’s eastern rainforests and in the Central Highlands, which makes them a valuable indicator species. Population divergences among forest-dwelling mammals can serve as a proxy for habitat fragmentation and patterns of post-divergence gene flow can reveal potential migration corridors consistent with a wooded grassland mosiac. We used RADseq data to infer phylogenetic relationships, population structure, demographic models of post-divergence gene flow, and population size change through time. These analyses offer evidence that open habitats are an ancient component of the Central Highlands, and that wide-spread forest fragmentation occurred naturally during a period of decreased precipitation near the last glacial maximum. Models of gene flow suggest that migration across the Central Highlands has been possible from the Pleistocene through the recent Holocene via riparian corridors. Notably, though our findings support the hypothesis that Central Highland grasslands predate human arrival, we also find evidence for human-mediated population declines. This highlights the extent to which species imminently threatened by human-mediated deforestation may be more vulnerable from paleoclimatic changes.

Unique roles of vaginal Megasphaera phylotypes in reproductive health

The composition of the human vaginal microbiome has been extensively studied and is known to influence reproductive health. However, the functional roles of individual taxa and their contributions to negative health outcomes have yet to be well characterized. Here, we examine two vaginal bacterial taxa grouped within the genus Megasphaera that have been previously associated with bacterial vaginosis (BV) and pregnancy complications. Phylogenetic analyses support the classification of these taxa as two distinct species. These two phylotypes, Megasphaera phylotype 1 (MP1) and Megasphaera phylotype 2 (MP2), differ in genomic structure and metabolic potential, suggestive of differential roles within the vaginal environment. Further, these vaginal taxa show evidence of genome reduction and changes in DNA base composition, which may be common features of host dependence and/or adaptation to the vaginal environment. In a cohort of 3870 women, we observed that MP1 has a stronger positive association with bacterial vaginosis whereas MP2 was positively associated with trichomoniasis. MP1, in contrast to MP2 and other common BV-associated organisms, was not significantly excluded in pregnancy. In a cohort of 52 pregnant women, MP1 was both present and transcriptionally active in 75.4 % of vaginal samples. Conversely, MP2 was largely absent in the pregnant cohort. This study provides insight into the evolutionary history, genomic potential and predicted functional role of two clinically relevant vaginal microbial taxa.