Integer programming for selecting set of informative markers in paternity inference

BACKGROUND: Parentage information is fundamental to various life sciences. Recent advances in sequencing technologies have made it possible to accurately infer parentage even in non-model species. The optimization of sets of genome-wide markers is valuable for cost-effective applications but requires extremely large amounts of computation, which presses for the development of new efficient algorithms. RESULTS: Here, for a closed half-sib population, we generalized the process of marker loci selection as a binary integer programming problem. The proposed systematic formulation considered marker localization and the family structure of the potential parental population, resulting in an accurate assignment with a small set of markers. We also proposed an efficient heuristic approach, which effectively improved the number of markers, localization, and tolerance to missing data of the set. Applying this method to the actual genotypes of apple (Malus × domestica) germplasm, we identified a set of 34 SNP markers that distinguished 300 potential parents crossed to a particular cultivar with a greater than 99% accuracy. CONCLUSIONS: We present a novel approach for selecting informative markers based on binary integer programming. Since the data generated by high-throughput sequencing technology far exceeds the requirement for parentage assignment, a combination of the systematic marker selection with targeted SNP genotyping, such as KASP, allows flexibly enlarging the analysis up to a scale that has been unrealistic in various species. The method developed in this study can be directly applied to unsolved large-scale problems in breeding, reproduction, and ecological research, and is expected to lead to novel knowledge in various biological fields. The implementation is available at https://github.com/SoNishiyama/IP-SIMPAT

Genomic insight into the developmental history of southern highbush blueberry populations

ハイブッシュブルーベリーに暖地適応性をもたらした遺伝要因を解明 --ゲノムに刻まれたブルーベリーの育種履歴--. 京都大学プレスリリース. 2020-09-07.Interspecific hybridization is a common breeding approach for introducing novel traits and genetic diversity to breeding populations. Southern highbush blueberry (SHB) is a blueberry cultivar group that has been intensively bred over the last 60 years. Specifically, it was developed by multiple interspecific crosses between northern highbush blueberry [NHB, Vaccinium corymbosum L. (2n = 4x = 48)] and low-chill Vaccinium species to expand the geographic limits of highbush blueberry production. In this study, we genotyped polyploid blueberries, including 105 SHB, 17 NHB, and 10 rabbiteye blueberry (RE) (Vaccinium virgatum Aiton), from the accessions planted at Poplarville, Mississippi, and accessions distributed in Japan, based on the double-digest restriction site-associated DNA sequencing. The genome-wide SNP data clearly indicated that RE cultivars were genetically distinct from SHB and NHB cultivars, whereas NHB and SHB were genetically indistinguishable. The population structure results appeared to reflect the differences in the allele selection strategies that breeders used for developing germplasm adapted to local climates. The genotype data implied that there are no or very few genomic segments that were commonly introgressed from low-chill Vaccinium species to the SHB genome. Principal component analysis-based outlier detection analysis found a few loci associated with a variable that could partially differentiate NHB and SHB. These SNP loci were detected in Mb-scale haplotype blocks and may be close to the functional genes related to SHB development. Collectively, the data generated in this study suggest a polygenic adaptation of SHB to the southern climate, and may be relevant for future population-scale genome-wide analyses of blueberry

Genomics-based discrimination of 2n gamete formation mechanisms in polyploids: a case study in nonaploid Diospyros kaki ‘Akiou’

Unreduced gametes (2n gametes), possessing double the haploid genome, whatever ploidy that happens to be, are a common source of ploidy variation in plant populations. First and second division restitution (FDR and SDR) are the dominant mechanisms of 2n gamete production; all else being equal, FDR gametes have a higher degree of heterozygosity, thus they are advantageous in breeding. The discrimination of these mechanisms from the consequence of hybridization is challenging, especially in higher polyploids, and usually requires information on centromere location. In this study, we propose a genotyping-based strategy to uncover the mechanisms of 2n gamete formation in progeny that has a higher ploidy than its parents. Simulation of 2n gamete production revealed that FDR and SDR pathways can be discriminated based on allele transmission patterns alone without information on centromere location. We applied this strategy to study the formation mechanism of a nonaploid Diospyros kaki ‘Akiou', which was bred via hybridization between D. kaki hexaploid cultivars. The result demonstrated that ‘Akiou' was derived from the fertilization of a normal female gamete by a 2n male gamete and that this 2n gamete was produced through FDR. Consequently, the distinct duplex transmission pattern in the FDR gamete enabled us to infer the genomic characteristics of polyploid persimmon. The method could be tested only for the plant being polypoid, which allows for the ability to discriminate causes of 2n gamete formation using allele dosage in progeny, and will be useful in future studies of polyploid genomics

HPF1-dependent PARP activation promotes LIG3-XRCC1-mediated backup pathway of Okazaki fragment ligation

