BAC-FISH Based Physical Map of Endangered Catfish Clarias magur for Chromosome Cataloguing and Gene Isolation through Positional Cloning.

Construction of a physical chromosome map of a species is important for positional cloning, targeted marker development, fine mapping of genes, selection of candidate genes for molecular breeding, as well as understanding the genome organization. The genomic libraries in the form of bacterial artificial chromosome (BAC) clones are also a very useful resource for physical mapping and identification and isolation of full-length genes and the related cis acting elements. Some BAC-FISH based studies reported in the past were gene based physical chromosome maps of Clarias magur (magur) to understand the genome organization of the species and to establish the relationships with other species in respect to genes' organization and evolution in the past. In the present study, we generated end sequences of the BAC clones and analyzed those end sequences within the scaffolds of the draft genome of magur to identify and map the genes bioinformatically for each clone. A total of 36 clones mostly possessing genes were identified and used in probe construction and their subsequent hybridization on the metaphase chromosomes of magur. This study successfully mapped all 36 specific clones on 16 chromosome pairs, out of 25 pairs of magur chromosomes. These clones are now recognized as chromosome-specific makers, which are an aid in individual chromosome identification and fine assembly of the genome sequence, and will ultimately help in developing anchored genes' map on the chromosomes of C. magur for understanding their organization, inheritance of important fishery traits and evolution of magur with respect to channel catfish, zebrafish and other species.

Multi-ancestry fine mapping implicates OAS1 splicing in risk of severe COVID-19.

The OAS1/2/3 cluster has been identified as a risk locus for severe COVID-19 among individuals of European ancestry, with a protective haplotype of approximately 75 kilobases (kb) derived from Neanderthals in the chromosomal region 12q24.13. This haplotype contains a splice variant of OAS1, which occurs in people of African ancestry independently of gene flow from Neanderthals. Using trans-ancestry fine-mapping approaches in 20,779 hospitalized cases, we demonstrate that this splice variant is likely to be the SNP responsible for the association at this locus, thus strongly implicating OAS1 as an effector gene influencing COVID-19 severity.

A high-resolution HLA reference panel capturing global population diversity enables multi-ancestry fine-mapping in HIV host response.

Fine-mapping to plausible causal variation may be more effective in multi-ancestry cohorts, particularly in the MHC, which has population-specific structure. To enable such studies, we constructed a large (n = 21,546) HLA reference panel spanning five global populations based on whole-genome sequences. Despite population-specific long-range haplotypes, we demonstrated accurate imputation at G-group resolution (94.2%, 93.7%, 97.8% and 93.7% in admixed African (AA), East Asian (EAS), European (EUR) and Latino (LAT) populations). Applying HLA imputation to genome-wide association study data for HIV-1 viral load in three populations (EUR, AA and LAT), we obviated effects of previously reported associations from population-specific HIV studies and discovered a novel association at position 156 in HLA-B. We pinpointed the MHC association to three amino acid positions (97, 67 and 156) marking three consecutive pockets (C, B and D) within the HLA-B peptide-binding groove, explaining 12.9% of trait variance.

Mining of a novel esterase (est3S) gene from a cow rumen metagenomic library with organosphosphorus insecticides degrading capability: Catalytic insights by site directed mutations, docking, and molecular dynamic simulations.

A novel esterase (est3S) gene, 1026 bp in size, was cloned from a metagenomic library made of uncultured microorganisms from the contents of cow rumen. The esterolytic enzyme (Est3S) is composed of 342 amino acids and shows the highest identity with EstGK1 (71.7%) and EstZ3 (63.78%) esterases from the uncultured bacterium. The Est3S did not cluster in any up-to-date classes (I to XVIII) of esterase and lipase. Est3S protein molecular weight was determined to be 38 kDa by gel electrophoresis and showed optimum activity at pH 7.0 and 40 °C and is partially resistant to organic solvents. Est3S activity was enhanced by K+, Na+, Mg2+, and Ca2+ and its highest activity was observed toward the short-chain p-nitrophenyl esters. Additionally, Est3S can degrade chlorpyrifos (CP) and methyl parathion (70% to 80%) in an hour. A mutated Est3S (Ser132-Ala132) did not show any activity toward CP and ester substrates. Notably, the GHS132QG motif is superimposed with the homolog esterase and cutinase-like esterase. Therefore, Ser132 is the critical amino acid like other esterases. The Est3S is relatively stable with ester compounds, and the methyl parathion complex was confirmed by molecular dynamics simulation. NOVELTY STATEMENT: A novel esterase gene (est3S) expressing esters and organophosphorus insecticide degradation traits was isolated from the uncultured bacterium in the contents of cow rumen. The Est3S protein did not cluster in any up-to-date classes (I to XVIII) of esterase/lipase proteins. Est3S was stable with the ligands up to 100 ns during the molecular dynamic simulations.

