Identification of Somatostatin Receptor Subtype 1 (SSTR1) Gene Polymorphism and Their Association with Growth Traits in Hulun Buir Sheep

This study was conducted to evaluate SSTR1 gene polymorphisms and their association with growth traits in Hulun Buir sheep. We followed 233 Hulun Buir sheep from birth to 16 months of age, born in the same pasture and on the same year under a consistent grazing conditions. The body weight (BW), body height (BH), body length (BL), chest circumference (ChC), chest depth (ChD), chest width (ChW), hip width (HW), and cannon circumference (CaC) were measured and recorded at birth, 4 months, 9 months, and 16 months of age. The polymorphisms of the SSTR1 gene in Hulun Buir sheep were excavated using exon sequencing, and association analyses of between SNPs and growth traits at each growth stage were conducted. The results showed that there were four SNPs in Exon 2 of the SSTR1 gene, SNP1, SNP2, and SNP3 were low mutation sites, and SNP4 was a moderate mutation site. Four SNPs were consistent with Hardy–Weinberg equilibrium, and all of them were synonymous mutations. The association analyses found that the genotypes of SNP2 were significantly associated with WW and BH at 4 months of age, BW, BL, ChC, and HW at 9 months of age (p < 0.05), and extremely significantly associated with ChD at 4 and 9 months of age (p < 0.01). There were significant associations between SNP3 and BH at 9 months of age, between SNP4 and ChD, ChW, and CaC at 9 months of age, and BW and ChC at 16 months of age (p < 0.05). There were no detectable associations with growth traits among the seven haplotypes between the SNP1, 3, and 4 of a strong linkage disequilibrium (p > 0.05). These results indicated that SNP2, SNP3, and SNP4 may be used as molecular markers for growth traits of Hulun Buir sheep.

Meta-Analysis of Transcriptome-Wide Association Studies across 13 Brain Tissues Identified Novel Clusters of Genes Associated with Nicotine Addiction

Genome-wide association studies (GWAS) have identified and reproduced thousands of diseases associated loci, but many of them are not directly interpretable due to the strong linkage disequilibrium among variants. Transcriptome-wide association studies (TWAS) incorporated expression quantitative trait loci (eQTL) cohorts as a reference panel to detect associations with the phenotype at the gene level and have been gaining popularity in recent years. For nicotine addiction, several important susceptible genetic variants were identified by GWAS, but TWAS that detected genes associated with nicotine addiction and unveiled the underlying molecular mechanism were still lacking. In this study, we used eQTL data from the Genotype-Tissue Expression (GTEx) consortium as a reference panel to conduct tissue-specific TWAS on cigarettes per day (CPD) over thirteen brain tissues in two large cohorts: UK Biobank (UKBB; number of participants (N) = 142,202) and the GWAS & Sequencing Consortium of Alcohol and Nicotine use (GSCAN; N = 143,210), then meta-analyzing the results across tissues while considering the heterogeneity across tissues. We identified three major clusters of genes with different meta-patterns across tissues consistent in both cohorts, including homogenous genes associated with CPD in all brain tissues; partially homogeneous genes associated with CPD in cortex, cerebellum, and hippocampus tissues; and, lastly, the tissue-specific genes associated with CPD in only a few specific brain tissues. Downstream enrichment analyses on each gene cluster identified unique biological pathways associated with CPD and provided important biological insights into the regulatory mechanism of nicotine dependence in the brain.

An HIV Vaccine Protective Allele in FCGR2C Associates With Increased Odds of Perinatal HIV Acquisition