DNA ligase 1 (LIG1) is known as the major DNA ligase responsible for Okazaki fragment joining. Recent studies have implicated LIG3 complexed with XRCC1 as an alternative player in Okazaki fragment joining in cases where LIG1 is not functional, although the underlying mechanisms are largely unknown. Here, using a cell-free system derived from Xenopus egg extracts, we demonstrated the essential role of PARP1-HPF1 in LIG3-dependent Okazaki fragment joining. We found that Okazaki fragments were eventually ligated even in the absence of LIG1, employing in its place LIG3-XRCC1, which was recruited onto chromatin. Concomitantly, LIG1 deficiency induces ADP-ribosylation of histone H3 in a PARP1-HPF1-dependent manner. The depletion of PARP1 or HPF1 resulted in a failure to recruit LIG3 onto chromatin and a subsequent failure in Okazaki fragment joining in LIG1-depleted extracts. Importantly, Okazaki fragments were not ligated at all when LIG1 and XRCC1 were co-depleted. Our results suggest that a unique form of ADP-ribosylation signaling promotes the recruitment of LIG3 on chromatin and its mediation of Okazaki fragment joining as a backup system for LIG1 perturbation

Functional and expressional analyses of apple FLC-like in relation to dormancy progress and flower bud development.

We previously identified the FLOWERING LOCUS C (FLC)-like gene, a MADS-box transcription factor gene that belongs to Arabidopsis thaliana L. FLC clade, in apple (Malus ×domestica Borkh.), and its expression in dormant flower buds is positively correlated with cumulative cold exposure. To elucidate the role of the MdFLC-like in the dormancy process and flower development, we first characterized the phenotypes of MdFLC-like overexpressing lines with the Arabidopsis Columbia-0 background. The overexpression of MdFLC-like significantly delayed the bolting date and reduced the plant size, but it did not significantly affect the number of rosette leaves or flower organ formation. Thus, MdFLC-like may affect vegetative growth and development rather than flowering when expressed in Arabidopsis, which is not like Arabidopsis FLC that affects development of flowering. We compared seasonal expression patterns of MdFLC-like in low-chill ‘Anna’ and high-chill ‘Fuji’ and ‘Tsugaru’ apples collected from trees grown in a cold winter region in temperate zone and found an earlier upregulation in ‘Anna’ compared with ‘Fuji’ and ‘Tsugaru’. Expression patterns were also compared in relation to developmental changes in the flower primordia during the chilling accumulation period. Overall, MdFLC-like was progressively upregulated during flower primordia differentiation and development in autumn to early winter and reached a maximum expression level at around the same time as the genotype-dependent chilling requirements were fulfilled in high-chill cultivars. Thus, we hypothesize MdFLC-like may be upregulated in response to cold exposure and flower primordia development during the progress of endodormancy. Our study also suggests MdFLC-like may have a growth-inhibiting function during the end of endodormancy and ecodormancy when the temperature is low and unfavorable for rapid bud outgrowth

Human intestinal spirochetosis accompanied by human immunodeficiency virus infection:a case report

We present a middle-aged, heterosexual Japanese man with mixed infections including human intestinal spirochetosis, which led us to the detection of human immunodeficiency virus (HIV) infection. The patient had syphilis without related physical or neurological findings. An examination for the serum antibody for HIV performed 9 years previously was negative. In a complete medical checkup at the present time, human intestinal spirochetosis and unspecified entamebic cysts were suggested by histological examination of colonic biopsy material and parasitic examination of the intestinal fluid, respectively. Moreover, a serological test for the antibody for HIV was positive. In specimens obtained by colonoscopy, Brachyspira aalborgi was diagnosed by ultrastructural study and the polymerase chain reaction method for bacterial 16S ribosomal deoxyribonucleic acid. Although HIV infection remains at low prevalence in Japan, we recommend examination for HIV infection in patients with human intestinal spirochetosis, especially when other co-infections are apparent.</p

Two distinct modes of DNMT1 recruitment ensure stable maintenance DNA methylation

Stable inheritance of DNA methylation is critical for maintaining differentiated phenotypes in multicellular organisms. We have recently identified dual mono-ubiquitylation of histone H3 (H3Ub2) by UHRF1 as an essential mechanism to recruit DNMT1 to chromatin. Here, we show that PCNA-associated factor 15 (PAF15) undergoes UHRF1-dependent dual mono-ubiquitylation (PAF15Ub2) on chromatin in a DNA replication-coupled manner. This event will, in turn, recruit DNMT1. During early S-phase, UHRF1 preferentially ubiquitylates PAF15, whereas H3Ub2 predominates during late S-phase. H3Ub2 is enhanced under PAF15 compromised conditions, suggesting that H3Ub2 serves as a backup for PAF15Ub2. In mouse ES cells, loss of PAF15Ub2 results in DNA hypomethylation at early replicating domains. Together, our results suggest that there are two distinct mechanisms underlying replication timing-dependent recruitment of DNMT1 through PAF15Ub2 and H3Ub2, both of which are prerequisite for high fidelity DNA methylation inheritance

Association between the presence of bacteria in prostate tissue and histopathology in biopsies from men not complaining of lower urinary tract symptoms