Multi-omics integration analysis identifies novel genes for alcoholism with potential overlap with neurodegenerative diseases.

Identification of causal variants and genes underlying genome-wide association study (GWAS) loci is essential to understand the biology of alcohol use disorder (AUD) and drinks per week (DPW). Multi-omics integration approaches have shown potential for fine mapping complex loci to obtain biological insights to disease mechanisms. In this study, we use multi-omics approaches, to fine-map AUD and DPW associations at single SNP resolution to demonstrate that rs56030824 on chromosome 11 significantly reduces SPI1 mRNA expression in myeloid cells and lowers risk for AUD and DPW. Our analysis also identifies MAPT as a candidate causal gene specifically associated with DPW. Genes prioritized in this study show overlap with causal genes associated with neurodegenerative disorders. Multi-omics integration analyses highlight, genetic similarities and differences between alcohol intake and disordered drinking, suggesting molecular heterogeneity that might inform future targeted functional and cross-species studies.

Tightly linked Rps12 and Rps13 genes provide broad-spectrum Phytophthora resistance in soybean.

The Phytophtora root and stem rot is a serious disease in soybean. It is caused by the oomycete pathogen Phytophthora sojae. Growing Phytophthora resistant cultivars is the major method of controlling this disease. Resistance is race- or gene-specific; a single gene confers immunity against only a subset of the P. sojae isolates. Unfortunately, rapid evolution of new Phytophthora sojae virulent pathotypes limits the effectiveness of an Rps ("resistance to Phytophthora sojae") gene to 8-15 years. The current study was designed to investigate the effectiveness of Rps12 against a set of P. sojae isolates using recombinant inbred lines (RILs) that contain recombination break points in the Rps12 region. Our study revealed a unique Rps gene linked to the Rps12 locus. We named this novel gene as Rps13 that confers resistance against P. sojae isolate V13, which is virulent to recombinants that contains Rps12 but lack Rps13. The genetic distance between the two Rps genes is 4 cM. Our study revealed that two tightly linked functional Rps genes with distinct race-specificity provide broad-spectrum resistance in soybean. We report here the molecular markers for incorporating the broad-spectrum Phytophthora resistance conferred by the two Rps genes in commercial soybean cultivars.

INFIMA leverages multi-omics model organism data to identify effector genes of human GWAS variants.

Genome-wide association studies reveal many non-coding variants associated with complex traits. However, model organism studies largely remain as an untapped resource for unveiling the effector genes of non-coding variants. We develop INFIMA, Integrative Fine-Mapping, to pinpoint causal SNPs for diversity outbred (DO) mice eQTL by integrating founder mice multi-omics data including ATAC-seq, RNA-seq, footprinting, and in silico mutation analysis. We demonstrate INFIMA's superior performance compared to alternatives with human and mouse chromatin conformation capture datasets. We apply INFIMA to identify novel effector genes for GWAS variants associated with diabetes. The results of the application are available at http://www.statlab.wisc.edu/shiny/INFIMA/ .

A single genetic locus controls both expression of DPEP1/CHMP1A and kidney disease development via ferroptosis.

