Inbreeding and purging at the genomic Level: the Chillingham cattle reveal extensive,non‐random SNP heterozygosity

Local breeds of livestock are of conservation significance as components of global biodiversity and as reservoirs of genetic variation relevant to the future sustainability of agriculture. One such rare historic breed, the Chillingham cattle of northern England, has a 350‐year history of isolation and inbreeding yet shows no diminution of viability or fertility. The Chillingham cattle have not been subjected to selective breeding. It has been suggested previously that the herd has minimal genetic variation. In this study, high‐density SNP genotyping with the 777K SNP chip showed that 9.1% of loci on the chip are polymorphic in the herd, compared with 62–90% seen in commercial cattle breeds. Instead of being homogeneously distributed along the genome, these loci are clustered at specific chromosomal locations. A high proportion of the Chillingham individuals examined were heterozygous at many of these polymorphic loci, suggesting that some loci are under balancing selection. Some of these frequently heterozygous loci have been implicated as sites of recessive lethal mutations in cattle. Linkage disequilibrium equal or close to 100% was found to span up to 1350 kb, and LD was above r² = 0.25 up to more than 5000 kb. This strong LD is consistent with the lack of polymorphic loci in the herd. The heterozygous regions in the Chillingham cattle may be the locations of genes relevant to fitness or survival, which may help elucidate the biology of local adaptation in traditional breeds and facilitate selection for such traits in commercial cattle.     

Exploring the unbinding of Leishmania (L.) amazonensis CPB derived‐epitopes from H2 MHC class I proteins

New strategies to control Leishmania disease demand an extensive knowledge about several aspects of infection including the understanding of its molecular events. In murine models, cysteine proteinase B from Leishmania amazonensis promotes regulation of immune response, and fragments from its C‐terminus extension (cyspep) can play a decisive role in the host‐parasite interaction. The interaction between cyspep‐derived peptides and major histocompatibility complex (MHC) proteins is a crucial factor in Leishmania infections. Seven cyspep‐derived peptides, previously identified as capable of interacting with H‐2 (murine) MHC class I proteins, were studied in this work. We established a protocol to simulate the unbinding of these peptides from the cleft of H‐2 receptors. From the simulations, we estimated the corresponding free energy of dissociation (ΔG_d) and described the molecular events that occur during the exit of peptides from the cleft. To test the reliability of this method, we first applied it to a calibration set of four crystallographic MHC/peptide complexes. Next, we explored the unbinding of the seven complexes mentioned above. Results were consistent with ΔG_d values obtained from surface plasmon resonance (SPR) experiments. We also identified some of the primary interactions between peptides and H‐2 receptors, and we detected three regions of influence for the interaction. This pattern was systematically observed for the peptides and helped determine a minimum distance for the real interaction between peptides and H‐2 proteins occurring at ～25 Å. Proteins 2016; 84:473–487. © 2016 Wiley Periodicals, Inc.     

Mitochondrial DNA diversity and phylogeography of endangered green turtle (Chelonia mydas) populations in Africa

We analysed the genetic structure of seven nesting sites of the endangered green turtle (Chelonia mydas) in Africa using mitochondrial DNA control region sequences. Tissue samples were collected from 188 nesting females at six sites in West Africa and one in the Indian Ocean. A 488 bp fragment of the control region revealed 14 different haplotypes, 10 of which are previously undescribed. The most common haplotype (CM8) was observed in 157 individuals. All other haplotypes were closely related, except two divergent lineages: CM38, removed by four substitutions, and the three Indian Ocean haplotypes, distinguished by 31 substitutions. Significant differences in haplotype and nucleotide diversity were observed between Atlantic rookeries and among ocean basins. Analysis of molecular variance revealed high levels of differentiation between the Atlantic and the Indian Ocean populations but a much shallower Atlantic substructuring. Green turtle population genetic structure is thought to have been shaped by a dynamic succession of extinction and recolonisation of rookeries, by natal homing and occasional breakdown in nest-site fidelity. Mismatch distributions of pairwise differences between haplotypes at each rookery were found to be consistent with recent population expansion. We argue that demographic histories can be explained by scenarios at several temporal scales, including geological events, sea level fluctuations and more recent patterns of exploitation. We discuss management and conservation implications of our results for these threatened populations, identifying two ESUs (one in the Atlantic and one in the Indian ocean) and three MUs within the Atlantic.     

Mitochondrial DNA polymorphisms in Phytophthora infestans: New haplotypes are identified and re-defined by PCR