Objective: To investigate the presence of bacteria in prostate tissue, and relationships between the bacteria and histopathological findings. Methods: Samples were collected from prostate biopsy patients with no obvious lower urinary tract symptoms (LUTS). Detection and identification of bacterial species in the prostate tissues were performed with PCR for 16SrDNA and DNA sequencing. Histopathology was also evaluated. LUTS and lower urinary tract function were assessed by questionnaires, uroflowmetry, and ultrasonography. Results: DNA was extracted from 97 prostate biopsies, with 5 bacterial species detected among samples from 7 patients (7.2%). The stroma-to-gland ratio in the prostate tissues from patients with bacteria was lower than in those without bacteria (p < 0.01). Glandular epithelial hyperplasia was also identified in the prostates harboring bacteria. International Prostate Symptom Score (IPSS), IPSS-quality of life (IPSS-QOL), Overactive Bladder Symptom Score (OABSS), maximum flow rate, urine volume by uroflowmetry, and post-voided residual urine were not significantly different when comparing patients with and without bacteria in their prostate samples. Conclusions: The present study demonstrated that 7.2% of men without obvious LUTS had bacteria in their prostate tissues. The presence of such bacteria might induce glandular hyperplasia and contribute to pathological changes in the early stages of benign prostate enlargement before affecting LUTS

The whole blood transcriptional regulation landscape in 465 COVID-19 infected samples from Japan COVID-19 Task Force

「コロナ制圧タスクフォース」COVID-19患者由来の血液細胞における遺伝子発現の網羅的解析 --重症度に応じた遺伝子発現の変化には、ヒトゲノム配列の個人差が影響する--. 京都大学プレスリリース. 2022-08-23.Coronavirus disease 2019 (COVID-19) is a recently-emerged infectious disease that has caused millions of deaths, where comprehensive understanding of disease mechanisms is still unestablished. In particular, studies of gene expression dynamics and regulation landscape in COVID-19 infected individuals are limited. Here, we report on a thorough analysis of whole blood RNA-seq data from 465 genotyped samples from the Japan COVID-19 Task Force, including 359 severe and 106 non-severe COVID-19 cases. We discover 1169 putative causal expression quantitative trait loci (eQTLs) including 34 possible colocalizations with biobank fine-mapping results of hematopoietic traits in a Japanese population, 1549 putative causal splice QTLs (sQTLs; e.g. two independent sQTLs at TOR1AIP1), as well as biologically interpretable trans-eQTL examples (e.g., REST and STING1), all fine-mapped at single variant resolution. We perform differential gene expression analysis to elucidate 198 genes with increased expression in severe COVID-19 cases and enriched for innate immune-related functions. Finally, we evaluate the limited but non-zero effect of COVID-19 phenotype on eQTL discovery, and highlight the presence of COVID-19 severity-interaction eQTLs (ieQTLs; e.g., CLEC4C and MYBL2). Our study provides a comprehensive catalog of whole blood regulatory variants in Japanese, as well as a reference for transcriptional landscapes in response to COVID-19 infection

DOCK2 is involved in the host genetics and biology of severe COVID-19

「コロナ制圧タスクフォース」COVID-19疾患感受性遺伝子DOCK2の重症化機序を解明 --アジア最大のバイオレポジトリーでCOVID-19の治療標的を発見--. 京都大学プレスリリース. 2022-08-10.Identifying the host genetic factors underlying severe COVID-19 is an emerging challenge. Here we conducted a genome-wide association study (GWAS) involving 2, 393 cases of COVID-19 in a cohort of Japanese individuals collected during the initial waves of the pandemic, with 3, 289 unaffected controls. We identified a variant on chromosome 5 at 5q35 (rs60200309-A), close to the dedicator of cytokinesis 2 gene (DOCK2), which was associated with severe COVID-19 in patients less than 65 years of age. This risk allele was prevalent in East Asian individuals but rare in Europeans, highlighting the value of genome-wide association studies in non-European populations. RNA-sequencing analysis of 473 bulk peripheral blood samples identified decreased expression of DOCK2 associated with the risk allele in these younger patients. DOCK2 expression was suppressed in patients with severe cases of COVID-19. Single-cell RNA-sequencing analysis (n = 61 individuals) identified cell-type-specific downregulation of DOCK2 and a COVID-19-specific decreasing effect of the risk allele on DOCK2 expression in non-classical monocytes. Immunohistochemistry of lung specimens from patients with severe COVID-19 pneumonia showed suppressed DOCK2 expression. Moreover, inhibition of DOCK2 function with CPYPP increased the severity of pneumonia in a Syrian hamster model of SARS-CoV-2 infection, characterized by weight loss, lung oedema, enhanced viral loads, impaired macrophage recruitment and dysregulated type I interferon responses. We conclude that DOCK2 has an important role in the host immune response to SARS-CoV-2 infection and the development of severe COVID-19, and could be further explored as a potential biomarker and/or therapeutic target