Genome-wide association studies (GWAS) have identified loci for kidney disease, but the causal variants, genes, and pathways remain unknown. Here we identify two kidney disease genes Dipeptidase 1 (DPEP1) and Charged Multivesicular Body Protein 1 A (CHMP1A) via the triangulation of kidney function GWAS, human kidney expression, and methylation quantitative trait loci. Using single-cell chromatin accessibility and genome editing, we fine map the region that controls the expression of both genes. Mouse genetic models demonstrate the causal roles of both genes in kidney disease. Cellular studies indicate that both Dpep1 and Chmp1a are important regulators of a single pathway, ferroptosis and lead to kidney disease development via altering cellular iron trafficking.

Genetic Linkage and Physical Mapping for an Oyster Mushroom (Pleurotus cornucopiae) and Quantitative Trait Locus Analysis for Cap Color.

Oyster mushrooms are grown commercially worldwide, especially in many developing countries, for their easy cultivation and high biological efficiency. Pleurotus cornucopiae is one of the main oyster mushroom species because of its gastronomic value and nutraceutical properties. Cap color is an important trait, since consumers prefer dark mushrooms, which are now represented by only a small portion of the commercial varieties. Breeding efforts are required to improve quality-related traits to satisfy various demands of consumers. Here, we present a saturated genetic linkage map of P. cornucopiae constructed by using a segregating population of 122 monokaryons and 3,449 single nucleotide polymorphism (SNP) markers generated by the 2b-RAD approach. The map contains 11 linkage groups covering 961.6 centimorgans (cM), with an average marker spacing of 0.27 cM. The genome of P. cornucopiae was de novo sequenced, resulting in 425 scaffolds (>1,000 bp) with a total genome size of 35.1 Mb. The scaffolds were assembled to the pseudochromosome level with the assistance of the genetic linkage map. A total of 97% SNP markers (3,357) were physically localized on 140 scaffolds that were assigned to 11 pseudochromosomes, with a total of 32.5 Mb, representing 92.5% of the whole genome. Six quantitative trait loci (QTL) controlling cap color of P. cornucopiae were detected, accounting for a total phenotypic variation of 65.6%, with the highest value for the QTL on pseudochromosome 5 (18%). The results of our study provide a solid base for marker-assisted breeding for agronomic traits and especially for studies on biological mechanisms controlling cap color in oyster mushrooms. IMPORTANCE Oyster mushrooms are produced and consumed all over the world. Pleurotus cornucopiae is one of the main oyster mushroom species. Dark-cap oyster mushrooms are becoming more and more popular with consumers, but dark varieties are rare on the market. Prerequisites for efficient breeding programs are the availability of high-quality whole genomes and genetic linkage maps. Genetic studies to fulfill some of these prerequisites have hardly been done for P. cornucopiae. In this study, we de novo sequenced the genome and constructed a saturated genetic linkage map for P. cornucopiae. The genetic linkage map was effectively used to assist the genome assembly and identify QTL that genetically control the trait cap color. As well, the genome characteristics of P. cornucopiae were compared to the closely related species Pleurotus ostreatus. The results provided a basis for understanding the genetic background and marker-assisted breeding of this economically important mushroom species.

qMrdd2, a novel quantitative resistance locus for maize rough dwarf disease.

BACKGROUND: Maize rough dwarf disease (MRDD), a widespread disease caused by four pathogenic viruses, severely reduces maize yield and grain quality. Resistance against MRDD is a complex trait that controlled by many quantitative trait loci (QTL) and easily influenced by environmental conditions. So far, many studies have reported numbers of resistant QTL, however, only one QTL have been cloned, so it is especially important to map and clone more genes that confer resistance to MRDD. RESULTS: In the study, a major quantitative trait locus (QTL) qMrdd2, which confers resistance to MRDD, was identified and fine mapped. qMrdd2, located on chromosome 2, was consistently identified in a 15-Mb interval between the simple sequence repeat (SSR) markers D184 and D1600 by using a recombinant inbred line (RIL) population derived from a cross between resistant ("80007") and susceptible ("80044") inbred lines. Using a recombinant-derived progeny test strategy, qMrdd2 was delineated to an interval of 577 kb flanked by markers N31 and N42. We further demonstrated that qMrdd2 is an incompletely dominant resistance locus for MRDD that reduced the disease severity index by 20.4%. CONCLUSIONS: A major resistance QTL (qMrdd2) have been identified and successfully refined into 577 kb region. This locus will be valuable for improving maize variety resistance to MRDD via marker-assisted selection (MAS).

