Effects of inbreeding on reproductive success in endangered North Atlantic right whales.

Only approximately 356 North Atlantic right whales (Eubalaena glacialis) remain. With extremely low levels of genetic diversity, limited options for mates, and variation in reproductive success across females, there is concern regarding the potential for genetic limitations of population growth from inbreeding depression. In this study, we quantified reproductive success of female North Atlantic right whales with a modified de-lifing approach using reproductive history information collected over decades of field observations. We used double-digest restriction site-associated sequencing to sequence approximately 2% of the genome of 105 female North Atlantic right whales and combined genomic inbreeding estimates with individual fecundity values to assess evidence of inbreeding depression. Inbreeding depression could not explain the variance in reproductive success of females, however we present evidence that inbreeding depression may be affecting the viability of inbred fetuses-potentially lowering the reproductive success of the species as a whole. Combined, these results allay some concerns that genetic factors are impacting species survival as genetic diversity is being retained through selection against inbred fetuses. While still far fewer calves are being born each year than expected, the small role of genetics underlying variance in female fecundity suggests that variance may be explained by external factors that can potentially be mitigated through protection measures designed to reduce serious injury and mortality from human activities.

Potato soup: analysis of cultivated potato gene bank populations reveals high diversity and little structure.

Cultivated potatoes are incredibly diverse, ranging from diploid to pentaploid and encompass four different species. They are adapted to disparate environments and conditions and carry unique alleles for resistance to pests and pathogens. Describing how diversity is partitioned within and among these populations is essential to understanding the potato genome and effectively utilizing landraces in breeding. This task is complicated by the difficulty of making comparisons across cytotypes and extensive admixture within section petota. We genotyped 730 accessions from the US Potato genebank including wild diploids and cultivated diploids and tetraploids using Genotype-by-sequencing. This data set allowed us to interrogate population structure and diversity as well as generate core subsets which will support breeders in efficiently screening genebank material for biotic and abiotic stress resistance alleles. We found that even controlling for ploidy, tetraploid material exhibited higher observed and expected heterozygosity than diploid accessions. In particular group chilotanum material was the most heterozygous and the only taxa not to exhibit any inbreeding. This may in part be because group chilotanum has a history of introgression not just from wild species, but landraces as well. All group chilotanum, exhibits introgression from group andigenum except clones from Southern South America near its origin, where the two groups are not highly differentiated. Moving north, we do not observe evidence for the same level of admixture back into group andigenum. This suggests that extensive history of admixture is a particular characteristic of chilotanum.

Genome-wide detection of runs of homozygosity and heterozygosity in Tunchang pigs.

Tunchang pigs, mainly distributed throughout Hainan Province of China, are well-known for their superior meat quality, crude feed tolerance, and adaptability to high temperatures and humidity. Runs of homozygosity (ROH) can provide valuable information about the inbreeding coefficient in individuals and selection signals that may reveal candidate genes associated with key functional traits. Runs of heterozygosity (ROHet) are commonly associated with balance selection, which can help us understand the adaptive evolutionary history of domestic animals. In this study, we investigated ROHs and ROHets in 88 Tunchang pigs. We also compared the estimates of inbreeding coefficients in individuals calculated based on four methods. In summary, we detected a total of 16 ROH islands in our study, and 100 genes were found within ROH regions. These genes were correlated with economically important traits such as reproduction (e.g., SERPIND1, HIRA), meat quality (e.g., PI4KA, TBX1), immunity (e.g., ESS2, RANBP1), adaption to heat stress (TXNRD2 and DGCR8), and crude food tolerance (TRPM6). Moreover, we discovered 18 ROHet islands harbouring genes associated with reproduction (e.g., ARHGEF12, BMPR2), immune system (e.g., BRD4, DNMT3B). These findings may help us design effective breeding and conservation strategies for this unique breed.

Prevalence of homologous recombination biomarkers in multiple tumor types: an observational study.

