Search

Knowing the genes involved in quantitative traits provides an entry point to understanding the biological bases of behavior, but there are very few examples where the pathway from genetic locus to behavioral change is known. To explore the role of specific genes in fear behavior, we mapped three fear-related traits, tested fourteen genes at six quantitative trait loci (QTLs) by quantitative complementation, and identified six genes. Four genes, Lamp, Ptprd, Nptx2, and Sh3gl, have known roles in synapse function; the fifth, Psip1, was not previously implicated in behavior; and the sixth is a long non-coding RNA, 4933413L06Rik, of unknown function. Variation in transcriptome and epigenetic modalities occurred preferentially in excitatory neurons, suggesting that genetic variation is more permissible in excitatory than inhibitory neuronal circuits. Our results relieve a bottleneck in using genetic mapping of QTLs to uncover biology underlying behavior and prompt a reconsideration of expected relationships between genetic and functional variation.

BackgroundStructural variations (SVs) have significant impacts on complex phenotypes by rearranging large amounts of DNA sequence.ResultsWe present a comprehensive SV catalog based on the whole-genome sequence of 1060 pigs (Sus scrofa) representing 101 breeds, covering 9.6% of the pig genome. This catalog includes 42,487 deletions, 37,913 mobile element insertions, 3308 duplications, 1664 inversions, and 45,184 break ends. Estimates of breed ancestry and hybridization using genotyped SVs align well with those from single nucleotide polymorphisms. Geographically stratified deletions are observed, along with known duplications of the KIT gene, responsible for white coat color in European pigs. Additionally, we identify a recent SINE element insertion in MYO5A transcripts of European pigs, potentially influencing alternative splicing patterns and coat color alterations. Furthermore, a Yorkshire-specific copy number gain within ABCG2 is found, impacting chromatin interactions and gene expression across multiple tissues over a stretch of genomic region of ~200 kb. Preliminary investigations into SV's impact on gene expression and traits using the Pig Genotype-Tissue Expression (PigGTEx) data reveal SV associations with regulatory variants and gene-trait pairs. For instance, a 51-bp deletion is linked to the lead eQTL of the lipid metabolism regulating gene FADS3, whose expression in embryo may affect loin muscle area, as revealed by our transcriptome-wide association studies.ConclusionsThis SV catalog serves as a valuable resource for studying diversity, evolutionary history, and functional shaping of the pig genome by processes like domestication, trait-based breeding, and adaptive evolution.

Genetic dissection of a QTL from wild emmer wheat, QGpc.huj.uh-5B.2, introgressed into bread wheat, identified candidate genes associated with tolerance to nitrogen deficiency, and potentially useful for improving nitrogen-use efficiency. Nitrogen (N) is an important macronutrient critical to wheat growth and development; its deficiency is one of the main factors causing reductions in grain yield and quality. N availability is significantly affected by drought or flooding, that are dependent on additional factors including soil type or duration and severity of stress. In a previous study, we identified a high grain protein content QTL (QGpc.huj.uh-5B.2) derived from the 5B chromosome of wild emmer wheat, that showed a higher proportion of explained variation under water-stress conditions. We hypothesized that this QTL is associated with tolerance to N deficiency as a possible mechanism underlying the higher effect under stress. To validate this hypothesis, we introgressed the QTL into the elite bread wheat var. Ruta, and showed that under N-deficient field conditions the introgression IL99 had a 33% increase in GPC (p < 0.05) compared to the recipient parent. Furthermore, evaluation of IL99 response to severe N deficiency (10% N) for 14 days, applied using a semi-hydroponic system under controlled conditions, confirmed its tolerance to N deficiency. Fine-mapping of the QTL resulted in 26 homozygous near-isogenic lines (BC4F5) segregating to N-deficiency tolerance. The QTL was delimited from - 28.28 to - 1.29 Mb and included 13 candidate genes, most associated with N-stress response, N transport, and abiotic stress responses. These genes may improve N-use efficiency under severely N-deficient environments. Our study demonstrates the importance of WEW as a source of novel candidate genes for sustainable improvement in tolerance to N deficiency in wheat.

BACKGROUND: Traditional recombinant inbred lines (RILs) are generated from repeated self-fertilization or brother-sister mating from the F1 hybrid of two inbred parents. Compared with the F2 population, RILs cumulate more crossovers between loci and thus increase the number of recombinants, resulting in an increased resolution of genetic mapping. Since they are inbred to the isogenic stage, another consequence of the heterozygosity reduction is the increased genetic variance and thus the increased power of QTL detection. Self-fertilization is the primary form of developing RILs in plants. Brother-sister mating is another way to develop RILs but in small laboratory animals. To ensure that the RILs have at least 98% of homozygosity, we need about seven generations of self-fertilization or 20 generations of brother-sister mating. Prior to homozygosity, these lines are called pre-recombinant inbred lines (PRERIL). Phenotypic values of traits in PRERILs are often collected but not used in QTL mapping. To perform QTL mapping in PRERILs, we need the recombination fraction between two markers at generation t for t < 7 (selfing) or t < 20 (brother-sister mating) so that the genotypes of QTL flanked by the markers can be inferred. RESULTS: In this study, we developed formulas to calculate the recombination fractions of PRERILs at generation t in self-fertilization, brother-sister mating, and random mating. In contrast to existing works in this topic, we used computer code to construct the transition matrix to form the Markov chain of genotype array between consecutive generations, the so-called recurrent equations. CONCLUSIONS: We provide R functions to calculate the recombination fraction using the newly developed recurrent equations of ordered genotype array. With the recurrent equations and the R code, users can perform QTL mapping in PRERILs. Substantial time and effort can be saved compared with QTL mapping in RILs.

Genome sequencing has revealed an incredible diversity of bacteria and archaea, but there are no fast and convenient tools for browsing across these genomes. It is cumbersome to view the prevalence of homologs for a protein of interest, or the gene neighborhoods of those homologs, across the diversity of the prokaryotes. We developed a web-based tool, fast.genomics, that uses two strategies to support fast browsing across the diversity of prokaryotes. First, the database of genomes is split up. The main database contains one representative from each of the 6,377 genera that have a high-quality genome, and additional databases for each taxonomic order contain up to 10 representatives of each species. Second, homologs of proteins of interest are identified quickly by using accelerated searches, usually in a few seconds. Once homologs are identified, fast.genomics can quickly show their prevalence across taxa, view their neighboring genes, or compare the prevalence of two different proteins. Fast.genomics is available at https://fast.genomics.lbl.gov.

