The third international hackathon for applying insights into large-scale genomic composition to use cases in a wide range of organisms.

In October 2021, 59 scientists from 14 countries and 13 U.S. states collaborated virtually in the Third Annual Baylor College of Medicine & DNANexus Structural Variation hackathon. The goal of the hackathon was to advance research on structural variants (SVs) by prototyping and iterating on open-source software. This led to nine hackathon projects focused on diverse genomics research interests, including various SV discovery and genotyping methods, SV sequence reconstruction, and clinically relevant structural variation, including SARS-CoV-2 variants. Repositories for the projects that participated in the hackathon are available at https://github.com/collaborativebioinformatics.

Measurement of SARS-CoV-2 Antibody Titers Improves the Prediction Accuracy of COVID-19 Maximum Severity by Machine Learning in Non-Vaccinated Patients.

Numerous studies have suggested that the titers of antibodies against SARS-CoV-2 are associated with the COVID-19 severity, however, the types of antibodies associated with the disease maximum severity and the timing at which the associations are best observed, especially within one week after symptom onset, remain controversial. We attempted to elucidate the antibody responses against SARS-CoV-2 that are associated with the maximum severity of COVID-19 in the early phase of the disease, and to investigate whether antibody testing might contribute to prediction of the disease maximum severity in COVID-19 patients. We classified the patients into four groups according to the disease maximum severity (severity group 1 (did not require oxygen supplementation), severity group 2a (required oxygen supplementation at low flow rates), severity group 2b (required oxygen supplementation at relatively high flow rates), and severity group 3 (required mechanical ventilatory support)), and serially measured the titers of IgM, IgG, and IgA against the nucleocapsid protein, spike protein, and receptor-binding domain of SARS-CoV-2 until day 12 after symptom onset. The titers of all the measured antibody responses were higher in severity group 2b and 3, especially severity group 2b, as early as at one week after symptom onset. Addition of data obtained from antibody testing improved the ability of analysis models constructed using a machine learning technique to distinguish severity group 2b and 3 from severity group 1 and 2a. These models constructed with non-vaccinated COVID-19 patients could not be applied to the cases of breakthrough infections. These results suggest that antibody testing might help physicians identify non-vaccinated COVID-19 patients who are likely to require admission to an intensive care unit.

Response kinetics of different classes of antibodies to SARS-CoV2 infection in the Japanese population: The IgA and IgG titers increased earlier than the IgM titers.

To better understand the immune responses to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in individuals with COVID-19, it is important to investigate the kinetics of the antibody responses and their associations with the clinical course in different populations, since there seem to be considerable differences between Western and Asian populations in the clinical features and spread of COVID-19. In this study, we serially measured the serum titers of IgM, IgG and IgA antibodies generated against the nucleocapsid protein (NCP), S1 subunit of the spike protein (S1), and receptor-binding domain in the S1 subunit (RBD) of SARS-CoV-2 in Japanese individuals with COVID-19. Among the IgM, IgG, and IgA antibodies, IgA antibodies against all of the aforementioned viral proteins were the first to appear after the infection, and IgG and/or IgA seroconversion often preceded IgM seroconversion. In regard to the timeline of the antibody responses to the different viral proteins (NCP, S1 and RBD), IgA against NCP appeared than IgA against S1 or RBD, while IgM and IgG against S1 appeared earlier than IgM/IgG against NCP or RBD. The IgG responses to all three viral proteins and responses of all three antibody classes to S1 and RBD were sustained for longer durations than the IgA/IgM responses to all three viral proteins and responses of all three antibody classes to NCP, respectively. The seroconversion of IgA against NCP occurred later and less frequently in patients with mild COVID-19. These results suggest possible differences in the antibody responses to SARS-CoV-2 antigens between the Japanese and Western populations.

Association of the Serum Levels of the Nucleocapsid Antigen of SARS-CoV-2 With the Diagnosis, Disease Severity, and Antibody Titers in Patients With COVID-19: A Retrospective Cross-Sectional Study.

