Genome-wide Single Nucleotide Polymorphism Data Research Articles

Drivers of genomic diversity and phenotypic development in early phases of domestication in Hermetia illucens.

The black soldier fly (BSF), Hermetia illucens, has the ability to efficiently bioremediate organic waste into usable bio-compounds. Understanding the impact of domestication and mass rearing on fitness and production traits is therefore important for sustainable production. This study aimed to assess patterns of genomic diversity and its association to phenotypic development across early generations of mass rearing under two selection strategies: selection for greater larval mass (SEL lines) and no direct artificial selection (NS lines). Genome-wide single nucleotide polymorphism (SNP) data were generated using 2bRAD sequencing, while phenotypic traits relating to production and population fitness were measured. Declining patterns of genomic diversity were observed across three generations of captive breeding, with the lowest diversity recorded for the F3 generation of both selection lines, most likely due to founder effects. The SEL cohort displayed statistically significantly greater larval weight com the NS lines with pronounced genetic and phenotypic directional changes across generations. Furthermore, lower genetic and phenotypic diversity, particularly for fitness traits, were evident for SEL lines, illustrating the trade-off between selecting for mass and the resulting decline in population fitness. SNP-based heritability was significant for growth, but was low or non-significant for fitness traits. Genotype-phenotype correlations were observed for traits, but individual locus effect sizes where small and very few of these loci demonstrated a signature for selection. Pronounced genetic drift, due to small effective population sizes, is likely overshadowing the impacts of selection on genomic diversity and consequently phenotypic development. The results hold particular relevance for genetic management and selective breeding for BSF in future.

O-271 Two component resolution of mosaic samples in PGT-A: primary template amplification (PTA) and independent analysis using copy-number variation and allelic balance

Abstract Study question How does accurate, unbiased whole genome amplification and tiered, independent, two component analysis impact PGT-A results? Summary answer When compared to previous methods a PGT-A workflow utilizing PTA followed by a CNV and allelic balance-based analysis generated more conclusive results. What is known already Accurate PGT-A results require unbiased error-free whole genome amplification from the trophectoderm biopsy and smart bioinformatics utilities to fully analyze and consider the large complex datasets. PTA generates fewer errors, greater coverage, and reduced allele drop-out offering the best starting point to interpret the genome of each embryo. Current PGT-A platforms leverage either normalized CNV data OR allelic balance at single nucleotide polymorphisms (SNPs), to detect alterations in the embryonic genome. A pipeline leveraging improved amplification technology and consolidated CNV and SNP data, particularly for putative mosaic and segmental mosaic calls, could reduce uncertainty and provide more accurate PGT-A results. Study design, size, duration A retrospective analysis of PGT-A results was undertaken with samples originating from multiple clinics and processed in a single laboratory, from January through December, 2023. Approximately 500 samples were processed simultaneously through two independent PGT-A pipelines. The primary objective was to compare the performance of CNV-only PGT-A with a test leveraging both CNV and SNP data. A secondary objective was to observe the resolution of mosaic embryos. Participants/materials, setting, methods De-identified samples were amplified with PTA and sequenced with both shotgun (coverage-based CNV) and capture (SNP-based)-approaches. Method specific analysis was applied to the resulting data. An initial bioinformatics approach analyzed CNV data only for variations in copy-number along the genome, proposing stretches of aneuploidy or mosaicism. A second, tiered analysis extended this data and incorporated SNP-based allelic balance to evaluate candidate regions and confirm the presence of aneuploidy and/or quantify mosaic levels. Main results and the role of chance Genome-wide read coverage and SNP data confirmed that PTA accurately and completely amplified the genetic material in each sample with broad distribution of informative data and a low allele drop-out rate. High concordance (95%) was observed for euploid and full aneuploid (whole chromosome and segmental) embryos classified with CNV-alone and CNV/SNV combined analysis. SNP data was processed, and B-allele ratios were used to confirm aneuploidies with particular utility in resolving the location of segmental changes. The classification of mosaic samples was altered in the pipeline leveraging both CNV and SNP data. Of 54 mosaic samples tested, 75% remained mosaic after independent consideration of SNP data, ∼15% and 6% of high-level mosaics were resolved to full aneuploid and euploid, respectively, and ∼25% of low-level mosaics reverted to a euploid classification. The deep SNP data produced during targeted sequencing confirmed the altered allele ratio predicted by CNV data and re-confirmed the presence and proportion of aneuploid cells. Mosaic levels were highly correlated in 90% of re-confirmed mosaic samples. This suggests that the classification of these samples as mosaic is precise and accurate and suggests a strong utility for a two component PGT-A pipeline leveraging CNV and SNV data generated from PTA embryos. Limitations, reasons for caution This was a retrospective study with limited sample numbers. Wider implications of the findings The study highlights the strong advantages of a hybrid PGT-A pipeline, leveraging PTA and two component analysis of coverage and allelic balance, particularly to resolve mosaic embryos more accurately. Overall, this could reduce uncertainty and deliver more confident PGT-A results for embryo transfer decision making. Trial registration number not applicable