Polymorphisms of mitochondrial DNA (mt-DNA) are particularly useful for monitoring specific pathogen populations like Phytophthora infestans. Basically type I and II of P. infestans mt-DNA were categorized by means of polymorphism lengths caused by an ~ 2 kb insertion, which can be detected via restriction enzyme digestion. In addition genome sequencing of haplotype Ib has been used as a simple Polymerase Chain Reaction–Restriction Fragment Length Polymorphism (PCR–RFLP) method to indirectly identify type I and II alterations through EcoR I restriction enzyme DNA fragment patterns of the genomic P4 area. However, with the common method, wrong mt-DNA typing occurs due to an EcoR I recognition site mutation in the P4 genomic area. Genome sequencing of the four haplotypes (Ia, Ib, IIa, and IIb) allowed us to thoroughly examine mt-DNA polymorphisms and we indentified two hypervariable regions (HVRs) named HVRi and HVRii. The HVRi length polymorphism caused by a 2 kb insertion/deletion was utilized to identify mt-DNA types I and II, while another length polymorphism in the HVRii region is caused by a variable number of tandem repeats (n = 1, 2, or 3) of a 36 bp sized DNA stretch and was further used to determine mt-DNA sub-types, which were described as R_{n=1, 2, or 3}. Finally, the P. infestans mt-DNA haplotypes were re-defined as IR₁ or IIR₂ according to PCR derived HVRi and HVRii length polymorphisms. Twenty-three isolates were chosen to verify the feasibility of our new approach for identifying mt-DNA haplotypes and a total of five haplotypes (IR₁, IR₂, IR₃, IIR₂ and IIR₃) were identified. Additionally, we found that six isolates determined as type I by our method were mistakenly identified as type II by the PCR–RFLP technique. In conclusion, we propose a simple and rapid PCR method for identification of mt-DNA haplotypes based on sequence analyses of the mitochondrial P. infestans genome.     

Genetic and morphological variation in two littorinid gastropods: evidence for recent population expansions along the East African coast

In marine species, population diversity and differentiation is affected by the population history and by the complex interaction between oceanographic dynamics and ecological traits. In the present study, we examined two species of marine gastropods (the mangrove periwinkle Littoraria scabra and the rocky shore Littoraria glabrata) along the East African coast, using both genetic and geometric morphometric methods. We report a greater variation of shell shape in L. scabra compared to the slightly smaller variation in L. glabrata. This variation was probably associated with variation of environmental factors along the coast, such as temperature and hydrodynamics. Despite morphological variation, we found low mitochondrial genetic differentiation among samples from different localities for both species, which is probably a consequence of the ongoing gene flow during the free‐swimming larval stage of these gastropods. Additionally, high levels of haplotype diversity, low nucleotide diversity, and ‘star‐like’ genealogies were found in both species. These observations and the results from mismatch distributions, indicate a possible signature of recent population expansions in both species, which probably started during interglacial periods of the Pleistocene and led to the colonization of the Indian Ocean coast. © 2013 The Linnean Society of London     

Genetic discontinuity revealed by chloroplast microsatellites in eastern North American Abies (Pinaceae)

Development of conservation strategies for Fraser fir (Abies fraseri) in the southern Appalachian Mountains depends in part on recognition of the extent to which Fraser fir is genetically distinct from the closely related balsam (A. balsamea) and intermediate (A. balsamea var. phanerolepis) fir. These sibling species have exhibited intergrading, clinal variation in morphological, chemical, and genetic characteristics in prior research. Chloroplast microsatellite markers were polymerase chain reaction amplified from genomic DNA samples of 78 individuals representing the geographic ranges of Fraser, balsam, and intermediate fir. Gene diversity levels at two loci ranged among taxa from 0.65 to 0.84. Allele frequencies demonstrated significant differentiation among taxa, with R(ST) values of 0.36 and 0.10. Haplotype diversity and D(SH) were highest for balsam fir and lowest for intermediate fir. A haplotype network analysis based on allele size distribution for the two loci revealed two distinct clusters of haplotypes and population-specific haplotypes. Ninety-two percent of the haplotypes in one cluster were from balsam fir and intermediate fir, and 84% of the haplotypes in the other cluster were from Fraser fir and intermediate fir. The genetic differentiation of chloroplast DNA markers provides justification for the recognition of Fraser fir as a distinct Management Unit (MU) for conservation purposes, regardless of its taxonomic classification.     

COMPETITION BETWEEN SEGREGATION DISTORTERS: COEXISTENCE OF “SUPERIOR” AND “INFERIOR” HAPLOTYPES AT THE t COMPLEX