In the Thai RV144 HIV-1 vaccine trial, a three-variant haplotype within the Fc gamma receptor 2C gene (FCGR2C) reduced the risk of HIV-1 acquisition. A follow-on trial, HVTN702, of a similar vaccine candidate found no efficacy in South Africa, where the predominant population is polymorphic for only a single variant in the haplotype, c.134-96C&gt;T (rs114945036). To investigate a role for this variant in HIV-1 acquisition in South Africans, we used the model of maternal-infant HIV-1 transmission. A nested case-control study was conducted of infants born to mothers living with HIV-1, comparing children with perinatally-acquired HIV-1 (cases, n = 176) to HIV-1-exposed uninfected children (controls, n = 349). All had received nevirapine for prevention of mother-to-child transmission. The FCGR2C copy number and expression variants (c.−386G&gt;C, c.−120A&gt;T c.169T&gt;C, and c.798+1A&gt;G) were determined using a multiplex ligation-dependent probe amplification assay and the c.134-96C&gt;T genotype with Sanger sequencing. The copy number, genotype and allele carriage were compared between groups using univariate and multivariate logistic regression. The FCGR2C c.134-96C&gt;T genotype distribution and copy number differed significantly between HIV-1 cases and exposed-uninfected controls (P = 0.002, PBonf = 0.032 and P = 0.010, PBonf = &gt; 0.05, respectively). The FCGR2C c.134-96T allele was overrepresented in the cases compared to the controls (58% vs 42%; P = 0.001, PBonf = 0.016). Adjusting for birthweight and FCGR2C copy number, perinatal HIV-1 acquisition was associated with the c.134-96C&gt;T (AOR = 1.89; 95% CI 1.25-2.87; P = 0.003, PBonf = 0.048) and c.169C&gt;T (AOR = 2.39; 95% CI 1.45-3.95; P = 0.001, PBonf = 0.016) minor alleles but not the promoter variant at position c.−386G&gt;C. The c.134-96C&gt;T variant was in strong linkage disequilibrium with the c.169C&gt;T variant, but remained significantly associated with perinatal acquisition when adjusted for c.169C&gt;T in multivariate analysis. In contrast to the protective effect observed in the Thai RV144 trial, we found the FCGR2C variant c.134-96T-allele associated with increased odds of perinatal HIV-1 acquisition in South African children. These findings, taken together with a similar deleterious association found with HIV-1 disease progression in South African adults, highlight the importance of elucidating the functional relevance of this variant in different populations and vaccination/disease contexts.

Meta-Analysis of Transcriptome-Wide Association Studies Across 13 Brain Tissues Identified Novel Clusters of Genes Associated with Nicotine Addiction

Genome-wide association studies (GWAS) have identified and reproduced thousands of diseases associated loci but many of them are not directly interpretable due to the strong linkage disequilibrium among variants. Transcriptome-wide association studies (TWAS) incorporated expression quantitative trait loci (eQTL) cohorts as reference panel to detect associations with the phenotype at the gene level and were gaining popularity in recent years. For nicotine addiction, several important susceptible genetic variants were identified by GWAS, but TWAS that detected genes associated with nicotine addiction and unveiled the underlying molecular mechanism were still lacking. In this study, we used eQTL data from the Genotype-Tissue Expression (GTEx) consortium as reference panel to conduct tissue specific TWAS on cigarettes per day (CPD) over 13 brain tissues in two large cohorts: UK Biobank (UKBB; N=142,202) and the GWAS &amp;amp; Sequencing Consortium of Alcohol and Nicotine use (GSCAN; N=143,210), and then meta-analyzed the results across tissues while considering the heterogeneity across tissues. We identified three major clusters of genes with different meta-patterns across tissues consistent in both cohorts, including homogenous genes associated with CPD in all brain tissues, partially homogeneous genes associated with CPD in cortex, cerebellum and hippocampus tissues, and lastly the tissue-specific genes associated with CPD in only few specific brain tissues. Downstream enrichment analyses on each gene cluster identified unique biological pathways associated with CPD and provided important biological insights into the regulatory mechanism of nicotine dependence in the brain.

Sequence analysis of tumor necrosis factor-α (TNF-α) promoter, 5'UTR and exon1 and association of rs361525 (-238 G>A) with BMI

Obesity is among the most common complex diseases with a high rate of morbidity and mortality globally and locally in Kuwait. Tumor Necrosis Factor-α(TNF-α) is a pro-inflammatory cytokine that is primarily secreted by monocytes/macrophage. Increased expression of TNF-α has been observed in the adipodse tissue of obese subjects that could disrupt lipid metabolism and lead to and lead to sustained obese state and obesity-related diseases. The human TNF-α promoter exhibits a high number of genetic variants, mainly single nucleotide polymorphisms (SNPs) that have been shown to influence the level of transcription in association with diseases. The human TNF-α genetic variants have never been fully reported in Arabs, therefore, we aimed to identify these variants by sequencing the TNF-α promoter, 5' UTR, and exon 1 in 290 Kuwaiti Arabs. As a result, we identified 14 genetic variants, including one novel SNP. Two promoter SNPs; rs1800750 (- 376G>A) and rs361525 (-238G>A) were found to be in strong linkage disequilibrium (r2 = 0.73) and (D′ = 1, LOD = 32). To investigate the association of rs361525 (-238G>A) with obesity, we genotyped an additional 573 samples of the general Kuwaiti population by Real-time PCR (total n=863). Linear and logistic regression analysis have not shown any significant association in carriers of the A allele of rs361525 with continuous and categorical BMI, respectively. This is the first study in the Middle East and Kuwait that has sequenced and identified the common, rare and novel genetic variants of TNF-α promoter, 5'UTR and exon 1 in Arabs.