On the road to losing connectivity: Faecal samples provide genome-wide insights into anthropogenic impacts on two large herbivore species in central India.

Humans have impacted most of the planet, and the ensuing fragmentation results in small, isolated habitat patches posing a risk of genetic diversity loss, inbreeding, and genetic load. Understanding how natural and anthropogenic landscape features affect gene flow among habitat patches is critical for maintaining connectivity. Genome-wide data are required to comprehend the impacts of recent fragmentation, which can be challenging when only non-invasive samples are available. Here, we build upon advancements in conservation genomics to address connectivity of two large herbivores, gaur (Bos gaurus) and sambar (Rusa unicolor) in central India. Given their habitat associations, we expected these species to respond similarly to habitat fragmentation. We used faecal-DNA and methylation-based host-DNA enrichment with modified ddRAD protocol to generate genome-wide single-nucleotide polymorphism (SNP) data for 124 gaur and 99 sambar individuals. Our findings reveal that gaur populations in central India are fragmented, displaying high genetic differentiation, with drift significantly affecting small populations like Umred Karhandla Wildlife Sanctuary. Although sambar shows low genetic structure, another small population, Bor Tiger Reserve is genetically differentiated. Our results suggest that although land cover change and roads restrict animal movement, the extent of this impact varies across the two species. We show that different species respond differently to landscape features, even with similar habitat associations. We highlight small and isolated populations requiring urgent conservation intervention. Such multi-species approaches enhance our understanding of cross-species connectivity patterns. We suggest shifting from single-species to multi-species holistic conservation approach in rapidly developing landscapes to better manage co-occurring endangered species.

Shared Genetic Determinants of Schizophrenia and Autism Spectrum Disorder Implicate Opposite Risk Patterns: A Genome-Wide Analysis of Common Variants.

The synaptic pruning hypothesis posits that schizophrenia (SCZ) and autism spectrum disorder (ASD) may represent opposite ends of neurodevelopmental disorders: individuals with ASD exhibit an overabundance of synapses and connections while SCZ was characterized by excessive pruning of synapses and a reduction. Given the strong genetic predisposition of both disorders, we propose a shared genetic component, with certain loci having differential regulatory impacts. Genome-Wide single nucleotide polymorphism (SNP) data of European descent from SCZ (N cases = 53386, N controls = 77258) and ASD (N cases = 18381, N controls = 27969) were analyzed. We used genetic correlation, bivariate causal mixture model, conditional false discovery rate method, colocalization, Transcriptome-Wide Association Study (TWAS), and Phenome-Wide Association Study (PheWAS) to investigate the genetic overlap and gene expression pattern. We found a positive genetic correlation between SCZ and ASD (rg = .26, SE = 0.01, P = 7.87e-14), with 11 genomic loci jointly influencing both conditions (conjFDR <0.05). Functional analysis highlights a significant enrichment of shared genes during early to mid-fetal developmental stages. A notable genetic region on chromosome 17q21.31 (lead SNP rs2696609) showed strong evidence of colocalization (PP.H4.abf = 0.85). This SNP rs2696609 is linked to many imaging-derived brain phenotypes. TWAS indicated opposing gene expression patterns (primarily pseudogenes and long noncoding RNAs [lncRNAs]) for ASD and SCZ in the 17q21.31 region and some genes (LRRC37A4P, LINC02210, and DND1P1) exhibit considerable variation in the cerebellum across the lifespan. Our findings support a shared genetic basis for SCZ and ASD. A common genetic variant, rs2696609, located in the Chr17q21.31 locus, may exert differential risk regulation on SCZ and ASD by altering brain structure. Future studies should focus on the role of pseudogenes, lncRNAs, and cerebellum in synaptic pruning and neurodevelopmental disorders.