High-throughput sequencing-based methods for bulked segregant analysis (BSA) allow for the rapid identification of genetic markers associated with traits of interest. BSA studies have successfully identified qualitative (binary) and quantitative trait loci (QTLs) using QTL mapping. However, most require population structures that fit the models available and a reference genome. Instead, high-throughput short-read sequencing can be combined with BSA of k-mers (BSA-k-mer) to map traits that appear refractory to standard approaches. This method can be applied to any organism and is particularly useful for species with genomes diverged from the closest sequenced genome. It is also instrumental when dealing with highly heterozygous and potentially polyploid genomes without phased haplotype assemblies and for which a single haplotype can control a trait. Finally, it is flexible in terms of population structure. Here, we apply the BSA-k-mer method for the rapid identification of candidate regions related to seed spot and seed size in diploid potato. Using a mixture of F1 and F2 individuals from a cross between 2 highly heterozygous parents, candidate sequences were identified for each trait using the BSA-k-mer approach. Using parental reads, we were able to determine the parental origin of the loci. Finally, we mapped the identified k-mers to a closely related potato genome to validate the method and determine the genomic loci underlying these sequences. The location identified for the seed spot matches with previously identified loci associated with pigmentation in potato. The loci associated with seed size are novel. Both loci are relevant in future breeding toward true seeds in potato.

Seeds are essential for plant reproduction, survival, and dispersal. Germination ability and successful establishment of young seedlings strongly depend on seed quality and on environmental factors such as nutrient availability. In tomato (Solanum lycopersicum) and many other species, seed quality and seedling establishment characteristics are determined by genetic variation, as well as the maternal environment in which the seeds develop and mature. The genetic contribution to variation in seed and seedling quality traits and environmental responsiveness can be estimated at transcriptome level in the dry seed by mapping genomic loci that affect gene expression (expression QTLs) in contrasting maternal environments. In this study, we applied RNA-sequencing to construct a linkage map and measure gene expression of seeds of a tomato recombinant inbred line (RIL) population derived from a cross between S. lycopersicum (cv. Moneymaker) and S. pimpinellifolium (G1.1554). The seeds matured on plants cultivated under different nutritional environments, i.e., on high phosphorus or low nitrogen. The obtained single-nucleotide polymorphisms (SNPs) were subsequently used to construct a genetic map. We show how the genetic landscape of plasticity in gene regulation in dry seeds is affected by the maternal nutrient environment. The combined information on natural genetic variation mediating (variation in) responsiveness to the environment may contribute to knowledge-based breeding programs aiming to develop crop cultivars that are resilient to stressful environments.

Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.

The population sizes of different retinal cell types vary between different strains of mice, and that variation can be mapped to genomic loci in order to identify its polygenic origin. In some cases, controlling genes act independently, whereas in other instances, they exhibit epistasis. Here, we identify an epistatic interaction revealed through the mapping of quantitative trait loci from a panel of recombinant inbred strains of mice. The population of retinal horizontal cells exhibits a twofold variation in number, mapping to quantitative trait loci on chromosomes 3 and 13, where these loci are shown to interact epistatically. We identify a prospective genetic interaction underlying this, mediated by the bHLH transcription factor Neurog2, at the chromosome 3 locus, functioning to repress the LIM homeodomain transcription factor Isl1, at the chromosome 13 locus. Using single and double conditional knockout mice, we confirm the countervailing actions of each gene, and validate in vitro a crucial role for two single nucleotide polymorphisms in the 5'UTR of Isl1, one of which yields a novel E-box, mediating the repressive action of Neurog2.

The short arms of the human acrocentric chromosomes 13, 14, 15, 21 and 22 (SAACs) share large homologous regions, including ribosomal DNA repeats and extended segmental duplications1,2. Although the resolution of these regions in the first complete assembly of a human genome—the Telomere-to-Telomere Consortium’s CHM13 assembly (T2T-CHM13)—provided a model of their homology3, it remained unclear whether these patterns were ancestral or maintained by ongoing recombination exchange. Here we show that acrocentric chromosomes contain pseudo-homologous regions (PHRs) indicative of recombination between non-homologous sequences. Utilizing an all-to-all comparison of the human pangenome from the Human Pangenome Reference Consortium4 (HPRC), we find that contigs from all of the SAACs form a community. A variation graph5 constructed from centromere-spanning acrocentric contigs indicates the presence of regions in which most contigs appear nearly identical between heterologous acrocentric chromosomes in T2T-CHM13. Except on chromosome 15, we observe faster decay of linkage disequilibrium in the pseudo-homologous regions than in the corresponding short and long arms, indicating higher rates of recombination6,7. The pseudo-homologous regions include sequences that have previously been shown to lie at the breakpoint of Robertsonian translocations8, and their arrangement is compatible with crossover in inverted duplications on chromosomes 13, 14 and 21. The ubiquity of signals of recombination between heterologous acrocentric chromosomes seen in the HPRC draft pangenome suggests that these shared sequences form the basis for recurrent Robertsonian translocations, providing sequence and population-based confirmation of hypotheses first developed from cytogenetic studies 50 years ago9., Our work depends on the HPRC draft human pangenome resource established in the accompanying Article4, and we thank the production and assembly groups for their efforts in establishing this resource. This work used the computational resources of the UTHSC Octopus cluster and NIH HPC Biowulf cluster. We acknowledge support in maintaining these systems that was critical to our analyses. The authors thank M. Miller for the development of a graphical synopsis of our study (Fig. 5); and R. Williams and N. Soranzo for support and guidance in the design and discussion of our work. This work was supported, in part, by National Institutes of Health/NIDA U01DA047638 (E.G.), National Institutes of Health/NIGMS R01GM123489 (E.G.), NSF PPoSS Award no. 2118709 (E.G. and C.F.), the Tennessee Governor’s Chairs programme (C.F. and E.G.), National Institutes of Health/NCI R01CA266339 (T.P., L.G.d.L. and J.L.G.), and the Intramural Research Program of the National Human Genome Research Institute, National Institutes of Health (A.R., S.K. and A.M.P.). We acknowledge support from Human Technopole (A.G.), Consiglio Nazionale delle Ricerche, Italy (S.B. and V.C.), and Stowers Institute for Medical Research (T.P., L.G.d.L., B.R. and J.L.G.)., Peer Reviewed, "Article signat per 13 autors/es: Andrea Guarracino, Silvia Buonaiuto, Leonardo Gomes de Lima, Tamara Potapova, Arang Rhie, Sergey Koren, Boris Rubinstein, Christian Fischer, Human Pangenome Reference Consortium, Jennifer L. Gerton, Adam M. Phillippy, Vincenza Colonna & Erik Garrison " Human Pangenome Reference Consortium: "Haley J. Abel, Lucinda L. Antonacci-Fulton, Mobin Asri, Gunjan Baid, Carl A. Baker, Anastasiya Belyaeva, Konstantinos Billis, Guillaume Bourque, Silvia Buonaiuto, Andrew Carroll, Mark J. P. Chaisson, Pi-Chuan Chang, Xian H. Chang, Haoyu Cheng, Justin Chu, Sarah Cody, Vincenza Colonna, Daniel E. Cook, Robert M. Cook-Deegan, Omar E. Cornejo, Mark Diekhans, Daniel Doerr, Peter Ebert, Jana Ebler, Evan E. Eichler, Jordan M. Eizenga, Susan Fairley, Olivier Fedrigo, Adam L. Felsenfeld, Xiaowen Feng, Christian Fischer, Paul Flicek, Giulio Formenti, Adam Frankish, Robert S. Fulton, Yan Gao, Shilpa Garg, Erik Garrison, Nanibaa’ A. Garrison, Carlos Garcia Giron, Richard E. Green, Cristian Groza, Andrea Guarracino, Leanne Haggerty, Ira Hall, William T. Harvey, Marina Haukness, David Haussler, Simon Heumos, Glenn Hickey, Kendra Hoekzema, Thibaut Hourlier, Kerstin Howe, Miten Jain, Erich D. Jarvis, Hanlee P. Ji, Eimear E. Kenny, Barbara A. Koenig, Alexey Kolesnikov, Jan O. Korbel, Jennifer Kordosky, Sergey Koren, HoJoon Lee, Alexandra P. Lewis, Heng Li, Wen-Wei Liao, Shuangjia Lu, Tsung-Yu Lu, Julian K. Lucas, Hugo Magalhães, Santiago Marco-Sola, Pierre Marijon, Charles Markello, Tobias Marschall, Fergal J. Martin, Ann McCartney, Jennifer McDaniel, Karen H. Miga, Matthew W. Mitchell, Jean Monlong, Jacquelyn Mountcastle, Katherine M. Munson, Moses Njagi Mwaniki, Maria Nattestad, Adam M. Novak, Sergey Nurk, Hugh E. Olsen, Nathan D. Olson, Benedict Paten, Trevor Pesout, Adam M. Phillippy, Alice B. Popejoy, David Porubsky, Pjotr Prins, Daniela Puiu, Mikko Rautiainen, Allison A. Regier, Arang Rhie, Samuel Sacco, Ashley D. Sanders, Valerie A. Schneider, Baergen I. S, Postprint (published version)