Broad-scale puma connectivity could restore genomic diversity to fine-scale coastal populations

Urbanization is decreasing wildlife habitat and connectivity worldwide, including for apex predators, such as the puma (Puma concolor). Puma populations along California's central and southern coastal habitats have experienced rapid fragmentation from development, leading to calls for demographic and genetic management. To address urgent conservation genomic concerns, we used double-digest restriction-site associated DNA (ddRAD) sequencing to analyze 16,285 genome-wide single-nucleotide polymorphisms (SNPs) from 401 broadly sampled pumas. Our analyses indicated support for 4–10 geographically nested, broad- to fine-scale genetic clusters. At the broadest scale, the 4 genetic clusters had high genetic diversity and exhibited low linkage disequilibrium, indicating pumas have retained statewide genomic diversity. However, multiple lines of evidence indicated substructure, including 10 fine-scale genetic clusters, some of which exhibited allelic fixation and linkage disequilibrium. Fragmented populations along the Southern Coast and Central Coast had particularly low genetic diversity and strong linkage disequilibrium, indicating genetic drift and close inbreeding. Our results demonstrate that genetically at-risk populations are typically nested within a broader-scale group of interconnected populations that collectively retains high genetic diversity and heterogeneous fixations. Thus, extant variation at the broader scale has potential to restore diversity to local populations if management actions can enhance vital gene flow and recombine locally sequestered genetic diversity. These state- and genome-wide results are critically important for science-based conservation and management practices. Our broad- and fine-scale population genomic analysis highlights the information that can be gained from population genomic studies aiming to provide guidance for fragmented population conservation management.

Five SNPs Within the FGF5 Gene Significantly Affect Both Wool Traits and Growth Performance in Fine-Wool Sheep (Ovis aries)

Fibroblast growth factor 5 (FGF5) gene, a member of fibroblast growth factor superfamily, plays significant roles in the regulation of the hair growth cycle during the development of mammalian hair follicles as well as the skeletal muscle development. In this study, DNA sequencing was used to scan the putative SNPs within the full-length of FGF5 gene, and SNPscan high-throughput technique was applied in the individual genotyping of 604 crossbred sheep. 10 SNPs were identified within FGF5 gene while five of them located in intron 1 could be genotyped, namely SNP1 (g. 105914953 G &gt; A), SNP2 (g. 105922232 T &gt; C), SNP3 (g. 105922244 A &gt; G), SNP4 (g. 105922334 A &gt; T) and SNP5 (g. 105922340 G &gt; T). All these SNPs were in accord with the Hardy-Weinberg equilibrium (P &gt; 0.05), and displayed the moderate polymorphism with PIC values ranging from 0.302 to 0.374. Thereafter, the correlation analysis between each SNP locus and economic traits including wool length, greasy wool weight and growth performance of sheep was systematically implemented. In our results, SNP1, SNP3, SNP4 and SNP5 were significantly associated with wool length, greasy wool weight and growth traits of SG sheep (P &lt; 0.05); SNP1, SNP2, SNP3, and SNP4 were significantly correlated with wool length and growth traits of SSG sheep (P &lt; 0.05). Meanwhile, our study revealed a strong linkage disequilibrium (LD) relationship among these SNPs (r2 &gt; 0.33), except for SNP3 and SNP4 sites (r2 = 0.30). Combination genotype analysis showed that combination genotypes were significantly associated with mean fiber diameter of SG (P &lt; 0.05), and body weight trait of SSG (P &lt; 0.01). The above findings suggested that these SNP loci might affect economic traits synergistically and could be regarded as potential molecular markers for improving both wool production and growth performance of fine-wool sheep, which lay a molecular foundation for the breeding of fine dual-purpose sheep thereby accelerating the pace of sheep breeding.

Assessment of linkage disequilibrium patterns between structural variants and single nucleotide polymorphisms in three commercial chicken populations