Aim: To determine the prevalence of deleterious mutations in BRCA1 and BRCA2 and in 13 genes involved in homologous recombination repair (HRR), the prevalence of genomic loss of heterozygosity and the allelic and hereditary status of BRCA1, BRCA2 and other HRR gene mutations in multiple solid tumor types. Patients & methods: This was a retrospective observational study of patients with an advanced/metastatic diagnosis in one of 15 solid tumor types, who were identified in a real-world clinico-genomic database. Results: Tumor tissue samples from 9457 patients were analyzed, among which 4.7% had known or suspected deleterious BRCA1/2 mutations. The prevalence (range) of mutations in HRR genes was 13.6% (2.4%-26.0%) and genomic loss of heterozygosity ≥16% was 20.6% (2.6-34.4%) across all tumor types. Conclusion: The prevalence of mutations varied significantly depending on the type of tumor.

The integrity of the human genome is maintained via multiple pathways of DNA repair, one of the most important of which is homologous recombination repair (HRR), which uses a sister chromatid as a template for high-fidelity restoration of altered DNA sequences. This study aimed to determine the prevalence of deleterious mutations, i.e., changes in the genetic code that interfere with proper cellular function, in the breast cancer genes BRCA1 and BRCA2 and in 13 other genes involved in HRR in various types of solid tumors in patients with advanced or metastatic cancer. The researchers found that 4.7% of tumor samples had BRCA1/2 mutations, 13.6% had mutations in any of the HRR genes and 20.6% had genomic loss of heterozygosity (gLOH) of at least 16% i.e., loss of sections of chromosomes affecting 16% or more of the genome. BRCA1/2 mutations were most common in ovarian cancer (13.1%), prostate cancer (9.3%), breast cancer (8.2%) and pancreatic cancer (4.9%). Prevalence for mutations in HRR genes ranges from 2.4 to 26.0% and gLOH ≥16% ranged from 2.6 to 34.4% depending on the tumor type. In conclusion, the prevalence of mutations in the BRCA1/2 genes, HRR genes and gLOH ≥16% varied widely across 15 tumor types.

Double Heterozygosity for Germline Mutations in Chinese Breast Cancer Patients.

Double pathogenic mutations occurring in an individual are considered a rare event. The introduction of a multiple-gene panel at Hong Kong Hereditary Breast Cancer Family Registry has allowed the identification of pathogenic variants in multiple genes, providing more information on clinical management and surveillance to the proband and their family members. Breast cancer patients who are double heterozygous (DH) for different hereditary breast and ovarian cancer syndrome (HBCO)-related genes were identified from a cohort of 3649 Chinese patients. Nine patients (0.25%) were observed to have germline DH mutations in ATM, BRCA1, BRCA2, BRIP1, CDH1, CHEK2, MSH6, PALB2, and TP53. Three probands were diagnosed with unilateral breast cancer, two patients were diagnosed with bilateral breast cancer, and four patients had multiple primary cancers. The median age for breast cancer diagnosis was an early age of 36 years. Chinese DH carriers did not show worse phenotypes or have a significantly downhill clinical presentation. However, seven out of nine (77.8%) of our DH carriers harbored a BRCA1 mutation, and four of them (44.4%) developed bilateral breast cancer, suggesting Chinese DH individuals may have a higher chance of having bilateral breast cancer than other populations (p = 0.0237).

Validation of a targeted next-generation sequencing panel for pancreatic ductal adenocarcinomas.

Pancreatic ductal adenocarcinoma (PDAC) is reported to be amongst the cancers with the lowest survival rate at 5 years. In the present study we aimed to validate a targeted next-generation sequencing (tNGS) panel to use in clinical routine, investigating genes important for PDAC diagnostic, prognostic and potential theragnostic aspect. In this NGS panel we also designed target regions to inquire about loss of heterozygosity (LOH) of chromosome 18 that has been described to be possibly linked to a worse disease progression. Copy number alteration has also been explored for a subset of genes. The last two methods are not commonly used for routine diagnostic with tNGS panels and we investigated their possible contribution to better characterize PDAC. A series of 140 formalin-fixed paraffin-embedded (FFPE) PDAC samples from 140 patients was characterized using this panel. Ninety-two % of patients showed alterations in at least one of the investigated genes (most frequent KRAS, TP53, SMAD4, CDKN2A and RNF43). Regarding LOH evaluation, we were able to detect chr18 LOH starting at 20% cell tumor percentage. The presence of LOH on chr18 is associated with a worse disease- and metastasis-free survival, in uni- and multivariate analyses. The present study validates the use of a tNGS panel for PDAC characterization, also evaluating chr18 LOH status for prognostic stratification.