The role of historical biogeography in shaping colour morph diversity in the common wall lizard.

The maintenance of polymorphisms often depends on multiple selective forces, but less is known on the role of stochastic or historical processes in maintaining variation. The common wall lizard (Podarcis muralis) is a colour polymorphic species in which local colour morph frequencies are thought to be modulated by natural and sexual selection. Here, we used genome-wide single-nucleotide polymorphism data to investigate the relationships between morph composition and population biogeography at a regional scale, by comparing morph composition with patterns of genetic variation of 54 populations sampled across the Pyrenees. We found that genetic divergence was explained by geographic distance but not by environmental features. Differences in morph composition were associated with genetic and environmental differentiation, as well as differences in sex ratio. Thus, variation in colour morph frequencies could have arisen via historical events and/or differences in the permeability to gene flow, possibly shaped by the complex topography and environment. In agreement with this hypothesis, colour morph diversity was positively correlated with genetic diversity and rates of gene flow and inversely correlated with the likelihood of the occurrence of bottlenecks. Concurrently, we did not find conclusive evidence for selection in the two colour loci. As an illustration of these effects, we observed that populations with higher proportions of the rarer yellow and yellow-orange morphs had higher genetic diversity. Our results suggest that processes involving a decay in overall genetic diversity, such as reduced gene flow and/or bottleneck events have an important role in shaping population-specific morph composition via non-selective processes.

Phylogenetic relationships of Quercus species (Fagaceae) in Vietnam based on multiplexed inter simple sequence repeat genotyping by sequencing

Quercus is the largest genus of the family Fagaceae in the world and the second-largest in Vietnam, with approximately 50 species. However, the phylogeny study of the Quercus genus in Vietnam has yet to be thoroughly explored. In this study, we utilize the genome-wide single-nucleotide polymorphisms (SNPs) data obtained through Multiplexed Inter Simple Sequence Repeat Genotyping by sequencing (MIG-seq) to explore the phylogenetic relationships among Quercus species in Vietnam. The results of this study reveal that all Quercus species in Vietnam belong to subgenus Cerris and the phylogenetic analysis strongly supports the recognition of two infrageneric sections: Quercus and Ilex section for the Vietnamese Quercus. These results also confirm that the monophyly of Quercus in Vietnam is strongly supported by both morphological and molecular data. The results of this study also align entirely with previous research, indicating that Quercus species in Asia fall under the subgenus Cerris, comprising two sections: section Ilex and section Cyclobalanopsis.

Genomic integrity of Phyciodes butterfly species in a region of contact (Lepidoptera: Nymphalidae)

Abstract Crescent butterflies of the Phyciodes tharos group have a long-standing reputation for taxonomic difficulty. We assessed species boundaries in the P. tharos group using genome-wide single-nucleotide polymorphism (SNP) data, focusing on sampling in western Canada, where 4 species (P. tharos, P. cocyta, P. batesii, and P. pulchella) have classically been recognized. Phylogenetic and cluster-based analyses confirm that there are indeed 4 species based on the maintenance of genomic integrity in the presence of occasional hybridization and low levels of gene flow. Mitochondrial cytochrome c oxidase I (COI) does not consistently distinguish species, with haplotypes extensively shared between species. Limited sampling of P. batesii indicates complex boundaries and interactions with other Phyciodes, especially with P. pulchella and P. cocyta. This study demonstrates the power of a genome-wide perspective to elucidate the integrity of boundaries between recently diverged species that remain capable of occasional hybridization.