By means of population genetical models, we investigate the competition between sex-specific segregation distorters. Although the models are quite general, they are motivated by a specific example, the t complex of the house mouse. Some variants at this gene complex, the t haplotypes, distort Mendelian segregation in heterozygous males in their favor. The selective advantage at the gamete level is counterbalanced by strong negative fitness effects at the individual level (male sterility or even lethality in both sexes). A plethora of different t haplotypes has been found, both in the field and in the lab. Up to now, however, models have focused on the equilibrium frequency of a single t haplotype. In contrast, we explicitly model the competition between several t haplotypes. A deterministic model for a large, well-mixed population predicts a surprisingly high degree of polymorphism. Haplotypes with seemingly inferior fitness characteristics may easily coexist with “superior” haplotypes. For instance, a lethal haplotype with a low segregation ratio may stably coexist with a sterile haplotype with a high segregation ratio. Stable coexistence is even possible for haplotypes with a segregation disadvantage. A simple stochastic model shows that the same principles apply in the context of a structured metapopulation. Although counterintuitive at first sight, all our results can be explained by the fact that segregation distorters have an inherent advantage when they are rare. We conclude that fitness comparisons are not sufficient to predict the outcome of competition when selective forces are acting at different levels.     

Novel alpha haemoglobin haplotypes in horses

Four minor haplotypes that produce abnormal haemoglobin phenotypes in horses have been characterized. Two of them, AIIb and V, are copy number variants with, respectively, one and three alpha genes instead of the normal complement of two. The AIIa and C haplotypes, on the other hand, each have two alpha genes but, as a result of probable gene conversions, they now encode identical, though haplotype specific, globins. Two out of 60 unrelated and phenotypically normal horses studied had an unusual triplicated rearrangement in the embryonic zeta-gene locus. Each of these variants appears to have been produced by aberrant recombination events.     

Population genetic structure of skipjack tuna Katsuwonus pelamis from the Indian coast using sequence analysis of the mitochondrial DNA D-loop region

Genetic structure of skipjack tuna Katsuwonus pelamis from the Indian region was investigated using sequence data of mitochondrial DNA (mtDNA) D-loop region. A total of 315 individuals were sampled from six major fishing grounds around the east and west coasts of India including the Andaman (Port Blair) and Lakshadweep (Minicoy) Islands. Nucleotide and gene diversities were high in all the sample collections. Significant genetic heterogeneity was observed for the mtDNA sequence data among sites (φ(ST) = 0·0273, P < 0·001). Analysis of molecular variance (AMOVA) showed significant genetic variation among four groups (φ(CT) = 0·0261, P < 0·05) which was also supported by spatial AMOVA results. The null hypothesis of single panmictic population of K. pelamis along the Indian coast can thus be rejected. Phylogenetic analysis of the mtDNA sequence data showed the presence of four clades of K. pelamis in the Indian waters. There was no clear pattern, however, of haplotypes and geographic location among samples. The results of this study suggest the occurrence of four genetically differentiated groups of K. pelamis across the coastal waters of India.     

Taiwanese aborigines: genetic heterogeneity and paternal contribution to Oceania

In the present study, for the first time, 293 Taiwanese aboriginal males from all nine major tribes (Ami, Atayal, Bunun, Rukai, Paiwan, Saisat, Puyuma, Tsou, Yami) were genotyped with 17 YSTR loci in a attend to reveal migrational patterns connected with the Austronesian expansion. We investigate the paternal genetic relationships of these Taiwanese aborigines to 42 Asia-Pacific reference populations, geographically selected to reflect various locations within the Austronesian domain. The Tsou and Puyuma tribes exhibit the lowest (0.1851) and the highest (0.5453) average total genetic diversity, respectively. Further, the fraction of unique haplotypes is also relatively high in the Puyuma (86.7%) and low in Tsou (33.3%) suggesting different demographic histories. Multidimensional scaling (MDS) and analysis of molecular variance (AMOVA) revealed several notable findings: 1) the Taiwan indigenous populations are highly diverse. In fact, the level of inter-population heterogeneity displayed by the Taiwanese aboriginal populations is close to that exhibited among all 51 Asia-Pacific populations examined; 2) the asymmetrical contribution of the Taiwanese aborigines to the Oceanic groups. Ami, Bunun and Saisiyat tribes exhibit the strongest paternal links to the Solomon and Polynesian island communities, whereas most of the remaining Taiwanese aboriginal groups are more genetically distant to these Oceanic inhabitants; 3) the present YSTR analyses does not reveal a strong paternal affinity of the nine Taiwanese tribes to their continental Asian neighbors. Overall, our current findings suggest that, perhaps, only a few of the tribes were involved in the migration out of Taiwan.