Bidirectional flow of information shapes the outcome of the host-pathogen interactions and depends on the genetics of each organism. Recent work has begun to use co-transcriptomic studies to shed light on this bidirectional flow, but it is unclear how plastic the co-transcriptome is in response to genetic variation in both the host and pathogen. To study co-transcriptome plasticity, we conducted transcriptomics using natural genetic variation in the pathogen, Botrytis cinerea, and large-effect genetic variation abolishing defense signaling pathways within the host, Arabidopsis thaliana. We show that genetic variation in the pathogen has a greater influence on the co-transcriptome than mutations that abolish defense signaling pathways in the host. Genome-wide association mapping using the pathogens' genetic variation and both organisms' transcriptomes allowed an assessment of how the pathogen modulates plasticity in response to the host. This showed that the differences in both organism's responses were linked to trans-expression quantitative trait loci (eQTL) hotspots within the pathogen's genome. These hotspots control gene sets in either the host or pathogen and show differential allele sensitivity to the host's genetic variation rather than qualitative host specificity. Interestingly, nearly all the trans-eQTL hotspots were unique to the host or pathogen transcriptomes. In this system of differential plasticity, the pathogen mediates the shift in the co-transcriptome more than the host.

Large-effect loci-those statistically significant loci discovered by genome-wide association studies or linkage mapping-associated with key traits segregate amidst a background of minor, often undetectable, genetic effects in wild and domesticated plants and animals. Accurately attributing mean differences and variance explained to the correct components in the linear mixed model analysis is vital for selecting superior progeny and parents in plant and animal breeding, gene therapy, and medical genetics in humans. Marker-assisted prediction and its successor, genomic prediction, have many advantages for selecting superior individuals and understanding disease risk. However, these two approaches are less often integrated to study complex traits with different genetic architectures. This simulation study demonstrates that the average semivariance can be applied to models incorporating Mendelian, oligogenic, and polygenic terms simultaneously and yields accurate estimates of the variance explained for all relevant variables. Our previous research focused on large-effect loci and polygenic variance separately. This work aims to synthesize and expand the average semivariance framework to various genetic architectures and the corresponding mixed models. This framework independently accounts for the effects of large-effect loci and the polygenic genetic background and is universally applicable to genetics studies in humans, plants, animals, and microbes.

The spiral gingers (Costus L.) are a pantropical genus of herbaceous perennial monocots; the Neotropical clade of Costus radiated rapidly in the past few million years into over 60 species. The Neotropical spiral gingers have a rich history of evolutionary and ecological research that can motivate and inform modern genetic investigations. Here, we present the first 2 chromosome-level genome assemblies in the genus, for C. pulverulentus and C. lasius, and briefly compare their synteny. We assembled the C. pulverulentus genome from a combination of short-read data, Chicago and Dovetail Hi-C chromatin-proximity sequencing, and alignment with a linkage map. We annotated the genome by mapping a C. pulverulentus transcriptome and querying mapped transcripts against a protein database. We assembled the C. lasius genome with Pacific Biosciences HiFi long reads and alignment to the C. pulverulentus genome. These 2 assemblies are the first published genomes for non-cultivated tropical plants. These genomes solidify the spiral gingers as a model system and will facilitate research on the poorly understood genetic basis of tropical plant diversification.

Key messageThe diploid wheat recessive stem rust resistance gene SrTm4 was fine-mapped to a 754-kb region on chromosome arm 2AmL and potential candidate genes were identified. Race Ug99 of Puccinia graminis f. sp. tritici (Pgt), the causal agent of wheat stem (or black) rust is one of the most serious threats to global wheat production. The identification, mapping, and deployment of effective stem rust resistance (Sr) genes are critical to reduce this threat. In this study, we generated SrTm4 monogenic lines and found that this gene confers resistance to North American and Chinese Pgt races. Using a large mapping population (9522 gametes), we mapped SrTm4 within a 0.06 cM interval flanked by marker loci CS4211 and 130K1519, which corresponds to a 1.0-Mb region in the Chinese Spring reference genome v2.1. A physical map of the SrTm4 region was constructed with 11 overlapping BACs from the resistant Triticum monococcum PI 306540. Comparison of the 754-kb physical map with the genomic sequence of Chinese Spring and a discontinuous BAC sequence of DV92 revealed a 593-kb chromosomal inversion in PI 306540. Within the candidate region, we identified an L-type lectin-domain containing receptor kinase (LLK1), which was disrupted by the proximal inversion breakpoint, as a potential candidate gene. Two diagnostic dominant markers were developed to detect the inversion breakpoints. In a survey of T. monococcum accessions, we identified 10 domesticated T. monococcum subsp. monococcum genotypes, mainly from the Balkans, carrying the inversion and showing similar mesothetic resistant infection types against Pgt races. The high-density map and tightly linked molecular markers developed in this study are useful tools to accelerate the deployment of SrTm4-mediated resistance in wheat breeding programs.

Ependymoma is a tumor of the brain or spinal cord. The two most common and aggressive molecular groups of ependymoma are the supratentorial ZFTA-fusion associated and the posterior fossa ependymoma group A. In both groups, tumors occur mainly in young children and frequently recur after treatment. Although molecular mechanisms underlying these diseases have recently been uncovered, they remain difficult to target and innovative therapeutic approaches are urgently needed. Here, we use genome-wide chromosome conformation capture (Hi-C), complemented with CTCF and H3K27ac ChIP-seq, as well as gene expression and DNA methylation analysis in primary and relapsed ependymoma tumors, to identify chromosomal conformations and regulatory mechanisms associated with aberrant gene expression. In particular, we observe the formation of new topologically associating domains ('neo-TADs') caused by structural variants, group-specific 3D chromatin loops, and the replacement of CTCF insulators by DNA hyper-methylation. Through inhibition experiments, we validate that genes implicated by these 3D genome conformations are essential for the survival of patient-derived ependymoma models in a group-specific manner. Thus, this study extends our ability to reveal tumor-dependency genes by 3D genome conformations even in tumors that lack targetable genetic alterations.