BackgroundStructural variants (SV) are causative for some prominent phenotypic traits of livestock as different comb types in chickens or color patterns in pigs. Their effects on production traits are also increasingly studied. Nevertheless, accurately calling SV remains challenging. It is therefore of interest, whether close-by single nucleotide polymorphisms (SNPs) are in strong linkage disequilibrium (LD) with SVs and can serve as markers. Literature comes to different conclusions on whether SVs are in LD to SNPs on the same level as SNPs to other SNPs. The present study aimed to generate a precise SV callset from whole-genome short-read sequencing (WGS) data for three commercial chicken populations and to evaluate LD patterns between the called SV and surrounding SNPs.ResultsThe final callset consisted of 12,294,329 bivariate SNPs, 4,301 deletions (DEL), 224 duplications (DUP), 218 inversions (INV) and 117 translocation breakpoints (BND). While average LD between DELs and SNPs was at the same level as between SNPs and SNPs, LD between other SVs and SNPs was strongly reduced (DUP: 40 %, INV: 27 %, BND: 19 % of between-SNP LD). A main factor for the reduced LD was the presence of local minor allele frequency differences, which accounted for 50 % of the difference between SNP – SNP and DUP – SNP LD. This was potentially accompanied by lower genotyping accuracies for DUP, INV and BND compared with SNPs and DELs. An evaluation of the presence of tag SNPs (SNP in highest LD to the variant of interest) further revealed DELs to be slightly less tagged by WGS SNPs than WGS SNPs by other SNPs. This difference, however, was no longer present when reducing the pool of potential tag SNPs to SNPs located on four different chicken genotyping arrays.ConclusionsThe results imply that genomic variance due to DELs in the chicken populations studied can be captured by different SNP marker sets as good as variance from WGS SNPs, whereas separate SV calling might be advisable for DUP, INV, and BND effects.

Whole Genome Sequencing Reveals Multiple Linked Genetic Variants on Canine Chromosome 12 Associated with Risk for Symmetrical Lupoid Onychodystrophy (SLO) in the Bearded Collie

In dogs, symmetrical lupoid onychodystrophy (SLO) results in nail loss and an abnormal regrowth of the claws. In Bearded Collies, an autoimmune nature has been suggested because certain dog leukocyte antigen (DLA) class II haplotypes are associated with the condition. A genome-wide association study of the Bearded Collie revealed two regions of association that conferred risk for disease: one on canine chromosome (CFA) 12 that encompasses the DLA genes, and one on CFA17. Case-control association was employed on whole genome sequencing data to uncover putative causative variants in SLO within the CFA12 and CFA17 associated regions. Genotype imputation was then employed to refine variants of interest. Although no SLO-associated protein-coding variants were identified on CFA17, multiple variants, many with predicted damaging effects, were identified within potential candidate genes on CFA12. Furthermore, many potentially damaging alleles were fully correlated with the presence of DLA class II risk haplotypes for SLO, suggesting that the variants may reflect DLA class II haplotype association with disease or vice versa. Strong linkage disequilibrium in the region precluded the ability to isolate and assess the individual or combined effect of variants on disease development. Nonetheless, all were predictive of risk for SLO and, with judicious assessment, their application in selective breeding may prove useful to reduce the incidence of SLO in the breed.

Mechanisms of ecological divergence with gene flow in a reef-building coral on an isolated atoll in Western Australia

Understanding the mechanisms driving phenotypic variation in traits facing intensified selection from climate change is a crucial step in developing effective conservation and restoration initiatives. This is particularly true for reef-building corals, which are among the most vulnerable to climate change and are in dramatic decline globally. At the Rowley Shoals in Western Australia, the prominent reef flat becomes exposed on low tide and the stagnant water in the shallow atoll lagoons heats up, creating a natural laboratory for characterising the mechanisms that control phenotypic responses to different environments. We combined whole genome re-sequencing, common garden heat stress experiments, transcriptome-wide gene expression analyses, and symbiont metabarcoding to explore the mechanisms that facilitate survival in contrasting habitat conditions. Our data show that, despite high gene flow between habitats, spatially varying selection drives subtle shifts in allele frequencies at hundreds of loci. These changes were concentrated into several islands of divergence spanning hundreds of SNPs that showed strong linkage disequilibrium and were associated with a coordinated increase in minor allele frequencies in corals taken from the lagoon habitat, where the range of environmental conditions is greatest. Common garden heat stress assays showed individuals from the lagoon exhibited higher bleaching resistance than colonies from the reef slope, and RNAseq identified pronounced physiological differences between the corals from the two habitats, primarily associated with molecular pathways including cell signalling, ion transport and metabolism. Despite the pronounced physioloigical and environmental differences between habitats, metabarcoding of the Symbiodiniaceae ITS2 region revealed all colonies to be associated exclusively with the genus Cladocopium, with no detectable differences between habitats. This study contributes to the growing number of studies documenting the complex mechanisms that facilitate coral survival in extreme environments, and showcases the utility of combining multiple sequencing techniques to unravel complex climate-related traits.