The trans-ancestral genomic architecture of glycemic traits.

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Integrative genomic analysis of blood pressure and related phenotypes in rats.

Despite remarkable progress made in human genome-wide association studies, there remains a substantial gap between statistical evidence for genetic associations and functional comprehension of the underlying mechanisms governing these associations. As a means of bridging this gap, we performed genomic analysis of blood pressure (BP) and related phenotypes in spontaneously hypertensive rats (SHR) and their substrain, stroke-prone SHR (SHRSP), both of which are unique genetic models of severe hypertension and cardiovascular complications. By integrating whole-genome sequencing, transcriptome profiling, genome-wide linkage scans (maximum n=1415), fine congenic mapping (maximum n=8704), pharmacological intervention and comparative analysis with transcriptome-wide association study (TWAS) datasets, we searched causal genes and causal pathways for the tested traits. The overall results validated the polygenic architecture of elevated BP compared with a non-hypertensive control strain, Wistar Kyoto rats (WKY); e.g. inter-strain BP differences between SHRSP and WKY could be largely explained by an aggregate of BP changes in seven SHRSP-derived consomic strains. We identified 26 potential target genes, including rat homologs of human TWAS loci, for the tested traits. In this study, we re-discovered 18 genes that had previously been determined to contribute to hypertension or cardiovascular phenotypes. Notably, five of these genes belong to the kallikrein-kinin/renin-angiotensin systems (KKS/RAS), in which the most prominent differential expression between hypertensive and non-hypertensive alleles could be detected in rat Klk1 paralogs. In combination with a pharmacological intervention, we provide in vivo experimental evidence supporting the presence of key disease pathways, such as KKS/RAS, in a rat polygenic hypertension model.

Genome-wide characterization and expression profiling of the PDR gene family in tobacco (Nicotiana tabacum).

The pleiotropic drug resistance (PDR) proteins of the ATP-binding cassette (ABC) family play essential roles in physiological processes and have been characterized in many plant species. However, no comprehensive investigation of tobacco (Nicotiana tabacum), an important economic crop and a useful model plant for scientific research, has been presented. We identified 32 PDR genes in the tobacco genome and explored their domain organization, chromosomal distribution and evolution, promoter cis-elements, and expression profiles. A phylogenetic analysis revealed that tobacco has a significantly expanded number of PDR genes involved in plant defense. It also revealed that two tobacco PDR proteins may function as strigolactone transporters to regulate shoot branching, and several NtPDR genes may be involved in cadmium transport. Moreover, tissue expression profiles of NtPDR genes and their responses to several hormones and abiotic stresses were assessed using quantitative real-time PCR. Most of the NtPDR genes were regulated by jasmonate or salicylic acid, suggesting the important regulatory roles of NtPDRs in plant defense and secondary metabolism. They were also responsive to abiotic stresses, like drought and cold, and there was a strong correlation between the presence of promoter cis-elements and abiotic/biotic stress responses. These results provide useful clues for further in-depth studies on the functions of the tobacco PDR genes.

A Modified Actin (Gly65Val Substitution) Expressed in Cotton Disrupts Polymerization of Actin Filaments Leading to the Phenotype of Ligon Lintless-1 (Li1) Mutant.

Dynamic remodeling of the actin cytoskeleton plays a central role in the elongation of cotton fibers, which are the most important natural fibers in the global textile industry. Here, a high-resolution mapping approach combined with comparative sequencing and a transgenic method revealed that a G65V substitution in the cotton actin Gh_D04G0865 (GhACT17D in the wild-type) is responsible for the Gossypium hirsutum Ligon lintless-1 (Li1) mutant (GhACT17DM). In the mutant GhACT17DM from Li1 plant, Gly65 is substituted with valine on the lip of the nucleotide-binding domain of GhACT17D, which probably affects the polymerization of F-actin. Over-expression of GhACT17DM, but not GhACT17D, driven by either a CaMV35 promoter or a fiber-specific promoter in cotton produced a Li1-like phenotype. Compared with the wild-type control, actin filaments in Li1 fibers showed higher growth and shrinkage rates, decreased filament skewness and parallelness, and increased filament density. Taken together, our results indicate that the incorporation of GhACT17DM during actin polymerization disrupts the establishment and dynamics of the actin cytoskeleton, resulting in defective fiber elongation and the overall dwarf and twisted phenotype of the Li1 mutant.