De-novo assembly of four rail (Aves: Rallidae) genomes: A resource for comparative genomics.

Rails are a phenotypically diverse family of birds that includes 130 species and displays a wide distribution around the world. Here we present annotated genome assemblies for four rails from Aotearoa New Zealand: two native volant species, pukeko Porphyrio melanotus and mioweka Gallirallus philippensis, and two endemic flightless species takahe Porphyrio hochstetteri and weka Gallirallus australis. Using the sequence read data, heterozygosity was found to be lowest in the endemic flightless species and this probably reflects their relatively small populations. The quality checks and comparison with other rallid genomes showed that the new assemblies were of good quality. This study significantly increases the number of available rallid genomes and will enable future genomic studies on the evolution of this family.

The evolutionary cost of homophily: social stratification facilitates stable variant coexistence and increased rates of evolution in host-associated pathogens.

Coexistence of multiple strains of a pathogen in a host population can present significant challenges to vaccine development or treatment efficacy. Here we discuss a novel mechanism that can increase rates of long-lived strain polymorphism, rooted in the presence of social structure in a host population. We show that social preference of interaction, in conjunction with differences in immunity between host subgroups, can exert varying selection pressure on pathogen strains, creating a balancing mechanism that supports stable viral coexistence, independent of other known mechanisms. We use population genetic models to study rates of pathogen heterozygosity as a function of population size, host population composition, mutant strain fitness differences and host social preferences of interaction. We also show that even small periodic epochs of host population stratification can lead to elevated strain coexistence. These results are robust to varying social preferences of interaction, overall differences in strain fitnesses, and spatial heterogeneity in host population composition. Our results highlight the role of host population social stratification in increasing rates of pathogen strain diversity, with effects that should be considered when designing policies or treatments with a long-term view of curbing pathogen evolution.

Characteristics of Chinese breast cancer patients with double heterozygosity for BRCA1 and BRCA2 germline pathogenic variants.

PURPOSE: Despite of very rare, breast cancer patients with double heterozygosity (DH) variants in BRCA1 and BRCA2 genes have been identified in other ethnic groups and seem to be associated with distinctive phenotypes. However, little is known about the frequency and clinical characteristics of Chinese breast cancer patients with BRCA1/2 DH variants. METHODS: Four hundred and eleven unrelated patients with BRCA1 or BRCA2 pathogenic variants (PVs) were identified in a large series of unselected breast cancer patients. Another two siblings with metachronous bilateral breast cancer were referred for genetic counseling, after which BRCA1/2 DH variants were detected. RESULTS: Four unrelated breast cancer patients with BRCA1/2 DH were identified in the cohort of 411 patients with BRCA1 or BRCA2 PVs, the frequency of BRCA1/2 DH was 0.97%. In total, six BRCA1/2 DH patients from five families were found in this study. In two families, the hereditary pattern of DH was speculated to have originated from both sides of the family. BRCA1/2 DH patients were more likely to have a family history of breast cancer than patients with a BRCA1 (100% vs. 29.2%, P = 0.004) or BRCA2 (100% vs. 29.6%, P = 0.004) single PV. BRCA1/2 DH patients were more likely to be triple-negative breast tumors than patients with single BRCA2 PVs (66.7% vs. 14.1%, P = 0.020), which was comparable to the findings in patients with single BRCA1 PVs (66.7% vs. 56.9%, P = 1.00). CONCLUSION: Chinese patients with BRCA1/2 DH exhibit a high percentage of family history of breast cancer. The tumor pathological features of BRCA1/2 DH carriers are similar to those of BRCA1 PV carriers.

Whole-genome resequencing reveals genomic variation and dynamics in Ethiopian indigenous goats.