BackgroundHigh sequence identity between segmental duplications (SDs) can facilitate copy number variants (CNVs) via non-allelic homologous recombination (NAHR). These CNVs are one of the fundamental causes of genomic disorders such as the 3q29 deletion syndrome (del3q29S). There are 21 protein-coding genes lost or gained as a result of such recurrent 1.6-Mbp deletions or duplications, respectively, in the 3q29 locus. While NAHR plays a role in CNV occurrence, the factors that increase the risk of NAHR at this particular locus are not well understood.MethodsWe employed an optical genome mapping technique to characterize the 3q29 locus in 161 unaffected individuals, 16 probands with del3q29S and their parents, and 2 probands with the 3q29 duplication syndrome (dup3q29S). Long-read sequencing-based haplotype resolved de novo assemblies from 44 unaffected individuals, and 1 trio was used for orthogonal validation of haplotypes and deletion breakpoints.ResultsIn total, we discovered 34 haplotypes, of which 19 were novel haplotypes. Among these 19 novel haplotypes, 18 were detected in unaffected individuals, while 1 novel haplotype was detected on the parent-of-origin chromosome of a proband with the del3q29S. Phased assemblies from 44 unaffected individuals enabled the orthogonal validation of 20 haplotypes. In 89% (16/18) of the probands, breakpoints were confined to paralogous copies of a 20-kbp segment within the 3q29 SDs. In one del3q29S proband, the breakpoint was confined to a 374-bp region using long-read sequencing. Furthermore, we categorized del3q29S cases into three classes and dup3q29S cases into two classes based on breakpoints. Finally, we found no evidence of inversions in parent-of-origin chromosomes.ConclusionsWe have generated the most comprehensive haplotype map for the 3q29 locus using unaffected individuals, probands with del3q29S or dup3q29S, and available parents, and also determined the deletion breakpoint to be within a 374-bp region in one prob

BACKGROUND: There is interest in lingering non-specific symptoms after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, referred to as Long coronavirus disease 2019 (Long COVID-19). It remains unknown whether the risk of Long COVID-19 is associated with pre-existing comorbidities or initial COVID-19 severity, including infections due to new Omicron lineages which predominated in 2023. OBJECTIVES: The aim of this case report was to characterize the clinical features of acute XBB.1.5 infection followed by Long COVID-19. METHODS: We followed a 73-year old female resident of Rio de Janeiro with laboratory-confirmed SARS-CoV-2 during acute infection and subsequent months. The SARS-CoV-2 lineage was determined by genome sequencing. FINDINGS: The participant denied comorbidities and had completed a two-dose vaccination schedule followed by two booster doses eight months prior to SARS-CoV-2 infection. Primary infection by viral lineage XBB.1.5. was clinically mild, but the participant subsequently reported persistent fatigue. MAIN CONCLUSIONS: This case demonstrates that Long COVID-19 may develop even after mild disease due to SARS-CoV-2 in fully vaccinated and boosted individuals without comorbidities. Continued monitoring of new SARS-CoV-2 lineages and associated clinical outcomes is warranted. Measures to prevent infection should continue to be implemented including development of new vaccines and antivirals effective against novel variants.

Allopolyploidy greatly expands the range of possible regulatory interactions among functionally redundant homoeologous genes. However, connection between the emerging regulatory complexity and expression and phenotypic diversity in polyploid crops remains elusive. Here, we use diverse wheat accessions to map expression quantitative trait loci (eQTL) and evaluate their effects on the population-scale variation in homoeolog expression dosage. The relative contribution of cis- and trans-eQTL to homoeolog expression variation is strongly affected by both selection and demographic events. Though trans-acting effects play major role in expression regulation, the expression dosage of homoeologs is largely influenced by cis-acting variants, which appear to be subjected to selection. The frequency and expression of homoeologous gene alleles showing strong expression dosage bias are predictive of variation in yield-related traits, and have likely been impacted by breeding for increased productivity. Our study highlights the importance of genomic variants affecting homoeolog expression dosage in shaping agronomic phenotypes and points at their potential utility for improving yield in polyploid crops.

The search for quantitative trait loci that explain complex traits such as yield and drought tolerance has been ongoing in all crops. Methods such as biparental quantitative trait loci mapping and genome-wide association studies each have their own advantages and limitations. Multiparent advanced generation intercross populations contain more recombination events and genetic diversity than biparental mapping populations and are better able to estimate effect sizes of rare alleles than association mapping populations. Here, we discuss the results of using a multiparent advanced generation intercross population of doubled haploid maize lines created from 16 diverse founders to perform quantitative trait loci mapping. We compare 3 models that assume bi-allelic, founder, and ancestral haplotype allelic states for quantitative trait loci. The 3 methods have differing power to detect quantitative trait loci for a variety of agronomic traits. Although the founder approach finds the most quantitative trait loci, all methods are able to find unique quantitative trait loci, suggesting that each model has advantages for traits with different genetic architectures. A closer look at a well-characterized flowering time quantitative trait loci, qDTA8, which contains vgt1, highlights the strengths and weaknesses of each method and suggests a potential epistatic interaction. Overall, our results reinforce the importance of considering different approaches to analyzing genotypic datasets, and shows the limitations of binary SNP data for identifying multiallelic quantitative trait loci.

Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp. parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern United States, Canada, and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time.

Structural variants (SVs) are implicated in the etiology of Mendelian diseases but have been systematically underascertained owing to sequencing technology limitations. Long-read sequencing enables comprehensive detection of SVs, but approaches for prioritization of candidate SVs are needed. Structural variant Annotation and analysis (SvAnna) assesses all classes of SVs and their intersection with transcripts and regulatory sequences, relating predicted effects on gene function with clinical phenotype data. SvAnna places 87% of deleterious SVs in the top ten ranks. The interpretable prioritizations offered by SvAnna will facilitate the widespread adoption of long-read sequencing in diagnostic genomics. SvAnna is available at https://github.com/TheJacksonLaboratory/SvAnn a .

The sexual strain of the planarian Schmidtea mediterranea, indigenous to Tunisia and several Mediterranean islands, is a hermaphrodite1,2. Here we isolate individual chromosomes and use sequencing, Hi-C3,4 and linkage mapping to assemble a chromosome-scale genome reference. The linkage map reveals an extremely low rate of recombination on chromosome 1. We confirm suppression of recombination on chromosome 1 by genotyping individual sperm cells and oocytes. We show that previously identified genomic regions that maintain heterozygosity even after prolonged inbreeding make up essentially all of chromosome 1. Genome sequencing of individuals isolated in the wild indicates that this phenomenon has evolved specifically in populations from Sardinia and Corsica. We find that most known master regulators5-13 of the reproductive system are located on chromosome 1. We used RNA interference14,15 to knock down a gene with haplotype-biased expression, which led to the formation of a more pronounced female mating organ. On the basis of these observations, we propose that chromosome 1 is a sex-primed autosome primed for evolution into a sex chromosome.