Background: Several types of laboratory tests for COVID-19 have been established to date; however, the clinical significance of the serum SARS-CoV-2 nucleocapsid (N) antigen levels remains to be fully elucidated. In the present study, we attempted to elucidate the usefulness and clinical significance of the serum N antigen levels. Methods: We measured the serum N antigen levels in 391 serum samples collected from symptomatic patients with a confirmed diagnosis of COVID-19 and 96 serum samples collected from patients with non-COVID-19, using a fully automated chemiluminescence immunoassay analyzer. Results: Receiver operating characteristic analysis identified the optimal cutoff value of the serum N antigen level (cutoff index, based on Youden's index) as 0.255, which yielded a sensitivity and specificity for the diagnosis of COVID-19 of 91.0 and 81.3%, respectively. The serum N antigen levels were significantly higher in the patient groups with moderate and severe COVID-19 than with mild disease. Moreover, a significant negative correlation was observed between the serum N antigen levels and the SARS-CoV-2 IgG antibody titers, especially in patients with severe COVID-19. Conclusion: Serum N antigen testing might be useful both for the diagnosis of COVID-19 and for obtaining a better understanding of the clinical features of the disease.

Effects of long-term intake of a yogurt fermented with Lactobacillus delbrueckii subsp. bulgaricus 2038 and Streptococcus thermophilus 1131 on mice.

The gut is an extremely complicated ecosystem where micro-organisms, nutrients and host cells interact vigorously. Although the function of the intestine and its barrier system weakens with age, some probiotics can potentially prevent age-related intestinal dysfunction. Lactobacillus delbrueckii subsp. bulgaricus 2038 and Streptococcus thermophilus 1131, which are the constituents of LB81 yogurt, are representative probiotics. However, it is unclear whether their long-term intake has a beneficial influence on systemic function. Here, we examined the gut microbiome, fecal metabolites and gene expression profiles of various organs in mice. Although age-related alterations were apparent in them, long-term LB81 yogurt intake led to an increased Bacteroidetes to Firmicutes ratio and elevated abundance of the bacterial family S24-7 (Bacteroidetes), which is known to be associated with butyrate and propanoate production. According to our fecal metabolite analysis to detect enrichment, long-term LB81 yogurt intake altered the intestinal metabolic pathways associated with propanoate and butanoate in the mice. Gene ontology analysis also revealed that long-term LB81 yogurt intake influenced many physiological functions related to the defense response. The profiles of various genes associated with antimicrobial peptides-, tight junctions-, adherens junctions- and mucus-associated intestinal barrier functions were also drastically altered in the LB81 yogurt-fed mice. Thus, long-term intake of LB81 yogurt has the potential to maintain systemic homeostasis, such as the gut barrier function, by controlling the intestinal microbiome and its metabolites.

Promoter-level transcriptome in primary lesions of endometrial cancer identified biomarkers associated with lymph node metastasis.

For endometrial cancer patients, lymphadenectomy is recommended to exclude rarely metastasized cancer cells. This procedure is performed even in patients with low risk of recurrence despite the risk of complications such as lymphedema. A method to accurately identify cases with no lymph node metastases (LN-) before lymphadenectomy is therefore highly required. We approached this clinical problem by examining primary lesions of endometrial cancers with CAGE (Cap Analysis Gene Expression), which quantifies promoter-level expression across the genome. Fourteen profiles delineated distinct transcriptional networks between LN + and LN- cases, within those classified as having the low or intermediate risk of recurrence. Subsequent quantitative reverse transcription polymerase chain reaction (qRT-PCR) analyses of 115 primary tumors showed SEMA3D mRNA and TACC2 isoforms expressed through a novel promoter as promising biomarkers with high accuracy (area under the receiver operating characteristic curve, 0.929) when used in combination. Our high-resolution transcriptome provided evidence of distinct molecular profiles underlying LN + /LN- status in endometrial cancers, raising the possibility of preoperative diagnosis to reduce unnecessary operations in patients with minimum recurrence risk.

Novel biomarkers that assist in accurate discrimination of squamous cell carcinoma from adenocarcinoma of the lung.

BACKGROUND: Targeted therapies based on the molecular and histological features of cancer types are becoming standard practice. The most effective regimen in lung cancers is different between squamous cell carcinoma (SCC) and adenocarcinoma (AD). Therefore a precise diagnosis is crucial, but this has been difficult, particularly for poorly differentiated SCC (PDSCC) and AD without a lepidic growth component (non-lepidic AD). Biomarkers enabling a precise diagnosis are therefore urgently needed. METHODS: Cap Analysis of Gene Expression (CAGE) is a method used to quantify promoter activities across the whole genome by determining the 5' ends of capped RNA molecules with next-generation sequencing. We performed CAGE on 97 frozen tissues from surgically resected lung cancers (22 SCC and 75 AD), and confirmed the findings by immunohistochemical analysis (IHC) in an independent group (29 SCC and 45 AD). RESULTS: Using the genome-wide promoter activity profiles, we confirmed that the expression of known molecular markers used in IHC for SCC (CK5, CK6, p40 and desmoglein-3) and AD (TTF-1 and napsin A) were different between SCC and AD. We identified two novel marker candidates, SPATS2 for SCC and ST6GALNAC1 for AD, as showing comparable performance and complementary utility to the known markers in discriminating PDSCC and non-lepidic AD. We subsequently confirmed their utility at the protein level by IHC in an independent group. CONCLUSIONS: We identified two genes, SPATS2 and ST6GALNAC1, as novel complemental biomarkers discriminating SCC and AD. These findings will contribute to a more accurate diagnosis of NSCLC, which is crucial for precision medicine for lung cancer.

Discovery of molecular markers to discriminate corneal endothelial cells in the human body.

The corneal endothelium is a monolayer of hexagonal corneal endothelial cells (CECs) on the inner surface of the cornea. CECs are critical in maintaining corneal transparency through their barrier and pump functions. CECs in vivo have a limited capacity in proliferation, and loss of a significant number of CECs results in corneal edema called bullous keratopathy which can lead to severe visual loss. Corneal transplantation is the most effective method to treat corneal endothelial dysfunction, where it suffers from donor shortage. Therefore, regeneration of CECs from other cell types attracts increasing interests, and specific markers of CECs are crucial to identify actual CECs. However, the currently used markers are far from satisfactory because of their non-specific expression in other cell types. Here, we explored molecular markers to discriminate CECs from other cell types in the human body by integrating the published RNA-seq data of CECs and the FANTOM5 atlas representing diverse range of cell types based on expression patterns. We identified five genes, CLRN1, MRGPRX3, HTR1D, GRIP1 and ZP4 as novel markers of CECs, and the specificities of these genes were successfully confirmed by independent experiments at both the RNA and protein levels. Notably none of them have been documented in the context of CEC function. These markers could be useful for the purification of actual CECs, and also available for the evaluation of the products derived from other cell types. Our results demonstrate an effective approach to identify molecular markers for CECs and open the door for the regeneration of CECs in vitro.

RECLU: a pipeline to discover reproducible transcriptional start sites and their alternative regulation using capped analysis of gene expression (CAGE).

BACKGROUND: Next generation sequencing based technologies are being extensively used to study transcriptomes. Among these, cap analysis of gene expression (CAGE) is specialized in detecting the most 5' ends of RNA molecules. After mapping the sequenced reads back to a reference genome CAGE data highlights the transcriptional start sites (TSSs) and their usage at a single nucleotide resolution. RESULTS: We propose a pipeline to group the single nucleotide TSS into larger reproducible peaks and compare their usage across biological states. Importantly, our pipeline discovers broad peaks as well as the fine structure of individual transcriptional start sites embedded within them. We assess the performance of our approach on a large CAGE datasets including 156 primary cell types and two cell lines with biological replicas. We demonstrate that genes have complicated structures of transcription initiation events. In particular, we discover that narrow peaks embedded in broader regions of transcriptional activity can be differentially used even if the larger region is not. CONCLUSIONS: By examining the reproducible fine scaled organization of TSS we can detect many differentially regulated peaks undetected by previous approaches.

Meta-analysis identifies nine new loci associated with rheumatoid arthritis in the Japanese population.

Rheumatoid arthritis is a common autoimmune disease characterized by chronic inflammation. We report a meta-analysis of genome-wide association studies (GWAS) in a Japanese population including 4,074 individuals with rheumatoid arthritis (cases) and 16,891 controls, followed by a replication in 5,277 rheumatoid arthritis cases and 21,684 controls. Our study identified nine loci newly associated with rheumatoid arthritis at a threshold of P < 5.0 × 10(-8), including B3GNT2, ANXA3, CSF2, CD83, NFKBIE, ARID5B, PDE2A-ARAP1, PLD4 and PTPN2. ANXA3 was also associated with susceptibility to systemic lupus erythematosus (P = 0.0040), and B3GNT2 and ARID5B were associated with Graves' disease (P = 3.5 × 10(-4) and 2.9 × 10(-4), respectively). We conducted a multi-ancestry comparative analysis with a previous meta-analysis in individuals of European descent (5,539 rheumatoid arthritis cases and 20,169 controls). This provided evidence of shared genetic risks of rheumatoid arthritis between the populations.

Common variants at CDKAL1 and KLF9 are associated with body mass index in east Asian populations.

Obesity is a disorder with a complex genetic etiology, and its epidemic is a worldwide problem. Although multiple genetic loci associated with body mass index, the most common measure of obesity, have been identified in European populations, few studies have focused on Asian populations. Here we report a genome-wide association study and replication studies with 62,245 east Asian subjects, which identified two new body mass index-associated loci in the CDKAL1 locus at 6p22 (rs2206734, P = 1.4 × 10(-11)) and the KLF9 locus at 9q21 (rs11142387, P = 1.3 × 10(-9)), as well as several previously reported loci (the SEC16B, BDNF, FTO, MC4R and GIPR loci, P < 5.0 × 10(-8)). We subsequently performed gene-gene interaction analyses and identified an interaction (P = 2.0 × 10(-8)) between a SNP in the KLF9 locus (rs11142387) and one in the MSTN (also known as GDF8) locus at 2q32 (rs13034723). These findings should provide useful insights into the etiology of obesity.

Identification of nine novel loci associated with white blood cell subtypes in a Japanese population.

White blood cells (WBCs) mediate immune systems and consist of various subtypes with distinct roles. Elucidation of the mechanism that regulates the counts of the WBC subtypes would provide useful insights into both the etiology of the immune system and disease pathogenesis. In this study, we report results of genome-wide association studies (GWAS) and a replication study for the counts of the 5 main WBC subtypes (neutrophils, lymphocytes, monocytes, basophils, and eosinophils) using 14,792 Japanese subjects enrolled in the BioBank Japan Project. We identified 12 significantly associated loci that satisfied the genome-wide significance threshold of P<5.0×10(-8), of which 9 loci were novel (the CDK6 locus for the neutrophil count; the ITGA4, MLZE, STXBP6 loci, and the MHC region for the monocyte count; the SLC45A3-NUCKS1, GATA2, NAALAD2, ERG loci for the basophil count). We further evaluated associations in the identified loci using 15,600 subjects from Caucasian populations. These WBC subtype-related loci demonstrated a variety of patterns of pleiotropic associations within the WBC subtypes, or with total WBC count, platelet count, or red blood cell-related traits (n = 30,454), which suggests unique and common functional roles of these loci in the processes of hematopoiesis. This study should contribute to the understanding of the genetic backgrounds of the WBC subtypes and hematological traits.

Selecting the desirable method for predicting 305-day lactation yields in Mongolia.

In Mongolia, milk production must be improved both quantitatively and qualitatively, and we must find the most desirable method for predicting 305-day milk yields for cows in Mongolia. Therefore, the predicted 305-day milk yield from the test interval method (TIM), multiple-trait prediction (MTP) and the random regression model (RRM) were compared. The daily milk production records during 305 days for cows calving from 1986 to 2007 from the National Livestock Breeding Center in Japan, and the test day records every month for cows calving from 1985 to 2005 from the Hokkaido Dairy Milk Recording and Testing Association were used. Wilmink's model for the average of lactation curve in MTP was adopted. A cubic Legendre polynomials and the parameters of Wilmink's function were used for RRM. The results suggested that RRM was the most desirable method for predicting 305-day milk yields in Mongolia because estimates by RRM were the most accurate when using over four records. Further analyses will be required to examine the properties of RRM when predicting 305-day milk yields using test day records in Mongolia.

Genome-wide association study for C-reactive protein levels identified pleiotropic associations in the IL6 locus.

C-reactive protein (CRP) is a hallmark acute-phase reactant and is widely used as a blood marker for inflammation. Substantial roles of serum CRP levels in the pathogenesis of diseases have been suggested, and investigation of the mechanisms that regulate serum CRP levels would have a substantial clinical impact. Here, through genome-wide association and replication studies performed using 12 854 Japanese subjects, we identified a significant association between serum CRP levels and a single nucleotide polymorphism in the promoter region of interleukin-6 (IL6) (rs2097677, P = 4.1 × 10(-11)), a typical pleiotropic pro-inflammatory cytokine. Our study also replicated the associations in the CRP (rs3093059, P = 3.5 × 10(-21)) and HNF1A loci (rs7310409, P = 2.7 × 10(-8)). Pleiotropic association analysis with hematological and biochemical traits using 30 466 Japanese subjects demonstrated that the CRP-increasing allele of rs2097677 in the IL6 locus was significantly associated with an increased white blood cell count, platelet count and serum globulin and a decreased mean corpuscular hemoglobin and mean corpuscular hemoglobin concentration (P < 5.0 × 10(-4)), although no pleiotropic association was observed in the CRP or HNF1A locus (α = 0.01). Our study demonstrated the pivotal role of the IL6 locus in the regulation of serum CRP levels and inflammatory pathways.

A genome-wide association study in 19 633 Japanese subjects identified LHX3-QSOX2 and IGF1 as adult height loci.

Previous genome-wide association studies (GWASs) have identified several loci associated with human height; however, such evidence was mostly reported in Caucasian populations. Since the different distributions of height between populations suggest their different genetic backgrounds, analysis in different populations would be useful. Here, we present the results of a GWAS for adult height in 19 633 Japanese subjects. We found eight significantly associated loci that satisfied the genome-wide significance level (P < 5.0 x 10(-8)). Of these, the association to the LHX3-QSOX2 locus was entirely novel (rs12338076, P = 2.2 x 10(-8)). We also identified the association to the IGF1 locus (rs17032362, P = 8.1 x 10(-9)). Conditional association analysis in the IGF1 locus with rs17032362 suggested the existence of an additional independent association with height to this locus (rs1457595, P = 1.2 x 10(-5)). We observed large differences in the allele frequencies of rs17032362 and rs1457595 between Japanese (34 and 9%, respectively) and Caucasian (1.7 and 0%, respectively) populations, thereby suggesting weak statistical powers for the IGF1 locus in the previous Caucasian GWASs for height. We extensively compared our results with those of previous reports on the Caucasian and Korean populations. We were able to replicate all four loci previously reported in Koreans (EFEMP1, ZBTB38, HMGA1 and PLAG1, P < 5.0 x 10(-8)) and 15 loci identified in Caucasians (P < 0.001). The combination of the height-associated loci identified in our study and the previous GWASs demonstrated an effect size of 1.26 cm (95% confidence interval: 1.18-1.34) per 1.0 increase of the normalized Z score for height-increasing alleles, explaining 4.6% of the total variance of adult height.

Common variations in PSMD3-CSF3 and PLCB4 are associated with neutrophil count.

Neutrophils are the most abundant subtype of white blood cells (WBCs). Although the regulation of the numbers of neutrophils would have substantial clinical impacts, the studies on the variations associated with neutrophil count had not been performed further. To investigate genetic variations that regulate neutrophil count, we performed a genome-wide association study in 5771 Japanese subjects and a replication study using independent 1894 Japanese subjects. We identified two genetic loci significantly associated with neutrophil count (rs4794822 in PSMD3-CSF3 at 17q21.1, P = 6.3 x 10(-10); rs2072910 in PLCB4 at 20p12, P = 3.1 x 10(-10)). As these loci did not indicate significant associations with the counts of the other subtypes of WBCs (lymphocytes, monocytes, eosinophils and basophils), their specific associations with neutrophils were suggested. The combination of the single nucleotide polymorphisms (SNPs) in these two loci explained 1.0% of the total variance of the log-transformed values of the neutrophil count in our study populations. The subjects who were homozygous for 'neutrophil-increasing alleles' in both of the SNPs (T alleles for rs4794822 and rs2072910) had 1.17-fold (95% confidence interval: 1.10-1.24) higher neutrophil count when compared with the subjects homozygous for 'neutrophil-decreasing alleles' (C alleles for rs4794822 and rs2072910). In conclusion, our study would demonstrate the significant contribution of PSMD3-CSF3 and PLCB4 loci to the regulation of neutrophil count.

Structural basis for gibberellin recognition by its receptor GID1.

Gibberellins (GAs) are phytohormones essential for many developmental processes in plants. A nuclear GA receptor, GIBBERELLIN INSENSITIVE DWARF1 (GID1), has a primary structure similar to that of the hormone-sensitive lipases (HSLs). Here we analyse the crystal structure of Oryza sativa GID1 (OsGID1) bound with GA(4) and GA(3) at 1.9 A resolution. The overall structure of both complexes shows an alpha/beta-hydrolase fold similar to that of HSLs except for an amino-terminal lid. The GA-binding pocket corresponds to the substrate-binding site of HSLs. On the basis of the OsGID1 structure, we mutagenized important residues for GA binding and examined their binding activities. Almost all of them showed very little or no activity, confirming that the residues revealed by structural analysis are important for GA binding. The replacement of Ile 133 with Leu or Val-residues corresponding to those of the lycophyte Selaginella moellendorffii GID1s-caused an increase in the binding affinity for GA(34), a 2beta-hydroxylated GA(4). These observations indicate that GID1 originated from HSL and was further modified to have higher affinity and more strict selectivity for bioactive GAs by adapting the amino acids involved in GA binding in the course of plant evolution.

Molecular interactions of a soluble gibberellin receptor, GID1, with a rice DELLA protein, SLR1, and gibberellin.

GIBBERELLIN INSENSITIVE DWARF1 (GID1) encodes a soluble gibberellin (GA) receptor that shares sequence similarity with a hormone-sensitive lipase (HSL). Previously, a yeast two-hybrid (Y2H) assay revealed that the GID1-GA complex directly interacts with SLENDER RICE1 (SLR1), a DELLA repressor protein in GA signaling. Here, we demonstrated, by pull-down and bimolecular fluorescence complementation (BiFC) experiments, that the GA-dependent GID1-SLR1 interaction also occurs in planta. GA(4) was found to have the highest affinity to GID1 in Y2H assays and is the most effective form of GA in planta. Domain analyses of SLR1 using Y2H, gel filtration, and BiFC methods revealed that the DELLA and TVHYNP domains of SLR1 are required for the GID1-SLR1 interaction. To identify the important regions of GID1 for GA and SLR1 interactions, we used many different mutant versions of GID1, such as the spontaneous mutant GID1s, N- and C-terminal truncated GID1s, and mutagenized GID1 proteins with conserved amino acids replaced with Ala. The amino acid residues important for SLR1 interaction completely overlapped the residues required for GA binding that were scattered throughout the GID1 molecule. When we plotted these residues on the GID1 structure predicted by analogy with HSL tertiary structure, many residues were located at regions corresponding to the substrate binding pocket and lid. Furthermore, the GA-GID1 interaction was stabilized by SLR1. Based on these observations, we proposed a molecular model for interaction between GA, GID1, and SLR1.