Ethiopia has about 52 million indigenous goats with marked phenotypic variability, which is the outcome of natural and artificial selection. Here, we obtained whole-genome sequence data of three Ethiopian indigenous goat populations (Arab, Fellata, and Oromo) from northwestern Ethiopia and analyzed their genome-wide genetic diversity, population structure, and signatures of selection. We included genotype data from four other Ethiopian goat populations (Abergelle, Keffa, Gumuz, and Woyto-Guji) and goats from Asia; Europe; and eastern, southern, western, and northern Africa to investigate the genetic predisposition of the three Ethiopian populations and performed comparative genomic analysis. Genetic diversity analysis showed that Fellata goats exhibited the lowest heterozygosity values (Ho = 0.288 ± 0.005 and He = 0.334 ± 0.0001). The highest values were observed in Arab goats (Ho = 0.310 ± 0.010 and He = 0.347 ± 4.35e-05). A higher inbreeding coefficient (FROH = 0.137 ± 0.016) was recorded for Fellata goats than the 0.105 ± 0.030 recorded for Arab and the 0.112 ± 0.034 recorded for Oromo goats. This indicates that the Fellata goat population should be prioritized in future conservation activities. The three goat populations showed the majority (∼63%) of runs of homozygosity in the shorter (100-150 Kb) length category, illustrating ancient inbreeding and/or small founder effects. Population relationship and structure analysis separated the Ethiopian indigenous goats into two distinct genetic clusters lacking phylogeographic structure. Arab, Fellata, Oromo, Abergelle, and Keffa represented one genetic cluster. Gumuz and Woyto-Guji formed a separate cluster and shared a common genetic background with the Kenyan Boran goat. Genome-wide selection signature analysis identified nine strongest regions spanning 163 genes influencing adaptation to arid and semi-arid environments (HOXC12, HOXC13, HOXC4, HOXC6, and HOXC9, MAPK8IP2), immune response (IL18, TYK2, ICAM3, ADGRG1, and ADGRG3), and production and reproduction (RARG and DNMT1). Our results provide insights into a thorough understanding of genetic architecture underlying selection signatures in Ethiopian indigenous goats in a semi-arid tropical environment and deliver valuable information for goat genetic improvement, conservation strategy, genome-wide association study, and marker-assisted breeding.

Molecular landscape of borderline ovarian tumours: A systematic review.

Borderline ovarian tumours (BOTs) show intriguing characteristics distinguishing them from other ovarian tumours. The aim of the systematic review was to analyse the spectrum of molecular changes found in BOTs and discuss their significance in the context of the overall therapeutic approach. The systematic review included articles published between 2000 and 2023 in the databases: PubMed, EMBASE, and Cochrane. After a detailed analysis of the available publications, we qualified for the systematic review: 28 publications on proto-oncogenes: BRAF, KRAS, NRAS, ERBB2, and PIK3CA, 20 publications on tumour suppressor genes: BRCA1/2, ARID1A, CHEK2, PTEN, 4 on adhesion molecules: CADM1, 8 on proteins: B-catenin, claudin-1, and 5 on glycoproteins: E-Cadherin. In addition, in the further part of the systematic review, we included eight publications on microsatellite instability and three describing loss of heterozygosity in BOT. Molecular changes found in BOTs can vary on a case-by-case basis, identifying carcinogenic mutations through molecular analysis and developing targeted therapies represent significant advancements in the diagnosis and treatment of ovarian malignancies. Molecular studies have contributed significantly to our understanding of BOT pathogenesis, but substantial research is still required to elucidate the relationship between ovarian neoplasms and extraneous disease, identify accurate prognostic indicators, and develop targeted therapeutic approaches.

Allelic loss of HLA class I facilitates evasion from immune surveillance in cervical intraepithelial neoplasia.

Loss of heterozygosity (LOH) has been reported to occur in HLA regions in cervical intraepithelial neoplasia (CIN) and cervical cancer. However, the details of how this is related to the progression of CIN have been unclear. In this study, we examined the human papillomavirus (HPV) antigen-presenting capacity of people with CIN and the significance of LOH of HLA class I in the progression of CIN. It was shown that differences in antigen-presenting capacity among each case depended on HLA types, not HPV genotypes. Focusing on the HLA type, there was a positive correlation between antigen-presenting capacity against HPV and the frequency of allelic loss. Furthermore, the lost HLA-B alleles had a higher HPV antigen-presenting capacity than intact alleles. In addition, frequency of LOH of HLA class I was significantly higher in advanced CIN (CIN2-3) than in cervicitis or early-stage CIN (CIN1): around half of CIN2-3 had LOH of any HLA class I. Moreover, the antigen-presenting capacity against E5, which is the HPV proteins that facilitate viral escape from this immune surveillance by suppressing HLA class I expression, had the most significant impact on the LOH in HLA-B. This study suggests that HPV evades immune surveillance mechanisms when host cells lose the capacity for antigen presentation by HLA class I molecules, resulting in long-term infection and progression to advanced lesions.