BackgroundEstimating and accounting for hidden variables is widely practiced as an important step in molecular quantitative trait locus (molecular QTL, henceforth "QTL") analysis for improving the power of QTL identification. However, few benchmark studies have been performed to evaluate the efficacy of the various methods developed for this purpose.ResultsHere we benchmark popular hidden variable inference methods including surrogate variable analysis (SVA), probabilistic estimation of expression residuals (PEER), and hidden covariates with prior (HCP) against principal component analysis (PCA)-a well-established dimension reduction and factor discovery method-via 362 synthetic and 110 real data sets. We show that PCA not only underlies the statistical methodology behind the popular methods but is also orders of magnitude faster, better-performing, and much easier to interpret and use.ConclusionsTo help researchers use PCA in their QTL analysis, we provide an R package PCAForQTL along with a detailed guide, both of which are freely available at https://github.com/heatherjzhou/PCAForQTL . We believe that using PCA rather than SVA, PEER, or HCP will substantially improve and simplify hidden variable inference in QTL mapping as well as increase the transparency and reproducibility of QTL research.

As one of the US Department of Agriculture-Agricultural Research Service flagship databases, GrainGenes (https://wheat.pw.usda.gov) serves the data and community needs of globally distributed small grains researchers for the genetic improvement of the Triticeae family and Avena species that include wheat, barley, rye and oat. GrainGenes accomplishes its mission by continually enriching its cross-linked data content following the findable, accessible, interoperable and reusable principles, enhancing and maintaining an intuitive web interface, creating tools to enable easy data access and establishing data connections within and between GrainGenes and other biological databases to facilitate knowledge discovery. GrainGenes operates within the biological database community, collaborates with curators and genome sequencing groups and contributes to the AgBioData Consortium and the International Wheat Initiative through the Wheat Information System (WheatIS). Interactive and linked content is paramount for successful biological databases and GrainGenes now has 2917 manually curated gene records, including 289 genes and 254 alleles from the Wheat Gene Catalogue (WGC). There are >4.8 million gene models in 51 genome browser assemblies, 6273 quantitative trait loci and >1.4 million genetic loci on 4756 genetic and physical maps contained within 443 mapping sets, complete with standardized metadata. Most notably, 50 new genome browsers that include outputs from the Wheat and Barley PanGenome projects have been created. We provide an example of an expression quantitative trait loci track on the International Wheat Genome Sequencing Consortium Chinese Spring wheat browser to demonstrate how genome browser tracks can be adapted for different data types. To help users benefit more from its data, GrainGenes created four tutorials available on YouTube. GrainGenes is executing its vision of service by continuously responding to the needs of the global small grains community b

R-loops play versatile roles in many physiological and pathological processes, and are of great interest to scientists in multiple fields. However, controversy about their genomic localization and incomplete understanding of their regulatory network raise great challenges for R-loop research. Here, we present R-loopBase (https://rloopbase.nju.edu.cn) to tackle these pressing issues by systematic integration of genomics and literature data. First, based on 107 high-quality genome-wide R-loop mapping datasets generated by 11 different technologies, we present a reference set of human R-loop zones for high-confidence R-loop localization, and spot conservative genomic features associated with R-loop formation. Second, through literature mining and multi-omics analyses, we curate the most comprehensive list of R-loop regulatory proteins and their targeted R-loops in multiple species to date. These efforts help reveal a global regulatory network of R-loop dynamics and its potential links to the development of cancers and neurological diseases. Finally, we integrate billions of functional genomic annotations, and develop interactive interfaces to search, visualize, download and analyze R-loops and R-loop regulators in a well-annotated genomic context. R-loopBase allows all users, including those with little bioinformatics background to utilize these data for their own research. We anticipate R-loopBase will become a one-stop resource for the R-loop community.

BackgroundPilocytic astrocytoma (PA) is the most common pediatric brain tumor. PA has at least a 50% higher incidence in populations of European ancestry compared to other ancestral groups, which may be due in part to genetic differences.MethodsWe first compared the global proportions of European, African, and Amerindian ancestries in 301 PA cases and 1185 controls of self-identified Latino ethnicity from the California Biobank. We then conducted admixture mapping analysis to assess PA risk with local ancestry.ResultsWe found PA cases had a significantly higher proportion of global European ancestry than controls (case median = 0.55, control median = 0.51, P value = 3.5x10-3). Admixture mapping identified 13 SNPs in the 6q14.3 region (SNX14) contributing to risk, as well as three other peaks approaching significance on chromosomes 7, 10 and 13. Downstream fine mapping in these regions revealed several SNPs potentially contributing to childhood PA risk.ConclusionsThere is a significant difference in genomic ancestry associated with Latino PA risk and several genomic loci potentially mediating this risk.

Existing human genome assemblies have almost entirely excluded repetitive sequences within and near centromeres, limiting our understanding of their organization, evolution, and functions, which include facilitating proper chromosome segregation. Now, a complete, telomere-to-telomere human genome assembly (T2T-CHM13) has enabled us to comprehensively characterize pericentromeric and centromeric repeats, which constitute 6.2% of the genome (189.9 megabases). Detailed maps of these regions revealed multimegabase structural rearrangements, including in active centromeric repeat arrays. Analysis of centromere-associated sequences uncovered a strong relationship between the position of the centromere and the evolution of the surrounding DNA through layered repeat expansions. Furthermore, comparisons of chromosome X centromeres across a diverse panel of individuals illuminated high degrees of structural, epigenetic, and sequence variation in these complex and rapidly evolving regions.

Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is responsible for significant yield losses worldwide, which can be minimized by the deployment of Pst resistance genes. Yr78 is an adult plant partial-resistance gene that has remained effective against the post-2000 virulent Pst races. In this study, we generated a high-resolution map of Yr78 based on 6,124 segregating chromosomes. We mapped Yr78 within a 0.05-cM interval on the short arm of chromosome 6B, which corresponds to an 11.16 Mb region between TraesCS6B02G116200 and TraesCS6B02G118000 in the 'Chinese Spring' Ref Seq. v1.1 genome. This interval is likely larger because it includes the unassembled NOR-B2 region, which may have contributed to the low recombination rate detected in this region. The Yr78 candidate region includes 15 genes that were prioritized for future functional studies based on their annotated function and polymorphisms between susceptible and resistant genotypes. Using exome capture data, we identified five major haplotypes in the candidate gene region, with the H1 haplotype associated with Yr78. The H1 haplotype was not detected in tetraploid wheat (Triticum turgidum L.) but was found in ∼30% of the common wheat cultivars (Triticum aestivum L.), suggesting that the associated resistance to stripe rust may have favored the selection of this haplotype. We developed two diagnostic molecular markers for the H1 haplotype that will facilitate the deployment of Yr78 in wheat breeding programs.

Improving our understanding of the genes regulating grain yield can contribute to the development of more productive wheat varieties. Previously, a highly significant QTL affecting spikelet number per spike (SNS), grain number per spike (GNS) and grain yield was detected on chromosome arm 7AL in multiple genome-wide association studies. Using a high-resolution genetic map, we established that the A-genome homeolog of WHEAT ORTHOLOG OF APO1 (WAPO-A1) was a leading candidate gene for this QTL. Using mutants and transgenic plants, we demonstrate in this study that WAPO-A1 is the causal gene underpinning this QTL. Loss-of-function mutants wapo-A1 and wapo-B1 showed reduced SNS in tetraploid wheat, and the effect was exacerbated in wapo1 combining both mutations. By contrast, spikes of transgenic wheat plants carrying extra copies of WAPO-A1 driven by its native promoter had higher SNS, a more compact spike apical region and a smaller terminal spikelet than the wild type. Taken together, these results indicate that WAPO1 affects SNS by regulating the timing of terminal spikelet formation. Both transgenic and wapo1 mutant plants showed a wide range of floral abnormalities, indicating additional roles of WAPO1 on wheat floral development. Previously, we found three widespread haplotypes in the QTL region (H1, H2 and H3), each associated with particular WAPO-A1 alleles. Results from this and our previous study show that the WAPO-A1 allele in the H1 haplotype (115-bp deletion in the promoter) is expressed at significantly lower levels in the developing spikes than the alleles in the H2 and H3 haplotypes, resulting in reduced SNS. Field experiments also showed that the H2 haplotype is associated with the strongest effects in increasing SNS and GNS (H2>H3>H1). The H2 haplotype is already present in most modern common wheat varieties but is rare in durum wheat, where it might be particularly useful to improve grain yield.

Key messageWe discovered a natural FT-A2 allele that increases grain number per spike in both pasta and bread wheat with limited effect on heading time. Increases in wheat grain yield are necessary to meet future global food demands. A previous study showed that loss-of-function mutations in FLOWERING LOCUS T2 (FT2) increase spikelet number per spike (SNS), an important grain yield component. However, these mutations were also associated with reduced fertility, offsetting the beneficial effect of the increases in SNS on grain number. Here, we report a natural mutation resulting in an aspartic acid to alanine change at position 10 (D10A) associated with significant increases in SNS and no negative effects on fertility. Using a high-density genetic map, we delimited the SNS candidate region to a 5.2-Mb region on chromosome 3AS including 28 genes. Among them, only FT-A2 showed a non-synonymous polymorphism (D10A) present in two different populations segregating for the SNS QTL on chromosome arm 3AS. These results, together with the known effect of the ft-A2 mutations on SNS, suggest that variation in FT-A2 is the most likely cause of the observed differences in SNS. We validated the positive effects of the A10 allele on SNS, grain number, and grain yield per spike in near-isogenic tetraploid wheat lines and in an hexaploid winter wheat population. The A10 allele is present at very low frequency in durum wheat and at much higher frequency in hexaploid wheat, particularly in winter and fall-planted spring varieties. These results suggest that the FT-A2 A10 allele may be particularly useful for improving grain yield in durum wheat and fall-planted common wheat varieties.

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

Wheat stem (or black) rust, caused by Puccinia graminis f. sp. tritici (Pgt), has been historically among the most devastating global fungal diseases of wheat. The recent occurrence and spread of new virulent races such as Ug99 have prompted global efforts to identify and isolate more effective stem rust resistance (Sr) genes. Here, we report the map-based cloning of the Ug99-effective SrTm5 gene from diploid wheat Triticum monococcum accession PI 306540 that encodes a typical coiled-coil nucleotide-binding leucine-rich repeat protein. This gene, designated as Sr22b, is a new allele of Sr22 with a rare insertion of a large (13.8-kb) retrotransposon into its second intron. Biolistic transformation of an ~112-kb circular bacterial artificial chromosome plasmid carrying Sr22b into the susceptible wheat variety Fielder was sufficient to confer resistance to stem rust. In a survey of 168 wheat genotypes, Sr22b was present only in cultivated T. monococcum subsp. monococcum accessions but absent in all tested tetraploid and hexaploid wheat lines. We developed a diagnostic molecular marker for Sr22b and successfully introgressed a T. monococcum chromosome segment containing this gene into hexaploid wheat to accelerate its deployment and pyramiding with other Sr genes in wheat breeding programmes. Sr22b can be a valuable component of gene pyramids or transgenic cassettes combining different resistance genes to control this devastating disease.

PurposeThere is limited payer coverage for genome sequencing (GS) relative to exome sequencing (ES) in the U.S. Our objective was to assess payers' considerations for coverage of GS versus coverage of ES and requirements payers have for coverage of GS. The study was conducted by the NIH-funded Clinical Sequencing Evidence-Generating Research Consortium (CSER).MethodsWe conducted semi-structured interviews with representatives of private payer organizations (payers, N = 12) on considerations and evidentiary and other needs for coverage of GS and ES. Data were analyzed using thematic analysis.ResultsWe described four categories of findings and solutions: demonstrated merits of GS versus ES, enhanced methods for evidence generation, consistent laboratory processes/sequencing methods, and enhanced implementation/care delivery. Payers see advantages to GS vs. ES and are open to broader GS coverage but need more proof of these advantages to consider them in coverage decision-making. Next steps include establishing evidence of benefits in specific clinical scenarios, developing quality standards, ensuring transparency of laboratory methods, developing clinical centers of excellence, and incorporating the role of genetic professionals.ConclusionBy comparing coverage considerations for GS and ES, we identified a path forward for coverage of GS. Future research should explicitly address payers' conditions for coverage.

Floral organs are produced by floral meristems (FMs), which harbor stem cells in their centers. Since each flower only has a finite number of organs, the stem cell activity of an FM will always terminate at a specific time point, a process termed floral meristem termination (FMT). Variation in the timing of FMT can give rise to floral morphological diversity, but how this process is fine-tuned at a developmental and evolutionary level is poorly understood. Flowers from the genus Aquilegia share identical floral organ arrangement except for stamen whorl number (SWN), making Aquilegia a well-suited system for investigation of this process: differences in SWN between species represent differences in the timing of FMT. By crossing A. canadensis and A. brevistyla, quantitative trait locus (QTL) mapping has revealed a complex genetic architecture with seven QTL. We explored potential candidate genes under each QTL and characterized novel expression patterns of select loci of interest using in situ hybridization. To our knowledge, this is the first attempt to dissect the genetic basis of how natural variation in the timing of FMT is regulated, and our results provide insight into how floral morphological diversity can be generated at the meristematic level.

The Stiff Stalk heterotic pool is a foundation of US maize seed parent germplasm and has been heavily utilized by both public and private maize breeders since its inception in the 1930s. Flowering time and plant height are critical characteristics for both inbred parents and their test crossed hybrid progeny. To study these traits, a 6-parent multiparent advanced generation intercross population was developed including maize inbred lines B73, B84, PHB47 (B37 type), LH145 (B14 type), PHJ40 (novel early Stiff Stalk), and NKH8431 (B73/B14 type). A set of 779 doubled haploid lines were evaluated for flowering time and plant height in 2 field replicates in 2016 and 2017, and a subset of 689 and 561 doubled haploid lines were crossed to 2 testers, respectively, and evaluated as hybrids in 2 locations in 2018 and 2019 using an incomplete block design. Markers were derived from a practical haplotype graph built from the founder whole genome assemblies and genotype-by-sequencing and exome capture-based sequencing of the population. Genetic mapping utilizing an update to R/qtl2 revealed differing profiles of significant loci for both traits between 635 of the DH lines and 2 sets of 570 and 471 derived hybrids. Genomic prediction was used to test the feasibility of predicting hybrid phenotypes based on the per se data. Predictive abilities were highest on direct models trained using the data they would predict (0.55-0.63), and indirect models trained using per se data to predict hybrid traits had slightly lower predictive abilities (0.49-0.55). Overall, this finding is consistent with the overlapping and nonoverlapping significant quantitative trait loci found within the per se and hybrid populations and suggests that selections for phenology traits can be made effectively on doubled haploid lines before hybrid data is available.

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Key messageWe propose the MADS-box transcription factor SVP-A1 as a promising candidate gene for the elongated glume locus P2, which maps to chromosome 6A instead of the previously proposed chromosome 7B. In rice and wheat, glume and floral organ length are positively correlated with grain size, making them an important target to increase grain size and potentially yield. The wheat subspecies Triticum ispahanicum is known to develop elongated glumes and floral organs as well as long grains. These multiple phenotypic effects are controlled by the P2 locus, which was previously mapped to wheat chromosome 7B. Using three mapping populations, we show that the long glume locus P2 does not map to chromosome 7B, but instead maps to a 1.68 Mbp interval on chromosome 6A. Within this interval, we identified SVP-A1, a MADS box transcription factor which is the direct ortholog of the maize gene underlying the 'pod corn' Tunicate locus and is a paralog to the T. polonicum elongated glume P1 gene. In T. ispahanicum, we identified a unique allele which has a 482-bp deletion in the SVP-A1 promoter and is associated with ectopic and higher expression of SVP-A1 in the elongated glumes and floral organs. We used near-isogenic lines (NILs) to show that P2 has a consistent positive effect on the length of glume, lemma, palea, spike and grain. Based on the mapping data, natural variation, biological function of SVP genes in cereals and expression analyses, we propose the MADS-box transcription factor SVP-A1 as a promising candidate for P2.

The freshwater snail Biomphalaria glabrata is an important intermediate host of the parasite Schistosoma mansoni that causes human intestinal schistosomiasis. To better understand vector snail biology and help advance innovative snail control strategies, we have developed a new snail model consisting of two homozygous B. glabrata lines (iM line and iBS90) with sharply contrasting schistosome-resistance phenotypes. We produced and compared high-quality genome sequences for iM line and iBS90 which were assembled from 255 (N50 = 22.7 Mb) and 346 (N50 = 19.4 Mb) scaffolds, respectively. Using F2 offspring bred from the two lines and the newly generated iM line genome, we constructed 18 linkage groups (representing the 18 haploid chromosomes) covering 96% of the genome and identified three new QTLs (quantitative trait loci), two involved in snail resistance/susceptibility and one relating to body pigmentation. This study provides excellent genomic resources for unveiling complex vector snail biology, reveals genomic difference between resistant and susceptible lines, and offers novel insights into genetic mechanism of the compatibility between snail and schistosome.

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

Trichechus manatus and Trichechus inunguis are the two Sirenia species that occur in the Americas. Despite their increasing extinction risk, many aspects of their biology remain understudied, including the repetitive DNA fraction of their genomes. Here we used the sequenced genome of T. manatus and TAREAN to identify satellite DNAs (satDNAs) in this species. We report the first description of TMAsat, a satDNA comprising ~0.87% of the genome, with ~684bp monomers and centromeric localization. In T. inunguis, TMAsat showed similar monomer length, chromosome localization and conserved CENP-B box-like motifs as in T. manatus. We also detected this satDNA in the Dugong dugon and in the now extinct Hydrodamalis gigas genomes. The neighbor-joining tree shows that TMAsat sequences from T. manatus, T. inunguis, D. dugon, and H. gigas lack species-specific clusters, which disagrees with the predictions of concerted evolution. We detected a divergent TMAsat-like homologous sequence in elephants and hyraxes, but not in other mammals, suggesting this sequence was already present in the common ancestor of Paenungulata, and later became a satDNA in the Sirenians. This is the first description of a centromeric satDNA in manatees and will facilitate the inclusion of Sirenia in future studies of centromeres and satDNA biology.

Platelet aggregation at the site of atherosclerotic vascular injury is the underlying pathophysiology of myocardial infarction and stroke. To build upon prior GWAS, here we report on 16 loci identified through a whole genome sequencing (WGS) approach in 3,855 NHLBI Trans-Omics for Precision Medicine (TOPMed) participants deeply phenotyped for platelet aggregation. We identify the RGS18 locus, which encodes a myeloerythroid lineage-specific regulator of G-protein signaling that co-localizes with expression quantitative trait loci (eQTL) signatures for RGS18 expression in platelets. Gene-based approaches implicate the SVEP1 gene, a known contributor of coronary artery disease risk. Sentinel variants at RGS18 and PEAR1 are associated with thrombosis risk and increased gastrointestinal bleeding risk, respectively. Our WGS findings add to previously identified GWAS loci, provide insights regarding the mechanism(s) by which genetics may influence cardiovascular disease risk, and underscore the importance of rare variant and regulatory approaches to identifying loci contributing to complex phenotypes.

Lima bean (Phaseolus lunatus L.), one of the five domesticated Phaseolus bean crops, shows a wide range of ecological adaptations along its distribution range from Mexico to Argentina. These adaptations make it a promising crop for improving food security under predicted scenarios of climate change in Latin America and elsewhere. In this work, we combine long and short read sequencing technologies with a dense genetic map from a biparental population to obtain the chromosome-level genome assembly for Lima bean. Annotation of 28,326 gene models show high diversity among 1917 genes with conserved domains related to disease resistance. Structural comparison across 22,180 orthologs with common bean reveals high genome synteny and five large intrachromosomal rearrangements. Population genomic analyses show that wild Lima bean is organized into six clusters with mostly non-overlapping distributions and that Mesomerican landraces can be further subdivided into three subclusters. RNA-seq data reveal 4275 differentially expressed genes, which can be related to pod dehiscence and seed development. We expect the resources presented here to serve as a solid basis to achieve a comprehensive view of the degree of convergent evolution of Phaseolus species under domestication and provide tools and information for breeding for climate change resiliency.

Corymbia citriodora is a member of the predominantly Southern Hemisphere Myrtaceae family, which includes the eucalypts (Eucalyptus, Corymbia and Angophora; ~800 species). Corymbia is grown for timber, pulp and paper, and essential oils in Australia, South Africa, Asia, and Brazil, maintaining a high-growth rate under marginal conditions due to drought, poor-quality soil, and biotic stresses. To dissect the genetic basis of these desirable traits, we sequenced and assembled the 408 Mb genome of Corymbia citriodora, anchored into eleven chromosomes. Comparative analysis with Eucalyptus grandis reveals high synteny, although the two diverged approximately 60 million years ago and have different genome sizes (408 vs 641 Mb), with few large intra-chromosomal rearrangements. C. citriodora shares an ancient whole-genome duplication event with E. grandis but has undergone tandem gene family expansions related to terpene biosynthesis, innate pathogen resistance, and leaf wax formation, enabling their successful adaptation to biotic/abiotic stresses and arid conditions of the Australian continent.

BackgroundDespite marked recent improvements in long-read sequencing technology, the assembly of diploid genomes remains a difficult task. A major obstacle is distinguishing between alternative contigs that represent highly heterozygous regions. If primary and secondary contigs are not properly identified, the primary assembly will overrepresent both the size and complexity of the genome, which complicates downstream analysis such as scaffolding.ResultsHere we illustrate a new method, which we call HapSolo, that identifies secondary contigs and defines a primary assembly based on multiple pairwise contig alignment metrics. HapSolo evaluates candidate primary assemblies using BUSCO scores and then distinguishes among candidate assemblies using a cost function. The cost function can be defined by the user but by default considers the number of missing, duplicated and single BUSCO genes within the assembly. HapSolo performs hill climbing to minimize cost over thousands of candidate assemblies. We illustrate the performance of HapSolo on genome data from three species: the Chardonnay grape (Vitis vinifera), with a genome of 490 Mb, a mosquito (Anopheles funestus; 200 Mb) and the Thorny Skate (Amblyraja radiata; 2650 Mb).ConclusionsHapSolo rapidly identified candidate assemblies that yield improvements in assembly metrics, including decreased genome size and improved N50 scores. Contig N50 scores improved by 35%, 9% and 9% for Chardonnay, mosquito and the thorny skate, respectively, relative to unreduced primary assemblies. The benefits of HapSolo were amplified by down-stream analyses, which we illustrated by scaffolding with Hi-C data. We found, for example, that prior to the application of HapSolo, only 52% of the Chardonnay genome was captured in the largest 19 scaffolds, corresponding to the number of chromosomes. After the application of HapSolo, this value increased to ~ 84%. The improvements for the mosquito's largest three scaffolds, representing

Variability in SARS-CoV-2 susceptibility and COVID-19 disease severity between individuals is partly due to genetic factors. Here, we identify 4 genomic loci with suggestive associations for SARS-CoV-2 susceptibility and 19 for COVID-19 disease severity. Four of these 23 loci likely have an ethnicity-specific component. Genome-wide association study (GWAS) signals in 11 loci colocalize with expression quantitative trait loci (eQTLs) associated with the expression of 20 genes in 62 tissues/cell types (range: 1:43 tissues/gene), including lung, brain, heart, muscle, and skin as well as the digestive system and immune system. We perform genetic fine mapping to compute 99% credible SNP sets, which identify 10 GWAS loci that have eight or fewer SNPs in the credible set, including three loci with one single likely causal SNP. Our study suggests that the diverse symptoms and disease severity of COVID-19 observed between individuals is associated with variants across the genome, affecting gene expression levels in a wide variety of tissue types.

We report the largest and most diverse genetic study of type 1 diabetes (T1D) to date (61,427 participants), yielding 78 genome-wide-significant (P < 5 × 10-8) regions, including 36 that are new. We define credible sets of T1D-associated variants and show that they are enriched in immune-cell accessible chromatin, particularly CD4+ effector T cells. Using chromatin-accessibility profiling of CD4+ T cells from 115 individuals, we map chromatin-accessibility quantitative trait loci and identify five regions where T1D risk variants co-localize with chromatin-accessibility quantitative trait loci. We highlight rs72928038 in BACH2 as a candidate causal T1D variant leading to decreased enhancer accessibility and BACH2 expression in T cells. Finally, we prioritize potential drug targets by integrating genetic evidence, functional genomic maps and immune protein-protein interactions, identifying 12 genes implicated in T1D that have been targeted in clinical trials for autoimmune diseases. These findings provide an expanded genomic landscape for T1D.

Bison are an icon of the American West and an ecologically, commercially, and culturally important species. Despite numbering in the hundreds of thousands today, conservation concerns remain for the species, including the impact on genetic diversity of a severe bottleneck around the turn of the 20th century and genetic introgression from domestic cattle. Genetic diversity and admixture are best evaluated at genome-wide scale, for which a high-quality reference is necessary. Here, we use trio binning of long reads from a bison-Simmental cattle (Bos taurus taurus) male F1 hybrid to sequence and assemble the genome of the American plains bison (Bison bison bison). The male haplotype genome is chromosome-scale, with a total length of 2.65 Gb across 775 scaffolds (839 contigs) and a scaffold N50 of 87.8 Mb. Our bison genome is ~13× more contiguous overall and ~3400× more contiguous at the contig level than the current bison reference genome. The bison genome sequence presented here (ARS-UCSC_bison1.0) will enable new research into the evolutionary history of this iconic megafauna species and provide a new tool for the management of bison populations in federal and commercial herds.

Sex chromosomes are classically predicted to stop recombining in the heterogametic sex, thereby enforcing linkage between sex-determining (SD) and sex-antagonistic (SA) genes. With the same rationale, a pre-existing sex asymmetry in recombination is expected to affect the evolution of heterogamety, for example, a low rate of male recombination might favor transitions to XY systems, by generating immediate linkage between SD and SA genes. Furthermore, the accumulation of deleterious mutations on nonrecombining Y chromosomes should favor XY-to-XY transitions (which discard the decayed Y), but disfavor XY-to-ZW transitions (which fix the decayed Y as an autosome). Like many anuran amphibians, Hyla tree frogs have been shown to display drastic heterochiasmy (males only recombine at chromosome tips) and are typically XY, which seems to fit the above expectations. Instead, here we demonstrate that two species, H. sarda and H. savignyi, share a common ZW system since at least 11 Ma. Surprisingly, the typical pattern of restricted male recombination has been maintained since then, despite female heterogamety. Hence, sex chromosomes recombine freely in ZW females, not in ZZ males. This suggests that heterochiasmy does not constrain heterogamety (and vice versa), and that the role of SA genes in the evolution of sex chromosomes might have been overemphasized.

Our assembly-free linkage analysis pipeline (AFLAP) identifies segregating markers as k-mers in the raw reads without using a reference genome assembly for calling variants and provides genotype tables for the construction of unbiased, high-density genetic maps without a genome assembly. AFLAP is validated and contrasted to a conventional workflow using simulated data. AFLAP is applied to whole genome sequencing and genotype-by-sequencing data of F1, F2, and recombinant inbred populations of two different plant species, producing genetic maps that are concordant with genome assemblies. The AFLAP-based genetic map for Bremia lactucae enables the production of a chromosome-scale genome assembly.