Fine Mapping and Identification of BnaC06.FtsH1, a Lethal Gene That Regulates the PSII Repair Cycle in Brassica napus.

Photosystem II (PSII) is an important component of the chloroplast. The PSII repair cycle is crucial for the relief of photoinhibition and may be advantageous when improving stress resistance and photosynthetic efficiency. Lethal genes are widely used in the efficiency detection and method improvement of gene editing. In the present study, we identified the naturally occurring lethal mutant 7-521Y with etiolated cotyledons in Brassica napus, controlled by double-recessive genes (named cyd1 and cyd2). By combining whole-genome resequencing and map-based cloning, CYD1 was fine-mapped to a 29 kb genomic region using 15,167 etiolated individuals. Through cosegregation analysis and functional verification of the transgene, BnaC06.FtsH1 was determined to be the target gene; it encodes an filamentation temperature sensitive protein H 1 (FtsH1) hydrolase that degrades damaged PSII D1 in Arabidopsis thaliana. The expression of BnaC06.FtsH1 was high in the cotyledons, leaves, and flowers of B. napus, and localized in the chloroplasts. In addition, the expression of EngA (upstream regulation gene of FtsH) increased and D1 decreased in 7-521Y. Double mutants of FtsH1 and FtsH5 were lethal in A. thaliana. Through phylogenetic analysis, the loss of FtsH5 was identified in Brassica, and the remaining FtsH1 was required for PSII repair cycle. CYD2 may be a homologous gene of FtsH1 on chromosome A07 of B. napus. Our study provides new insights into lethal mutants, the findings may help improve the efficiency of the PSII repair cycle and biomass accumulation in oilseed rape.

Ph2 encodes the mismatch repair protein MSH7-3D that inhibits wheat homoeologous recombination.

Meiotic recombination is a critical process for plant breeding, as it creates novel allele combinations that can be exploited for crop improvement. In wheat, a complex allohexaploid that has a diploid-like behaviour, meiotic recombination between homoeologous or alien chromosomes is suppressed through the action of several loci. Here, we report positional cloning of Pairing homoeologous 2 (Ph2) and functional validation of the wheat DNA mismatch repair protein MSH7-3D as a key inhibitor of homoeologous recombination, thus solving a half-century-old question. Similar to ph2 mutant phenotype, we show that mutating MSH7-3D induces a substantial increase in homoeologous recombination (up to 5.5 fold) in wheat-wild relative hybrids, which is also associated with a reduction in homologous recombination. These data reveal a role for MSH7-3D in meiotic stabilisation of allopolyploidy and provides an opportunity to improve wheat's genetic diversity through alien gene introgression, a major bottleneck facing crop improvement.

GWAS of thyroid stimulating hormone highlights pleiotropic effects and inverse association with thyroid cancer.

Thyroid stimulating hormone (TSH) is critical for normal development and metabolism. To better understand the genetic contribution to TSH levels, we conduct a GWAS meta-analysis at 22.4 million genetic markers in up to 119,715 individuals and identify 74 genome-wide significant loci for TSH, of which 28 are previously unreported. Functional experiments show that the thyroglobulin protein-altering variants P118L and G67S impact thyroglobulin secretion. Phenome-wide association analysis in the UK Biobank demonstrates the pleiotropic effects of TSH-associated variants and a polygenic score for higher TSH levels is associated with a reduced risk of thyroid cancer in the UK Biobank and three other independent studies. Two-sample Mendelian randomization using TSH index variants as instrumental variables suggests a protective effect of higher TSH levels (indicating lower thyroid function) on risk of thyroid cancer and goiter. Our findings highlight the pleiotropic effects of TSH-associated variants on thyroid function and growth of malignant and benign thyroid tumors.