What does effective population size tell us about loss of allelic variation?

There are two primary measures of the amount of genetic variation in a population at a locus: heterozygosity and the number of alleles. Effective population size (N e) provides both an expectation of the amount of heterozygosity in a population at drift-mutation equilibrium and the rate of loss of heterozygosity because of genetic drift. In contrast, the number of alleles in a population at drift-mutation equilibrium is a function of both N e and census size (N C). In addition, populations with the same N e can lose allelic variation at very different rates. Allelic variation is generally much more sensitive to bottlenecks than heterozygosity. Expressions used to adjust for the effects of violations of the ideal population on N e do not provide good predictions of the loss of allelic variation. These effects are much greater for loci with many alleles, which are often important for adaptation. We show that there is a linear relationship between the reduction of N C and the corresponding reduction of the expected number of alleles at drift-mutation equilibrium. This makes it possible to predict the expected effect of a bottleneck on allelic variation. Heterozygosity provides good estimates of the rate of adaptive change in the short-term, but allelic variation provides important information about long-term adaptive change. The guideline of long-term N e being greater than 500 is often used as a primary genetic metric for evaluating conservation status. We recommend that this guideline be expanded to take into account allelic variation as well as heterozygosity.

Cancer Sample Analysis Utilizing Single-Nucleotide Polymorphism Array and Array Comparative Genomic Hybridization.

Chromosomal microarray, including single-nucleotide polymorphism (SNP) array and array comparative genomic hybridization (aCGH), enables the detection of DNA copy-number loss and/or gain associated with unbalanced chromosomal aberrations. In addition, SNP array and aCGH with SNP component also detect copy-neutral loss of heterozygosity (CN-LOH). Here we describe the chromosomal microarray procedure from the sample preparation using extracted DNA to the scanning of the array chip.

Genomic alteration discordance in the paired primary-recurrent ovarian cancers: based on the comprehensive genomic profiling (CGP) analysis.

PURPOSE: Ovarian cancer (OC) is characterized by a high recurrence rate, and homologous recombination deficiency (HRD) is an important biomarker in the clinical management of OC. We investigated the differences in clinical genomic profiles between the primary and platinum-sensitive recurrent OC (PSROC), focusing on HRD status. MATERIALS AND METHODS: A total of 40 formalin-fixed paraffin-embedded (FFPE) tissues of primary tumors and their first platinum-sensitive recurrence from 20 OC patients were collected, and comprehensive genomic profiling (CGP) analysis of FoundationOne®CDx (F1CDx) was applied to explore the genetic (dis)similarities of the primary and recurrent tumors. RESULTS: By comparing between paired samples, we found that genomic loss of heterozygosity (gLOH) score had a high intra-patient correlation (r2 = 0.79) and that short variants (including TP53, BRCA1/2 and NOTCH1 mutations), tumor mutational burden (TMB) and microsatellite stability status remained stable. The frequency of (likely) pathological BRCA1/2 mutations was 30% (12/40) in all samples positively correlated with gLOH scores, but the proportion of gLOH-high status (score > 16%) was 50% (10/20) and 55% (11/20) in the primary and recurrent samples, respectively. An additional 20% (4/20) of patients needed attention, a quarter of which carried the pathological BRCA1 mutation but had a gLOH-low status (gLOH < 16%), and three-quarters had different gLOH status in primary-recurrent pairs. Furthermore, we observed the PSROC samples had higher gLOH scores (16.1 ± 9.24 vs. 19.4 ± 11.1, p = 0.007), more CNVs (36.1% vs. 15.1% of discordant genomic alternations), and significant enrichment of altered genes in TGF-beta signaling and Hippo signaling pathways (p < 0.05 for all) than their paired primaries. Lastly, mutational signature and oncodrive gene analyses showed that the computed mutational signature similarity in the primary and recurrent tumors were best matched the COSMI 3 signature (Aetiology of HRD) and had consistent candidate cancer driver genes of MSH2, NOTCH1 and MSH6. CONCLUSION: The high genetic concordance of the short variants remains stable along OC recurrence. However, the results reveal significantly higher gLOH scores in the recurrent setting than in paired primaries, supporting further clinically instantaneity HRD assay strategy.

