Associations between genetic loci, environment factors and mental disorders: a genome-wide survival analysis using the UK Biobank data

It is well-accepted that both environment and genetic factors contribute to the development of mental disorders (MD). However, few genetic studies used time-to-event data analysis to identify the susceptibility genetic variants associated with MD and explore the role of environment factors in these associations. In order to detect novel genetic loci associated with MD based on the time-to-event data and identify the role of environmental factors in them, this study recruited 376,806 participants from the UK Biobank cohort. The MD outcomes (including overall MD status, anxiety, depression and substance use disorders (SUD)) were defined based on in-patient hospital, self-reported and death registry data collected in the UK Biobank. SPACOX approach was used to identify the susceptibility loci for MD using the time-to-event data of the UK Biobank cohort. And then we estimated the associations between identified candidate loci, fourteen environment factors and MD through a phenome-wide association study and mediation analysis. SPACOX identified multiple candidate loci for overall MD status, depression and SUD, such as rs139813674 (P value = 8.39 × 10–9, ZNF684) for overall MD status, rs7231178 (DCC, P value = 2.11 × 10–9) for depression, and rs10228494 (FOXP2, P value = 6.58 × 10–10) for SUD. Multiple environment factors could influence the associations between identified loci and MD, such as confide in others and felt hated. Our study identified novel candidate loci for MD, highlighting the strength of time-to-event data based genetic association studies. We also observed that multiple environment factors could influence the association between susceptibility loci and MD.

A replication study separates polymorphisms behind migraine with and without depression

The largest migraine genome-wide association study identified 38 candidate loci. In this study we assessed whether these results replicate on a gene level in our European cohort and whether effects are altered by lifetime depression. We tested SNPs of the loci and their vicinity with or without interaction with depression in regression models. Advanced analysis methods such as Bayesian relevance analysis and a neural network based classifier were used to confirm findings. Main effects were found for rs2455107 of PRDM16 (OR = 1.304, p = 0.007) and five intergenic polymorphisms in 1p31.1 region: two of them showed risk effect (OR = 1.277, p = 0.003 for both rs11209657 and rs6686879), while the other three variants were protective factors (OR = 0.4956, p = 0.006 for both rs12090642 and rs72948266; OR = 0.4756, p = 0.005 for rs77864828). Additionally, 26 polymorphisms within ADGRL2, 2 in REST, 1 in HPSE2 and 33 mostly intergenic SNPs from 1p31.1 showed interaction effects. Among clumped results representing these significant regions, only rs11163394 of ADGRL2 showed a protective effect (OR = 0.607, p = 0.002), all other variants were risk factors (rs1043215 of REST with the strongest effect: OR = 6.596, p = 0.003). Bayesian relevance analysis confirmed the relevance of intergenic rs6660757 and rs12128399 (p31.1), rs1043215 (REST), rs1889974 (HPSE2) and rs11163394 (ADGRL2) from depression interaction results, and the moderate relevance of rs77864828 and rs2455107 of PRDM16 from main effect analysis. Both main and interaction effect SNPs could enhance predictive power with the neural network based classifier. In summary, we replicated p31.1, PRDM16, REST, HPSE2 and ADGRL2 genes with classic genetic and advanced analysis methods. While the p31.1 region and PRDM16 are worthy of further investigations in migraine in general, REST, HPSE2 and ADGRL2 may be prime candidates behind migraine pathophysiology in patients with comorbid depression.

Genetic Diversity and Population Structure for Resistance and Susceptibility to Mastitis in Braunvieh Cattle

Mastitis is a disease that causes significant economic losses, since resistance to mastitis is a difficult trait to be improved due to its multifactorial occurrence. Therefore, our objective was to characterize a Mexican Braunvieh cattle population for genetic resistance and susceptibility to mastitis. We used 66 SNP markers for 45 candidate genes in 150 animals. The average heterozygosity was 0.445 ± 0.076, a value higher than those reported for some European breeds. The inbreeding coefficient was slightly negative for resistance to subclinical (−0.058 ± 0.055) and clinical (−0.034 ± 0.076) mastitis, possibly due to low selection for the immunological candidate genes that influence these traits. The genotypic profiles for the candidate loci per K-means group were obtained, as well as the group distribution through the graphics of the principal component analysis. The genotypic profiles showed high genetic diversity among groups. Resistance to clinical mastitis had the lowest presence of the heterozygous genotypes. Although the percentage of highly inbred animals (>50%) is up to 13.3%, there are highly heterozygous groups in terms of the studied traits, a favorable indicator of the presence of genetic diversity. The results of this study constitute evidence of the genetic potential of the Mexican Braunvieh population to improve mastitis-related traits.