Karyotypic Evolution of Sauropsid Vertebrates Illuminated by Optical and Physical Mapping of the Painted Turtle and Slider Turtle Genomes.

Recent sequencing and software enhancements have advanced our understanding of the evolution of genomic structure and function, especially addressing novel evolutionary biology questions. Yet fragmentary turtle genome assemblies remain a challenge to fully decipher the genetic architecture of adaptive evolution. Here, we use optical mapping to improve the contiguity of the painted turtle (Chrysemys picta) genome assembly and use de novo fluorescent in situ hybridization (FISH) of bacterial artificial chromosome (BAC) clones, BAC-FISH, to physically map the genomes of the painted and slider turtles (Trachemys scripta elegans). Optical mapping increased C. picta's N50 by ~242% compared to the previous assembly. Physical mapping permitted anchoring ~45% of the genome assembly, spanning 5544 genes (including 20 genes related to the sex determination network of turtles and vertebrates). BAC-FISH data revealed assembly errors in C. picta and T. s. elegans assemblies, highlighting the importance of molecular cytogenetic data to complement bioinformatic approaches. We also compared C. picta's anchored scaffolds to the genomes of other chelonians, chicken, lizards, and snake. Results revealed a mostly one-to-one correspondence between chromosomes of painted and slider turtles, and high homology among large syntenic blocks shared with other turtles and sauropsids. Yet, numerous chromosomal rearrangements were also evident across chelonians, between turtles and squamates, and between avian and non-avian reptiles.

Prospects of GWAS and predictive breeding for European winter wheat's grain protein content, grain starch content, and grain hardness.

Grain quality traits determine the classification of registered wheat (Triticum aestivum L.) varieties. Although environmental factors and crop management practices exert a considerable influence on wheat quality traits, a significant proportion of the variance is attributed to the genetic factors. To identify the underlying genetic factors of wheat quality parameters viz., grain protein content (GPC), grain starch content (GSC), and grain hardness (GH), we evaluated 372 diverse European wheat varieties in replicated field trials in up to eight environments. We observed that all of the investigated traits hold a wide and significant genetic variation, and a significant negative correlation exists between GPC and GSC plus grain yield. Our association analyses based on 26,694 high-quality single nucleotide polymorphic markers revealed a strong quantitative genetic nature of GPC and GSC with associations on groups 2, 3, and 6 chromosomes. The identification of known Puroindoline-b gene for GH provided a positive analytic proof for our studies. We report that a locus QGpc.ipk-6A controls both GPC and GSC with opposite allelic effects. Based on wheat's reference and pan-genome sequences, the physical characterization of two loci viz., QGpc.ipk-2B and QGpc.ipk-6A facilitated the identification of the candidate genes for GPC. Furthermore, by exploiting additive and epistatic interactions of loci, we evaluated the prospects of predictive breeding for the investigated traits that suggested its efficient use in the breeding programs.

3D mapping and accelerated super-resolution imaging of the human genome using in situ sequencing.

There is a need for methods that can image chromosomes with genome-wide coverage, as well as greater genomic and optical resolution. We introduce OligoFISSEQ, a suite of three methods that leverage fluorescence in situ sequencing (FISSEQ) of barcoded Oligopaint probes to enable the rapid visualization of many targeted genomic regions. Applying OligoFISSEQ to human diploid fibroblast cells, we show how four rounds of sequencing are sufficient to produce 3D maps of 36 genomic targets across six chromosomes in hundreds to thousands of cells, implying a potential to image thousands of targets in only five to eight rounds of sequencing. We also use OligoFISSEQ to trace chromosomes at finer resolution, following the path of the X chromosome through 46 regions, with separate studies showing compatibility of OligoFISSEQ with immunocytochemistry. Finally, we combined OligoFISSEQ with OligoSTORM, laying the foundation for accelerated single-molecule super-resolution imaging of large swaths of, if not entire, human genomes.