Genomics Research of Lifetime Depression in the Netherlands: The BIObanks Netherlands Internet Collaboration (BIONIC) Project.

In this cohort profile article we describe the lifetime major depressive disorder (MDD) database that has been established as part of the BIObanks Netherlands Internet Collaboration (BIONIC). Across the Netherlands we collected data on Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) lifetime MDD diagnosis in 132,850 Dutch individuals. Currently, N = 66,684 of these also have genomewide single nucleotide polymorphism (SNP) data. We initiated this project because the complex genetic basis of MDD requires large population-wide studies with uniform in-depth phenotyping. For standardized phenotyping we developed the LIDAS (LIfetime Depression Assessment Survey), which then was used to measure MDD in 11 Dutch cohorts. Data from these cohorts were combined with diagnostic interview depression data from 5 clinical cohorts to create a dataset of N = 29,650 lifetime MDD cases (22%) meeting DSM-5 criteria and 94,300 screened controls. In addition, genomewide genotype data from the cohorts were assembled into a genomewide association study (GWAS) dataset of N = 66,684 Dutch individuals (25.3% cases). Phenotype data include DSM-5-based MDD diagnoses, sociodemographic variables, information on lifestyle and BMI, characteristics of depressive symptoms and episodes, and psychiatric diagnosis and treatment history. We describe the establishment and harmonization of the BIONIC phenotype and GWAS datasets and provide an overview of the available information and sample characteristics. Our next step is the GWAS of lifetime MDD in the Netherlands, with future plans including fine-grained genetic analyses of depression characteristics, international collaborations and multi-omics studies.

Modification of Huntington's disease by short tandem repeats.