Barcoded reciprocal hemizygosity analysis via sequencing illuminates the complex genetic basis of yeast thermotolerance

Decades of successes in statistical genetics have revealed the molecular underpinnings of traits as they vary across individuals of a given species. But standard methods in the field can’t be applied to divergences between reproductively isolated taxa. Genome-wide reciprocal hemizygosity mapping (RH-seq), a mutagenesis screen in an inter-species hybrid background, holds promise as a method to accelerate the progress of interspecies genetics research. Here we describe an improvement to RH-seq in which mutants harbor barcodes for cheap and straightforward sequencing after selection in a condition of interest. As a proof of concept for the new tool, we carried out genetic dissection of the difference in thermotolerance between two reproductively isolated budding yeast species. Experimental screening identified dozens of candidate loci at which variation between the species contributed to the thermotolerance trait. Hits were enriched for mitosis genes and other housekeeping factors, and among them were multiple loci with robust sequence signatures of positive selection. Together, these results shed new light on the mechanisms by which evolution solved the problems of cell survival and division at high temperature in the yeast clade, and they illustrate the power of the barcoded RH-seq approach.

Identification, Analysis, and Confirmation of Seed Storability-Related Loci in Dongxiang Wild Rice (Oryza rufipogon Griff.)

Dongxiang wild rice (Oryza rufipogon Griff.) (DXWR) has strong seed storability and identifying its elite gene resources may facilitate genetic improvements in rice seed storability. In this study, we developed two backcross inbred lines (BILs) populations, with DXWR as a common donor parent and two rice varieties (F6 and R974) as recipient parents. Bulked segregant analysis via whole genome sequencing (BSA-seq) was used to identify seed storability-related loci in the DXWR and F6 population. Two main genomic regions containing 18,550,000–20,870,000 bp on chromosome 4 and 7,860,000–9,780,000 bp on chromosome 9 were identified as candidate loci of DXWR seed storability; these overlapped partially with seed storability-related quantitative trait loci (QTLs) discovered in previous studies, suggesting that these loci may provide important regions for isolating the responsible genes. In total, 448 annotated genes were predicted within the identified regions, of which 274 and 82 had nonsynonymous and frameshift mutations, respectively. We detected extensive metabolic activities and cellular processes during seed storability and confirmed the effects of the seed storability-related candidate loci using four BILs from DXWR and R974. These results may facilitate the cloning of DXWR seed storability-related genes, thereby elucidating rice seed storability and its improvement potential.

Identification, analysis, and confirmation of seed storability-related loci in Dongxiang wild rice (Oryza rufipogon Griff.)

Dongxiang wild rice (Oryza rufipogon Griff.) (DXWR) has strong seed storability and identifying its elite gene resources may facilitate genetic improvements in rice seed storability. In this study, we developed two backcross inbred lines (BILs) populations, with DXWR as a common donor parent and two rice varieties (F6 and R974) as recipient parents. Bulked segregant analysis via whole genome sequencing (BSAseq) was used to identify seed storability related loci in the DXWR and F6 population. Two main genomic regions containing 18,550,000 to 20,870,000 bp on chromosome 4 and 7,860,000 to 9,780,000 bp on chromosome 9 were identified as candidate loci of DXWR seed storability; these overlapped partially with seed storability related quantitative trait loci (QTLs) discovered in previous studies, suggesting that these loci may provide important regions for isolating the responsible genes. In total, 448 annotated genes were predicted within the identified regions, of which 274 and 82 had nonsynonymous and frameshift mutations, respectively. We detected extensive metabolic activities and cellular processes during seed storability and confirmed the effects of the seed storability related candidate loci using four BILs from DXWR and R974. These results may facilitate the cloning of DXWR seed storability-related genes, thereby elucidating rice seed storability and its improvement potential.