Digenic Inheritance of Mutations in Homologous Recombination Genes in Cancer Patients.

BACKGROUND/OBJECTIVES: BRCA1, BRCA2, ATM, and CHEK2 are known cancer predisposition genes (CPGs), but tumor risk in patients with simultaneous pathogenic variants (PVs) in CPGs remains largely unknown. In this study, we describe six patients from five families with multiple cancers who coinherited a combination of PVs in these genes. METHODS: PVs were identified using NGS DNA sequencing and were confirmed by Sanger. RESULTS: Families 1, 2, and 3 presented PVs in BRCA2 and ATM, family 4 in BRCA2 and BRCA1, and family 5 in BRCA2 and CHEK2. PVs were identified using NGS DNA sequencing and were confirmed by Sanger. The first family included patients with kidney, prostate, and breast cancer, in addition to pancreatic adenocarcinomas. In the second family, a female had breast cancer, while a male from the third family had prostate, gastric, and pancreatic cancer. The fourth family included a male with pancreatic cancer, and the fifth family a female with breast cancer. CONCLUSIONS: The early age of diagnosis and the development of multiple cancers in the reported patients indicate a very high risk of cancer in double-heterozygous patients associated with PVs in HR-related CPGs. Therefore, in families with patients who differ from other family members in terms of phenotype, age of diagnosis, or type of cancer, the cascade testing needs to include the study of other CPGs.

Application of Chromosomal Microarray Analysis in Genetic Reasons of Miscarriage Tissues.

Background: The potential causes of miscarriage are very complex, including genetic, immune, infectious, and endocrine factors. 50%-60% of miscarriages are caused by chromosomal abnormalities. Chromosomal microarray analysis (CMA) is a key tool in this context, capable of detecting not only copy number variations (CNV) but also loss of heterozygosity (LOH). CMA has been used as a tool to investigate the genetic reasons for miscarriage. Methods: In our study, chromosomal microarray analysis (CMA) conducted 1220 miscarriage villous tissues. The results from this technology were used to identify the genetic reasons for miscarriage and evaluated strategies for subsequent pre-pregnancy planning. Results: Here, the abnormality rate of miscarriage was 56.07%(684/1220). The aneuploidy rate accounted for 81.14%(555/684), and was significantly higher in group >35-year-old age. The second most common genetic reason for miscarriage was polyploidy, accounting for 10.09%(69/684). Additionally, we discovered loss of heterozygosity (LOH) in a small percentage of cases, accounting for 2.20%(15/684) reason for miscarriage genetic reasons, due to the advantage of CMA can detect isodisomy (a kind of uniparental disomy). 45 cases (6.58%) with copy number variants, which due to the CMA can detect copy number variations. Conclusion: Our study indicated that miscarriage villous tissues should be performed genetic analysis, seek help from professional genetic counseling.

Probing gene function in Candida albicans wild-type strains by Cas9-facilitated one-step integration of two dominant selection markers: a systematic analysis of recombination events at the target locus.