Expansions of glutamine-coding CAG trinucleotide repeats cause a number of neurodegenerative diseases, including Huntington's disease and several of spinocerebellar ataxias. In general, age-at-onset of the polyglutamine diseases is inversely correlated with the size of the respective inherited expanded CAG repeat. Expanded CAG repeats are also somatically unstable in certain tissues, and age-at-onset of Huntington's disease corrected for individual HTT CAG repeat length (i.e. residual age-at-onset), is modified by repeat instability-related DNA maintenance/repair genes as demonstrated by recent genome-wide association studies. Modification of one polyglutamine disease (e.g. Huntington's disease) by the repeat length of another (e.g. ATXN3, CAG expansions in which cause spinocerebellar ataxia 3) has also been hypothesized. Consequently, we determined whether age-at-onset in Huntington's disease is modified by the CAG repeats of other polyglutamine disease genes. We found that the CAG measured repeat sizes of other polyglutamine disease genes that were polymorphic in Huntington's disease participants but did not influence Huntington's disease age-at-onset. Additional analysis focusing specifically on ATXN3 in a larger sample set (n = 1388) confirmed the lack of association between Huntington's disease residual age-at-onset and ATXN3 CAG repeat length. Additionally, neither our Huntington's disease onset modifier genome-wide association studies single nucleotide polymorphism data nor imputed short tandem repeat data supported the involvement of other polyglutamine disease genes in modifying Huntington's disease. By contrast, our genome-wide association studies based on imputed short tandem repeats revealed significant modification signals for other genomic regions. Together, our short tandem repeat genome-wide association studies show that modification of Huntington's disease is associated with short tandem repeats that do not involve other polyglutamine disease-causing genes, refining the landscape of Huntington's disease modification and highlighting the importance of rigorous data analysis, especially in genetic studies testing candidate modifiers.

Introgression and incomplete lineage sorting blurred phylogenetic relationships across the genomes of sclerophyllous oaks from southwest China.

Resolving evolutionary relationships among closely related species with interspecific gene flow is challenging. Genome-scale data provide opportunities to clarify complex evolutionary relationships in closely related species and to observe variations in species relationships across the genomes of such species. The Himalayan-Hengduan subalpine oaks have a nearly completely sympatric distribution in southwest China and probably constitute a syngameon. In this study, we mapped resequencing data from different species in this group to the Quercus aquifolioides reference genome to obtain a high-quality filtered single nucleotide polymorphism (SNP) dataset. We also assembled their plastomes. We reconstructed their phylogenetic relationships, explored the level and pattern of introgression among these species and investigated gene tree variation in the genomes of these species using sliding windows. The same or closely related plastomes were found to be shared extensively among different species within a specific geographical area. Phylogenomic analyses of genome-wide SNP data found that most oaks in the Himalayan-Hengduan subalpine clade showed genetic coherence, but several species were found to be connected by introgression. The gene trees obtained using sliding windows showed that the phylogenetic relationships in the genomes of oaks are highly heterogeneous and therefore highly obscured. Our study found that all the oaks of the Himalayan-Hengduan subalpine clade from southwest China form a syngameon. The obscured phylogenetic relationships observed empirically across the genome are best explained by interspecific gene flow in conjunction with incomplete lineage sorting.

Contrasting patterns of genetic structure and population demography in two Dicraeus species feeding on bamboo flowers in Japan

Abstract Population outbreaks of insects provide an excellent opportunity to elucidate their evolutionary history. In this study, we uncovered the phylogenetic relationships of Dicraeus spp. in Japan and compared the population genetic structure of two bamboo-feeding Dicraeus spp., Dicraeus phyllostachyus and Dicraeus nartshukae, by leveraging a recent nationwide mass flowering event of several bamboo species in Japan. We performed analysis using mitochondrial COI gene sequencing and genome-wide single-nucleotide polymorphisms detected by multiplexed inter-simple sequence repeat genotyping by sequencing (MIG-seq). We found genetic similarity in COI haplotypes over a large geographical range in the two Dicraeus spp., but, for both species, no concordance was detected between the collected host plants and genetic structures. In the genome-wide single-nucleotide polymorphism data, no clear genetic clusters were detected for D. phyllostachyus populations, whereas D. nartshukae populations were divided into two groups, in the north and south of Japan. Demographic analysis indicated that D. phyllostachyus had greatly expanded its population recently compared with D. nartshukae. These contrasting patterns between the two species might reflect the differences in the demographic history between the species. Our study suggested that the mass flowering event might cascade to the genetic structure of the florivorous insects through trophic interactions.

The demographic history and adaptation of Canarian goat breeds to environmental conditions through the use of genome-wide SNP data

BackgroundThe presence of goats in the Canary Islands dates back to the late 1st millennium BC, which coincides with the colonization by the Amazigh settlers. However, the exact geographic origin of Canarian goats is uncertain since the Amazigh peoples were distributed over a wide spatial range. Nowadays, three Canarian breeds (Palmera, Majorera and Tinerfeña) are officially recognized, along with two distinct South and North Tinerfeña ecotypes, with the South Tinerfeña and Majorera goats thriving in arid and dry semi-desertic environments and the Palmera and North Tinerfeña goats are adapted to humid and temperate areas that are influenced by trade winds. Genotypes for 224 Canarian goats were generated using the Illumina Goat single nucleotide polymorphism (SNP)50 BeadChip. By merging these data with the genotypes from 1007 individuals of African and Southern European ancestry, our aim was to ascertain the geographic origin of the Canarian goats and identify genes associated with adaptation to diverse environmental conditions.ResultsThe diversity indices of the Canarian breeds align with most of those of the analyzed local breeds from Africa and Europe, except for the Palmera goats that showed lower levels of genetic variation. The Canarian breeds demonstrate a significant genetic differentiation compared to other populations, which indicates a history of prolonged geographic isolation. Moreover, the phylogenetic reconstruction indicated that the ancestry of the Canarian goats is fundamentally North African rather than West African. The ADMIXTURE and the TreeMix analyses showed no evidence of gene flow between Canarian goats and other continental breeds. The analysis of runs of homozygosity (ROH) identified 13 ROH islands while the window-based FST method detected 25 genomic regions under selection. Major signals of selection were found on Capra hircus (CHI) chromosomes 6, 7, and 10 using various comparisons and methods.ConclusionsThis genome-wide analysis sheds new light on the evolutionary history of the four breeds that inhabit the Canary Islands. Our findings suggest a North African origin of the Canarian goats. In addition, within the genomic regions highlighted by the ROH and FST approaches, several genes related to body size and heat tolerance were identified.

Genome-wide SNP analysis provides insights into the XX/XY sex-determination system in silver barb (Barbonymus gonionotus).

Silver barb (Barbonymus gonionotus) is among the most economically important freshwater fish species in Thailand. It ranks fourth in economic value and third in production weight for fisheries and culture in Thailand. An XX/XY sex-determination system based on gynogenesis was previously reported for this fish. In this study, the molecular basis underlying the sex-determination system was further investigated. Genome-wide single-nucleotide polymorphism data were generated for 32 captive-bred silver barb individuals, previously scored by phenotypic sex, to identify sex-linked regions associated with sex determination. Sixty-three male-linked loci, indicating putative XY chromosomes, were identified. Male-specific loci were not observed, which indicates that the putative Y chromosome is young and the sex determination region is cryptic. A homology search revealed that most male-linked loci were homologous to the Mariner/Tc1 and Gypsy transposable elements and are probably the remnants of an initial accumulation of repeats on the Y chromosome from the early stages of sex chromosome differentiation. This research provides convincing insights into the mechanism of sex determination and reveals the potential sex determination regions in silver barb. The study provides the basic data necessary for increasing the commercial value of silver barbs through genetic improvements.

Evidence for selfing in a vertebrate from whole-genome sequencing.

A growing number of recent genomic studies report asexual parthenogenetic reproduction in a wide range of taxa, including vertebrate species from the reptile, bird, and fish lineages. Yet, self-fertilization (selfing) has been recorded only in a single vertebrate, the mangrove killifish Kryptolebias marmoratus In cichlid fishes, sex determination is notably diverse and can be influenced by the environment, and sequential hermaphroditism has been reported for some species. Here, we present evidence for a case of facultative selfing in the cichlid fish Benitochromis nigrodorsalis, which is otherwise known as biparentally reproducing ovophilic mouthbrooder from Western Africa. Our laboratory observations revealed that a wild-caught individual produced repeatedly viable offspring in absence of a mating partner. By analyzing genome-wide single-nucleotide polymorphism (SNP) data, we compare that individual and two of its offspring to shed light on its reproductive mode. First, our results confirm uniparental reproduction. Second, overall heterozygosity is reduced in the offspring compared with outbred individuals. Retained maternal heterozygosity in the offspring is ∼51%, which is close to the theoretically expected value of a heterozygosity reduction of 50% by selfing. Heterozygosity patterns along individual chromosomes do not point to alternative parthenogenetic reproductive mechanisms like automixis by terminal or central fusion. Facultative selfing may represent an adaptive strategy ensuring reproduction when mating partners are absent and, hence, contribute to the cichlids' enormous evolutionary success.

Revisiting a cryptic species complex: interspecific gene flow among woodland salamanders in the Blue Ridge Mountains of northern Georgia, USA

Abstract The complex topography of mountainous landscapes can create biogeographical barriers and promote allopatric speciation—even among morphologically cryptic species. However, these biogeographical barriers may not be permanent, and secondary contact between previously independently evolving populations can result in hybridization and introgression. In this study, we used genome-wide single nucleotide polymorphism data to re-examine a classic example of cryptic species in the Appalachian Mountains: the slimy salamander (Plethodon glutinosus) species complex. This group of morphologically similar terrestrial salamanders includes 14 species with parapatric distributions across the eastern USA. We focused on the Chattahoochee slimy salamander (P. chattahoochee) in northern Georgia and used a series of complementary phylogenomic and population genomic analyses to evaluate spatial genetic structure within this species and admixture with at least four other species of parapatric woodland salamanders. Our results highlight the utility of genomic data in clarifying species boundaries and underline the difficulty of species delimitation in organisms with complex evolutionary histories.

Discovering SNP-disease relationships in genome-wide SNP data using an improved harmony search based on SNP locus and genetic inheritance patterns.

Advances in high-throughput sequencing technologies have made it possible to access millions of measurements from thousands of people. Single nucleotide polymorphisms (SNPs), the most common type of mutation in the human genome, have been shown to play a significant role in the development of complex and multifactorial diseases. However, studying the synergistic interactions between different SNPs in explaining multifactorial diseases is challenging due to the high dimensionality of the data and methodological complexities. Existing solutions often use a multi-objective approach based on metaheuristic optimization algorithms such as harmony search. However, previous studies have shown that using a multi-objective approach is not sufficient to address complex disease models with no or low marginal effect. In this research, we introduce a locus-driven harmony search (LDHS), an improved harmony search algorithm that focuses on using SNP locus information and genetic inheritance patterns to initialize harmony memories. The proposed method integrates biological knowledge to improve harmony memory initialization by adding SNP combinations that are likely candidates for interaction and disease causation. Using a SNP grouping process, LDHS generates harmonies that include SNPs with a higher potential for interaction, resulting in greater power in detecting disease-causing SNP combinations. The performance of the proposed algorithm was evaluated on 200 synthesized datasets for disease models with and without marginal effect. The results show significant improvement in the power of the algorithm to find disease-related SNP sets while decreasing computational cost compared to state-of-the-art algorithms. The proposed algorithm also demonstrated notable performance on real breast cancer data, showing that integrating prior knowledge can significantly improve the process of detecting disease-related SNPs in both real and synthesized data.

Genetic characterisation of the Connemara pony and the Warmblood horse using a within-breed clustering approach

BackgroundThe Connemara pony (CP) is an Irish breed that has experienced varied selection by breeders over the last fifty years, with objectives ranging from the traditional hardy pony to an agile athlete. We compared these ponies with well-studied Warmblood (WB) horses, which are also selectively bred for athletic performance but with a much larger census population. Using genome-wide single nucleotide polymorphism (SNP) and whole-genome sequencing data from 116 WB (94 UK WB and 22 European WB) and 36 CP (33 UK CP and 3 US CP), we studied the genomic diversity, inbreeding and population structure of these breeds.ResultsThe k-means clustering approach divided both the CP and WB populations into four genetic groups, among which the CP genetic group 1 (C1) associated with non-registered CP, C4 with US CP, WB genetic group 1 (W1) with Holsteiners, and W3 with Anglo European and British WB. Maximum and mean linkage disequilibrium (LD) varied significantly between the two breeds (mean from 0.077 to 0.130 for CP and from 0.016 to 0.370 for WB), but the rate of LD decay was generally slower in CP than WB. The LD block size distribution peaked at 225 kb for all genetic groups, with most of the LD blocks not exceeding 1 Mb. The top 0.5% harmonic mean pairwise fixation index (FST) values identified ontology terms related to cancer risk when the four CP genetic groups were compared. The four CP genetic groups were less inbred than the WB genetic groups, but C2, C3 and C4 had a lower proportion of shorter runs of homozygosity (ROH) (74 to 76% < 4 Mb) than the four WB genetic groups (80 to 85% < 4 Mb), indicating more recent inbreeding. The CP and WB genetic groups had a similar ratio of effective number of breeders (Neb) to effective population size (Ne).ConclusionsDistinct genetic groups of individuals were revealed within each breed, and in WB these genetic groups reflected population substructure better than studbook or country of origin. Ontology terms associated with immune and inflammatory responses were identified from the signatures of selection between CP genetic groups, and while CP were less inbred than WB, the evidence pointed to a greater degree of recent inbreeding. The ratio of Neb to Ne was similar in CP and WB, indicating the influence of popular sires is similar in CP and WB.

Genome-wide single-nucleotide polymorphism data and mitochondrial hypervariable region 1 nucleotide sequence reveal the origin of the Akhal-Teke horse.

The study investigated the origin of the Akhal-Teke horse using genome-wide single-nucleotide polymorphism (SNP) data and mitochondrial hypervariable region 1 (HVR-1) nucleotide sequences. Genome-wide SNP data from 22 breeds (481 horses) and mitochondrial HVR-1 sequences from 24 breeds (544 sequences) worldwide to examine the origin of the Akhal- Teke horse. The data were analyzed using principal component analysis, linkage disequilibrium analysis, neighbor-joining dendrograms, and ancestry inference to determine the population relationships, ancestral source, genetic structure, and relationships with other varieties. A close genetic relationship between the Akhal-Teke horse and horses from the Middle East was found. Analysis of mitochondrial HVR-1 sequences showed that there were no shared haplotypes between the Akhal-Teke and Tarpan horses, and the mitochondrial data indicated that the Akhal-Teke horse has not historically expanded its group. Ancestral inference suggested that Arabian and Caspian horses were the likely ancestors of the Akhal- Teke horse. The Akhal-Teke horse originated in the Middle East.

First genome-wide association study for lymphatic filariasis in a West African population points to a human leukocyte antigen-mediated disease pathophysiology

ObjectivesLymphatic filariasis (LF) represents a parasitic disease caused by filarial nematodes. Although some infected individuals present an asymptomatic course, others suffer severe chronic lymphatic pathology, including lymphedema, hydrocele, and elephantiasis. Several studies have shown that host genetic factors influence LF susceptibility and chronic pathology. The current study aimed to conduct the first genome-wide association study to systematically determine LF susceptibility. MethodsWe analyzed genome-wide single-nucleotide polymorphism data from 1459 LF cases and 1492 asymptomatic controls of West African (Ghanaian) descent. ResultsWe identified two independent genome-wide significant associated genetic variants near the genes HLA-DQB2 (rs7742085) and HLA-DQA1 (rs4959107) contributing to LF and/or lymphedema susceptibility (P <5.0 × 10−8, odds ratios [ORs] >1.30). We also observed suggestive evidence of LF associations (P <1.0 × 10−6) at two non-HLA loci, near the genes ZFHX4-AS1 (rs79562145) and CHP2 (rs12933387). In contrast, we could not replicate any previously reported LF associations drawn from candidate gene association studies. On the polygenic level, we show that our genome-wide association study data explain 24-42% of LF heritability, depending on an assumed population prevalence of 0.5-5.0%. ConclusionOur findings point to an involvement of HLA-mediated immune mechanisms in LF pathophysiology.

Population genetic structure of the maize weevil, Sitophilus zeamais, in southern Mexico.

The maize weevil, Sitophilus zeamais, is a ubiquitous pest of maize and other cereal crops worldwide and remains a threat to food security in subsistence communities. Few population genetic studies have been conducted on the maize weevil, but those that exist have shown that there is very little genetic differentiation between geographically dispersed populations and that it is likely the species has experienced a recent range expansion within the last few hundred years. While the previous studies found little genetic structure, they relied primarily on mitochondrial and nuclear microsatellite markers for their analyses. It is possible that more fine-scaled population genetic structure exists due to local adaptation, the biological limits of natural species dispersal, and the isolated nature of subsistence farming communities. In contrast to previous studies, here, we utilized genome-wide single nucleotide polymorphism data to evaluate the genetic population structure of the maize weevil from the southern and coastal Mexican states of Oaxaca and Chiapas. We employed strict SNP filtering to manage large next generation sequencing lane effects and this study is the first to find fine-scale genetic population structure in the maize weevil. Here, we show that although there continues to be gene flow between populations of maize weevil, that fine-scale genetic structure exists. It is possible that this structure is shaped by local adaptation of the insects, the movement and trade of maize by humans in the region, geographic barriers to gene flow, or a combination of these factors.