Genomic investigations provide insights into the mechanisms of resilience to heterogeneous habitats of the Indian Ocean in a pelagic fish

The adaptive genetic variation in response to heterogeneous habitats of the Indian Ocean was investigated in the Indian oil sardine using ddRAD sequencing to understand the subpopulation structure, stock complexity, mechanisms of resilience, and vulnerability in the face of climate change. Samples were collected from different ecoregions of the Indian ocean and ddRAD sequencing was carried out. Population genetic analyses revealed that samples from the Gulf of Oman significantly diverged from other Indian Ocean samples. SNP allele-environment correlation revealed the presence of candidate loci correlated with the environmental variables like annual sea surface temperature, chlorophyll-a, and dissolved oxygen concentration which might represent genomic regions allegedly diverging as a result of local adaptation. Larval dispersal modelling along the southwest coast of India indicated a high dispersal rate. The two major subpopulations (Gulf of Oman and Indian) need to be managed regionally to ensure the preservation of genetic diversity, which is crucial for climatic resilience.

Rapid Mining of Candidate Genes for Verticillium Wilt Resistance in Cotton Based on BSA-Seq Analysis

Cotton is a globally important cash crop. Verticillium wilt (VW) is commonly known as “cancer” of cotton and causes serious loss of yield and fiber quality in cotton production around the world. Here, we performed a BSA-seq analysis using an F2:3 segregation population to identify the candidate loci involved in VW resistance. Two QTLs (qvw-D05-1 and qvw-D05-2) related to VW resistance in cotton were identified using two resistant/susceptible bulks from the F2 segregation population constructed by crossing the resistant cultivar ZZM2 with the susceptible cultivar J11. A total of 30stop-lost SNPs and 42 stop-gained SNPs, which included 17 genes, were screened in the qvw-D05-2 region by SnpEff analysis. Further analysis of the transcriptome data and qRT-PCR revealed that the expression level of Ghir_D05G037630 (designated as GhDRP) varied significantly at certain time points after infection with V. dahliae. The virus-induced gene silencing of GhDRP resulted in higher susceptibility of the plants to V. dahliae than the control, suggesting that GhDRP is involved in the resistance to V. dahlia infection. This study provides a method for rapid mining of quantitative trait loci and screening of candidate genes, as well as enriches the genomic information and gene resources for the molecular breeding of disease resistance in cotton.

Genetic diversity of candidate loci linked to Mycobacterium tuberculosis resistance to bedaquiline, delamanid and pretomanid

Tuberculosis (TB), caused by Mycobacterium tuberculosis, is one of the deadliest infectious diseases worldwide. Multidrug and extensively drug-resistant strains are making disease control difficult, and exhausting treatment options. New anti-TB drugs bedaquiline (BDQ), delamanid (DLM) and pretomanid (PTM) have been approved for the treatment of multi-drug resistant TB, but there is increasing resistance to them. Nine genetic loci strongly linked to resistance have been identified (mmpR5, atpE, and pepQ for BDQ; ddn, fgd1, fbiA, fbiB, fbiC, and fbiD for DLM/PTM). Here we investigated the genetic diversity of these loci across >33,000 M. tuberculosis isolates. In addition, epistatic mutations in mmpL5-mmpS5 as well as variants in ndh, implicated for DLM/PTM resistance in M. smegmatis, were explored. Our analysis revealed 1,227 variants across the nine genes, with the majority (78%) present in isolates collected prior to the roll-out of BDQ and DLM/PTM. We identified phylogenetically-related mutations, which are unlikely to be resistance associated, but also high-impact variants such as frameshifts (e.g. in mmpR5, ddn) with likely functional effects, as well as non-synonymous mutations predominantly in MDR-/XDR-TB strains with predicted protein destabilising effects. Overall, our work provides a comprehensive mutational catalogue for BDQ and DLM/PTM associated genes, which will assist with establishing associations with phenotypic resistance; thereby, improving the understanding of the causative mechanisms of resistance for these drugs, leading to better treatment outcomes.