The adaptation of gene deletion methods based on the CRISPR-Cas9 system has facilitated the genetic manipulation of the pathogenic yeast Candida albicans, because homozygous mutants of this diploid fungus can now be generated in a single step, allowing the rapid screening of candidate genes for their involvement in a phenotype of interest. However, the Cas9-mediated double-strand breaks at the target site may result in an undesired loss of heterozygosity (LOH) on the affected chromosome and cause phenotypic alterations that are not related to the function of the investigated gene. In our present study, we harnessed Cas9-facilitated gene deletion to probe a set of genes that are constitutively overexpressed in strains containing hyperactive forms of the transcription factor Mrr1 for a possible contribution to the fluconazole resistance of such strains. To this aim, we used gene deletion cassettes containing two different dominant selection markers, caSAT1 and HygB, which confer resistance to nourseothricin and hygromycin, respectively, for simultaneous genomic integration in a single step, hypothesizing that this would minimize undesired LOH events at the target locus. We found that selection for resistance to both nourseothricin and hygromycin strongly increased the proportion of homozygous deletion mutants among the transformants compared with selection on media containing only one of the antibiotics, but it did not avoid undesired LOH events. Our results demonstrate that LOH on the target chromosome is a significant problem when using Cas9 for the generation of C. albicans gene deletion mutants, which demands a thorough examination of recombination events at the target site. IMPORTANCE: Candida albicans is one of the medically most important fungi and a model organism to study fungal pathogenicity. Investigating gene function in this diploid yeast has been facilitated by the adaptation of gene deletion methods based on the bacterial CRISPR-Cas9 system, because they enable the generation of homozygous mutants in a single step. We found that, in addition to increasing the efficiency of gene replacement by selection markers, the Cas9-mediated double-strand breaks also result in frequent loss of heterozygosity on the same chromosome, even when two different selection markers were independently integrated into the two alleles of the target gene. Since loss of heterozygosity for other genes can result in phenotypic alterations that are not caused by the absence of the target gene, these findings show that it is important to thoroughly analyze recombination events at the target locus when using Cas9 to generate gene deletion mutants in C. albicans.

SLRfinder: A method to detect candidate sex-linked regions with linkage disequilibrium clustering.

Despite their critical roles in genetic sex determination, sex chromosomes remain unknown in many non-model organisms, especially those having recently evolved sex-linked regions (SLRs). These evolutionarily young and labile sex chromosomes are important for understanding early sex chromosome evolution but are difficult to identify due to the lack of Y/W degeneration and SLRs limited to small genomic regions. Here, we present SLRfinder, a method to identify candidate SLRs using linkage disequilibrium (LD) clustering, heterozygosity and genetic divergence. SLRfinder does not rely on specific sequencing methods or a specific type of reference genome (e.g., from the homomorphic sex). In addition, the input of SLRfinder does not require phenotypic sexes, which may be unknown from population sampling, but sex information can be incorporated and is necessary to validate candidate SLRs. We tested SLRfinder using various published datasets and compared it to the local principal component analysis (PCA) method and the depth-based method Sex Assignment Through Coverage (SATC). As expected, the local PCA method could not be used to identify unknown SLRs. SATC works better on conserved sex chromosomes, whereas SLRfinder outperforms SATC in analysing labile sex chromosomes, especially when SLRs harbour inversions. Power analyses showed that SLRfinder worked better when sampling more populations that share the same SLR. If analysing one population, a relatively larger sample size (around 50) is needed for sufficient statistical power to detect significant SLR candidates, although true SLRs are likely always top-ranked. SLRfinder provides a novel and complementary approach for identifying SLRs and uncovering additional sex chromosome diversity in nature.

Exome-wide evidence of compound heterozygous effects across common phenotypes in the UK Biobank.

The phenotypic impact of compound heterozygous (CH) variation has not been investigated at the population scale. We phased rare variants (MAF ∼0.001%) in the UK Biobank (UKBB) exome-sequencing data to characterize recessive effects in 175,587 individuals across 311 common diseases. A total of 6.5% of individuals carry putatively damaging CH variants, 90% of which are only identifiable upon phasing rare variants (MAF < 0.38%). We identify six recessive gene-trait associations (p < 1.68 × 10-7) after accounting for relatedness, polygenicity, nearby common variants, and rare variant burden. Of these, just one is discovered when considering homozygosity alone. Using longitudinal health records, we additionally identify and replicate a novel association between bi-allelic variation in ATP2C2 and an earlier age at onset of chronic obstructive pulmonary disease (COPD) (p < 3.58 × 10-8). Genetic phase contributes to disease risk for gene-trait pairs: ATP2C2-COPD (p = 0.000238), FLG-asthma (p = 0.00205), and USH2A-visual impairment (p = 0.0084). We demonstrate the power of phasing large-scale genetic cohorts to discover phenome-wide consequences of compound heterozygosity.