Associations between low serum levels of ANRIL and some common gene SNPs in Iranian patients with premature coronary artery disease

Coronary artery disease (CAD) is the major cause of mortality in the world. Premature development of CAD can be attributed to women under 55 and men under 45. Many genetic factors play a part in premature CAD. Among them, ANRIL, a long noncoding RNA is located at the 9p21 risk locus, and its expression seems to be correlated with CAD. In the current study, premature CAD and control blood samples, with and without Type 2 Diabetes (T2D), were genotyped for six SNPs at the 9p21 locus. Additionally, ANRIL serum expression was assessed in both groups using real-time PCR. It was performed using different primers targeting exons 1, 5-6, and 19. The χ2 test for association, along with t-tests and ANOVA, was employed for statistical analysis. In this study, we did not find any significant correlation between premature coronary artery disease and rs10757274, rs2383206, rs2383207, rs496892, rs10757278 and rs10738605. However, a lower ANRIL expression was correlated with each SNP risk genotype. Despite the correlation between lower ANRIL expression and CAD, Type 2 diabetes was associated with higher ANRIL expression. Altogether, the correlation between ANRIL expression and the genotypes of the studied SNPs indicated that genetic variants, even those in intronic regions, affect long noncoding RNA expression levels. In conclusion, we recommend combining genetic variants with expression analysis when developing screening strategies for families with premature CAD. To prevent the devastating outcomes of CAD in young adults, it is crucial to discover noninvasive genetic-based screening tests.

Fine mapping of the Cepaea nemoralis shell colour and mid-banded loci using a high-density linkage map

Molluscs are a highly speciose phylum that exhibits an astonishing array of colours and patterns, yet relatively little progress has been made in identifying the underlying genes that determine phenotypic variation. One prominent example is the land snail Cepaea nemoralis for which classical genetic studies have shown that around nine loci, several physically linked and inherited together as a 'supergene', control the shell colour and banding polymorphism. As a first step towards identifying the genes involved, we used whole-genome resequencing of individuals from a laboratory cross to construct a high-density linkage map, and then trait mapping to identify 95% confidence intervals for the chromosomal region that contains the supergene, specifically the colour locus (C), and the unlinked mid-banded locus (U). The linkage map is made up of 215,593 markers, ordered into 22 linkage groups, with one large group making up ~27% of the genome. The C locus was mapped to a ~1.3 cM region on linkage group 11, and the U locus was mapped to a ~0.7 cM region on linkage group 15. The linkage map will serve as an important resource for further evolutionary and population genomic studies of C. nemoralis and related species, as well as the identification of candidate genes within the supergene and for the mid-banding phenotype.

A family with ulcerative colitis maps to 7p21.1 and comprises a region with regulatory activity for the aryl hydrocarbon receptor gene

We have mapped a locus on chromosome 7p22.3-7p15.3 spanning a 22.4 Mb region for ulcerative colitis (UC) by whole genome linkage analyses of a large Danish family. The family represent three generations with UC segregating as an autosomal dominant trait with variable expressivity. The whole-genome scan resulted in a logarithm of odds score (LOD score) of Z = 3.31, and a whole genome sequencing (WGS) of two affected excluded disease-causing mutations in the protein coding genes. Two rare heterozygote variants, rs182281985:G>A and rs541426369:G>A, both with low allele frequencies (MAF A:0.0001, gnomAD ver3.1.2), were found in clusters of ChiP-seq transcription factors binding sites close to the AHR (aryl hydrocarbon receptor) gene and the UC associated SNP rs1077773:G>A. Testing the two SNPs in a promoter reporter assay for regulatory activity revealed that rs182281985:G>A influenced the AHR promoter. These results suggest a regulatory region that include rs182281985:G>A close to the UC GWAS SNP rs1077773:G>A and further demonstrate evidence that the AHR gene on the 7p-tel region is a candidate susceptible gene for UC.

A missense mutation in Lama3 causes androgen alopecia

Hair loss disorders such as androgenetic alopecia have caused serious disturbances to normal human life. Animal models play an important role in exploring pathogenesis of disease and evaluating new therapies. NIH hairless mice are a spontaneous hairless mouse discovered and bred in our laboratory. In this study, we resequenced the genomes of NIH normal mice and NIH hairless mice and obtained 3,575,560 high-quality, plausible SNP loci and 995,475 InDels. The Euclidean distance algorithm was used to assess the association of SNP loci with the hairless phenotype, at a threshold of 0.62. Two regions of chromosome 18 having the highest association with the phenotype contained 345 genes with a total length of 13.98 Mb. The same algorithm was used to assess the association of InDels with the hairless phenotype at a threshold of 0.54 and revealed a region of 25.45 Mb in length, containing 518 genes. The mutation candidate gene Lama3 (NM_010680.2: c.652C>T; NP_034810.1: p. Arg217Cys) was selected based on the results of functional gene analysis and mutation prediction screening. Lama3 (R217C) mutant mice were further constructed using CRISPR/Cas9 technology, and the relationship between Lama3 point mutations and the hairless phenotype were clarified by phenotypic observation. The results showed that male Lama3 point mutation mice started to lose hair on the 80th day after birth, and the hair loss area gradually expanded over time. H&E staining of skin sections showed that the point mutation mice had increased sebaceous glands in the dermis and missing hair follicle structure (i.e., typical symptoms of androgenetic alopecia). This study is a good extension of the current body of knowledge about the function of Lama3, and the constructed Lama3 (R217C) mutant mice may be a good animal model for studying androgenetic alopecia.

Genetic architecture of individual meiotic crossover rate and distribution in Atlantic Salmon

Meiotic recombination through chromosomal crossovers ensures proper segregation of homologous chromosomes during meiosis, while also breaking down linkage disequilibrium and shuffling alleles at loci located on the same chromosome. Rates of recombination can vary between species, but also between and within individuals, sex and chromosomes within species. Indeed, the Atlantic salmon genome is known to have clear sex differences in recombination with female biased heterochiasmy and markedly different landscapes of crossovers between males and females. In male meiosis, crossovers occur strictly in the telomeric regions, whereas in female meiosis crossovers tend to occur closer to the centromeres. However, little is known about the genetic control of these patterns and how this differs at the individual level. Here, we investigate genetic variation in individual measures of recombination in > 5000 large full-sib families of a Norwegian Atlantic salmon breeding population with high-density SNP genotypes. We show that females had 1.6 × higher crossover counts (CC) than males, with autosomal linkage maps spanning a total of 2174 cM in females and 1483 cM in males. However, because of the extreme telomeric bias of male crossovers, female recombination is much more important for generation of new haplotypes with 8 × higher intra-chromosomal genetic shuffling than males. CC was heritable in females (h2 = 0.11) and males (h2 = 0.10), and shuffling was also heritable in both sex but with a lower heritability in females (h2 = 0.06) than in males (h2 = 0.11). Inter-sex genetic correlations for both traits were close to zero, suggesting that rates and distribution of crossovers are genetically distinct traits in males and females, and that there is a potential for independent genetic change in both sexes in the Atlantic Salmon. Together, these findings give novel insights into the genetic architecture of recombination in salmonids and contribute to a better understanding of how rates and distribution of recombination may evolve in eukaryotes more broadly.

Recombinant inbred lines derived from wide crosses in Pisum

Genomic resources are becoming available for Pisum but to link these to phenotypic diversity requires well marked populations segregating for relevant traits. Here we describe two such resources. Two recombinant inbred populations, derived from wide crosses in Pisum are described. One high resolution mapping population involves cv Caméor, for which the first pea whole genome assembly was obtained, crossed to JI0281, a basally divergent P. sativum sativum landrace from Ethiopia. The other is an inter sub-specific cross between P. s. sativum and the independently domesticated P. s. abyssinicum. The corresponding genetic maps provide information on chromosome level sequence assemblies and identify structural differences between the genomes of these two Pisum subspecies. In order to visualise chromosomal translocations that distinguish the mapping parents, we created a simplified version of Threadmapper to optimise it for interactive 3-dimensional display of multiple linkage groups. The genetic mapping of traits affecting seed coat roughness and colour, plant height, axil ring pigmentation, leaflet number and leaflet indentation enabled the definition of their corresponding genomic regions. The consequence of structural rearrangement for trait analysis is illustrated by leaf serration. These analyses pave the way for identification of the underlying genes and illustrate the utility of these publicly available resources. Segregating inbred populations derived from wide crosses in Pisum, together with the associated marker data, are made publicly available for trait dissection. Genetic analysis of these populations is informative about chromosome scale assemblies, structural diversity in the pea genome and has been useful for the fine mapping of several discrete and quantitative traits.

Deciphering genetic causality between inflammatory bowel disease and periodontitis through bi-directional two-sample Mendelian randomization

Inflammatory bowel disease (IBD) and periodontitis are reported to be closely associated; however, whether there is a causal association between them remains unclear. To explore the existence of this causality, this study applied a bidirectional two-sample Mendelian randomization (MR). The genetic variants were obtained from the summary statistics of genome-wide association studies of IBD, including its subtypes CD and UC, and periodontitis. 175, 148, 113, and six single-nucleotide polymorphisms were selected as instrumental variables for IBD, CD, UC, and periodontitis, respectively. In MR analysis, random-effects inverse-variance weighted was used as the primary method, and weighted median and MR Egger regression were applied as the complementary method. A series of sensitivity analyses were also conducted to ensure the reliability of the results. None of these analyses found a significant effect of genetically proxied IBD and its subtypes on periodontitis, and vice versa. Subsequent sensitivity analyses did not detect any horizontal pleiotropy and heterogeneity. Caution should be exerted when it comes to clinical relevance and further studies are needed to clarify the relationship between IBD and periodontitis.

Molecular markers of type II alveolar epithelial cells in acute lung injury by bioinformatics analysis

In this study, we aimed to identify molecular markers associated with type II alveolar epithelial cell injury in acute lung injury (ALI) models using bioinformatics methods. The objective was to provide new insights for the diagnosis and treatment of ALI/ARDS. We downloaded RNA SEQ datasets (GSE109913, GSE179418, and GSE119123) from the Gene Expression Omnibus (GEO) and used R language package to screen differentially expressed genes (DEGs). DEGs were annotated using Gene Ontology (GO), and their pathways were analyzed using Kyoto Encyclopedia of Genes and Genomes (KEGG). DEGs were imported into the STRING database and analyzed using Cytoscape software to determine the protein network of DEGs and calculate the top 10 nodes for the hub genes. Finally, potential therapeutic drugs for the hub genes were predicted using the DGIdb database. We identified 78 DEGs, including 70 up-regulated genes and 8 down-regulated genes. GO analysis revealed that the DEGs were mainly involved in biological processes such as granulocyte migration, response to bacterial-derived molecules, and cytokine-mediated signaling pathways. Additionally, they had cytokine activity, chemokine activity, and receptor ligand activity, and functioned in related receptor binding, CXCR chemokine receptor binding, G protein-coupled receptor binding, and other molecular functions. KEGG analysis indicated that the DEGs were mainly involved in TNF signaling pathway, IL-17 signaling pathway, NF-κB signal pathway, chemokine signal pathway, cytokine-cytokine receptor interaction signal pathway, and others. We identified eight hub genes, including IRF7, IFIT1, IFIT3, PSMB8, PSMB9, BST2, OASL2, and ZBP1, which were all up-regulated genes. We identified several hub genes of type II alveolar epithelial cells in ALI mouse models using bioinformatics analysis. These results provide new targets for understanding and treating of ALI.

Genetic mapping of quantitative trait loci associated with drought tolerance in chickpea (Cicer arietinum L.)

Elucidation of the genetic basis of drought tolerance is vital for genomics-assisted breeding of drought tolerant crop varieties. Here, we used genotyping-by-sequencing (GBS) to identify single nucleotide polymorphisms (SNPs) in recombinant inbred lines (RILs) derived from a cross between a drought tolerant chickpea variety, Pusa 362 and a drought sensitive variety, SBD 377. The GBS identified a total of 35,502 SNPs and subsequent filtering of these resulted in 3237 high-quality SNPs included in the eight linkage groups. Fifty-one percent of these SNPs were located in the genic regions distributed throughout the genome. The high density linkage map has total map length of 1069 cm with an average marker interval of 0.33 cm. The linkage map was used to identify 9 robust and consistent QTLs for four drought related traits viz. membrane stability index, relative water content, seed weight and yield under drought, with percent variance explained within the range of 6.29%-90.68% and LOD scores of 2.64 to 6.38, which were located on five of the eight linkage groups. A genomic region on LG 7 harbors quantitative trait loci (QTLs) explaining > 90% phenotypic variance for membrane stability index, and > 10% PVE for yield. This study also provides the first report of major QTLs for physiological traits such as membrane stability index and relative water content for drought stress in chickpea. A total of 369 putative candidate genes were identified in the 6.6 Mb genomic region spanning these QTLs. In-silico expression profiling based on the available transcriptome data revealed that 326 of these genes were differentially expressed under drought stress. KEGG analysis resulted in reduction of candidate genes from 369 to 99, revealing enrichment in various signaling pathways. Haplotype analysis confirmed 5 QTLs among the initially identified 9 QTLs. Two QTLs, qRWC1.1 and qYLD7.1, were chosen based on high SNP density. Candidate gene-based analysis revealed distinct haplotypes in qYLD7.1 associated with significant phenotypic differences, potentially linked to pathways for secondary metabolite biosynthesis. These identified candidate genes bolster defenses through flavonoids and phenylalanine-derived compounds, aiding UV protection, pathogen resistance, and plant structure.The study provides novel genomic regions and candidate genes which can be utilized in genomics-assisted breeding of superior drought tolerant chickpea cultivars.

The first linkage map for Australo-Papuan Treefrogs (family: Pelodryadidae) reveals the sex-determination system of the Green-eyed Treefrog (Litoria serrata)

Amphibians represent a useful taxon to study the evolution of sex determination because of their highly variable sex-determination systems. However, the sex-determination system for many amphibian families remains unknown, in part because of a lack of genomic resources. Here, using an F1 family of Green-eyed Treefrogs (Litoria serrata), we produce the first genetic linkage map for any Australo-Papuan Treefrogs (family: Pelodryadidae). The resulting linkage map contains 8662 SNPs across 13 linkage groups. Using an independent set of sexed adults, we identify a small region in linkage group 6 matching an XY sex-determination system. These results suggest Litoria serrata possesses a male heterogametic system, with a candidate sex-determination locus on linkage group 6. Furthermore, this linkage map represents the first genomic resource for Australo-Papuan Treefrogs, an ecologically diverse family of over 220 species.

Functional genomics identify causal variant underlying the protective CTSH locus for Alzheimer's disease

Alzheimer's disease (AD) is the most prevalent age-related neurodegenerative disease, which has a high heritability of up to 79%. Exploring the genetic basis is essential for understanding the pathogenic mechanisms underlying AD development. Recent genome-wide association studies (GWASs) reported an AD-associated signal in the Cathepsin H (CTSH) gene in European populations. However, the exact functional/causal variant(s), and the genetic regulating mechanism of CTSH in AD remain to be determined. In this study, we carried out a comprehensive study to characterize the role of CTSH variants in the pathogenesis of AD. We identified rs2289702 in CTSH as the most significant functional variant that is associated with a protective effect against AD. The genetic association between rs2289702 and AD was validated in independent cohorts of the Han Chinese population. The CTSH mRNA expression level was significantly increased in AD patients and AD animal models, and the protective allele T of rs2289702 was associated with a decreased expression level of CTSH through the disruption of the binding affinity of transcription factors. Human microglia cells with CTSH knockout showed a significantly increased phagocytosis of Aβ peptides. Our study identified CTSH as being involved in AD genetic susceptibility and uncovered the genetic regulating mechanism of CTSH in pathogenesis of AD.

Cloning of the wheat leaf rust resistance gene Lr47 introgressed from Aegilops speltoides

Leaf rust, caused by Puccinia triticina Eriksson (Pt), is one of the most severe foliar diseases of wheat. Breeding for leaf rust resistance is a practical and sustainable method to control this devastating disease. Here, we report the identification of Lr47, a broadly effective leaf rust resistance gene introgressed into wheat from Aegilops speltoides. Lr47 encodes a coiled-coil nucleotide-binding leucine-rich repeat protein that is both necessary and sufficient to confer Pt resistance, as demonstrated by loss-of-function mutations and transgenic complementation. Lr47 introgression lines with no or reduced linkage drag are generated using the Pairing homoeologous1 mutation, and a diagnostic molecular marker for Lr47 is developed. The coiled-coil domain of the Lr47 protein is unable to induce cell death, nor does it have self-protein interaction. The cloning of Lr47 expands the number of leaf rust resistance genes that can be incorporated into multigene transgenic cassettes to control this devastating disease.

Linkage analysis using whole exome sequencing data implicates SLC17A1, SLC17A3, TATDN2 and TMEM131L in type 1 diabetes in Kuwaiti families

Type 1 diabetes (T1D) is characterized by the progressive destruction of pancreatic β-cells, leading to insulin deficiency and lifelong dependency on exogenous insulin. Higher estimates of heritability rates in monozygotic twins, followed by dizygotic twins and sib-pairs, indicate the role of genetics in the pathogenesis of T1D. The incidence and prevalence of T1D are alarmingly high in Kuwait. Consanguineous marriages account for 50-70% of all marriages in Kuwait, leading to an excessive burden of recessive allele enrichment and clustering of familial disorders. Thus, genetic studies from this Arab region are expected to lead to the identification of novel gene loci for T1D. In this study, we performed linkage analyses to identify the recurrent genetic variants segregating in high-risk Kuwaiti families with T1D. We studied 18 unrelated Kuwaiti native T1D families using whole exome sequencing data from 86 individuals, of whom 37 were diagnosed with T1D. The study identified three potential loci with a LOD score of ≥ 3, spanning across four candidate genes, namely SLC17A1 (rs1165196:pT269I), SLC17A3 (rs942379: p.S370S), TATDN2 (rs394558:p.V256I), and TMEM131L (rs6848033:p.R190R). Upon examination of missense variants from these genes in the familial T1D dataset, we observed a significantly increased enrichment of the genotype homozygous for the minor allele at SLC17A3 rs56027330_p.G279R accounting for 16.2% in affected children from 6 unrelated Kuwaiti T1D families compared to 1000 genomes Phase 3 data (0.9%). Data from the NephQTL database revealed that the rs1165196, rs942379, rs394558, and rs56027330 SNPs exhibited genotype-based differential expression in either glomerular or tubular tissues. Data from the GTEx database revealed rs942379 and rs394558 as QTL variants altering the expression of TRIM38 and IRAK2 respectively. Global genome-wide association studies indicated that SLC17A1 rs1165196 and other variants from SLC17A3 are associated with uric acid concentrations and gout. Further evidence from the T1D Knowledge portal supported the role of shortlisted variants in T1D pathogenesis and urate metabolism. Our study suggests the involvement of SLC17A1, SLC17A3, TATDN2, and TMEM131L genes in familial T1D in Kuwait. An enrichment selection of genotype homozygous for the minor allele is observed at SLC17A3 rs56027330_p.G279R variant in affected members of Kuwaiti T1D families. Future studies may focus on replicating the findings in a larger T1D cohort and delineate the mechanistic details of the impact of these novel candidate genes on the pathophysiology of T1D.

Dynamic QTL mapping revealed primarily the genetic structure of photosynthetic traits in castor (Ricinus communis L.)

High photosynthetic efficiency is the basis of high biomass and high harvest index in castor (Ricinus communis L.). Understanding the genetic law of photosynthetic traits will facilitate the breeding for high photosynthetic efficiency. In this study, the dynamic QTL mapping was performed with the populations F2 and BC1 derived from 2 parents with significant difference in net photosynthetic rate (Pn) at 3 stages, in order to reveal the genetic structure of photosynthetic traits. In F2 population, 26 single-locus QTLs were identified, including 3/3/1 (the QTL number at stage I/II/III, the same below), 1/2/0, 1/2/2, 1/3/1, 0/1/1, and 1/1/2 QTLs conferring Pn, water use efficiency (Wue), transpiration rate (Tr), stomatal conductance (Gs), intercellular CO2 concentration (Ci) and chlorophyll content (Cc), with a phenotypic variation explained (PVE) of 8.40%/8.91%/6.17%, 5.36%/31.74%/0, 7.31%/12.80%/15.15%, 1.60%/6.44%/0.02%, 0/1.10%/0.70% and 2.77%/3.96%/6.50% respectively. And 53 epistatic QTLs (31 pairs) were identified, including 2/2/5, 5/6/3, 4/4/2, 6/3/2, 3/2/0 and 4/0/0 ones conferring the above 6 traits, with a PVE of 6.52%/6.47%/19.04%, 16.72%/15.67%/14.12%, 18.57%/15.58%/7.34%, 21.72%/8.52%/7.13%, 13.33%/4.94%/0 and 7.84%/0/0 respectively. The QTL mapping results in BC1 population were consistent with those in F2 population, except fewer QTLs detected. Most QTLs identified were minor-effect ones, only a few were main-effect ones (PVE > 10%), focused on 2 traits, Wue and Tr, such as qWue1.1, qWue1.2, FqTr1.1, FqTr6, BqWue1.1 and BqTr3; The epistatic effects, especially those related to the dominance effects were the main genetic component of photosynthetic traits, and all the epistatic QTLs had no single-locus effects except qPn1.2, FqGs1.2, FqCi1.2 and qCc3.2; The detected QTLs underlying each trait varied at different stages except stable QTLs qGs1.1, detected at 3 stages, qWue2, qTr1.2 and qCc3.2, detected at 2 stages; 6 co-located QTLs were identified, each of which conferring 2-5 different traits, demonstrated the gene pleiotropy between photosynthetic traits; 2 QTL clusters, located within the marker intervals RCM1842-RCM1335 and RCM523-RCM83, contained 15/5 (F2/BC1) and 4/4 (F2/BC1) QTLs conferring multiple traits, including co-located QTLs and main-effect QTLs. The above results provided new insights into the genetic structure of photosynthetic traits and important references for the high photosynthetic efficiency breeding in castor plant.

IHH enhancer variant within neighboring NHEJ1 intron causes microphthalmia anophthalmia and coloboma

Genomic sequences residing within introns of few genes have been shown to act as enhancers affecting expression of neighboring genes. We studied an autosomal recessive phenotypic continuum of microphthalmia, anophthalmia and ocular coloboma, with no apparent coding-region disease-causing mutation. Homozygosity mapping of several affected Jewish Iranian families, combined with whole genome sequence analysis, identified a 0.5 Mb disease-associated chromosome 2q35 locus (maximal LOD score 6.8) harboring an intronic founder variant in NHEJ1, not predicted to affect NHEJ1. The human NHEJ1 intronic variant lies within a known specifically limb-development enhancer of a neighboring gene, Indian hedgehog (Ihh), known to be involved in eye development in mice and chickens. Through mouse and chicken molecular development studies, we demonstrated that this variant is within an Ihh enhancer that drives gene expression in the developing eye and that the identified variant affects this eye-specific enhancer activity. We thus delineate an Ihh enhancer active in mammalian eye development whose variant causes human microphthalmia, anophthalmia and ocular coloboma. The findings highlight disease causation by an intronic variant affecting the expression of a neighboring gene, delineating molecular pathways of eye development.

Predicting recombination suppression outside chromosomal inversions in Drosophila melanogaster using crossover interference theory

Recombination suppression in chromosomal inversion heterozygotes is a well-known but poorly understood phenomenon. Surprisingly, recombination suppression extends far outside of inverted regions where there are no intrinsic barriers to normal chromosome pairing, synapsis, double-strand break formation, or recovery of crossover products. The interference hypothesis of recombination suppression proposes heterozygous inversion breakpoints possess chiasma-like properties such that recombination suppression extends from these breakpoints in a process analogous to crossover interference. This hypothesis is qualitatively consistent with chromosome-wide patterns of recombination suppression extending to both inverted and uninverted regions of the chromosome. The present study generated quantitative predictions for this hypothesis using a probabilistic model of crossover interference with gamma-distributed inter-event distances. These predictions were then tested with experimental genetic data (>40,000 meioses) on crossing-over in intervals that are external and adjacent to four common inversions of Drosophila melanogaster. The crossover interference model accurately predicted the partially suppressed recombination rates in euchromatic intervals outside inverted regions. Furthermore, assuming interference does not extend across centromeres dramatically improved model fit and partially accounted for excess recombination observed in pericentromeric intervals. Finally, inversions with breakpoints closest to the centromere had the greatest excess of recombination in pericentromeric intervals, an observation that is consistent with negative crossover interference previously documented near Drosophila melanogaster centromeres. In conclusion, the experimental data support the interference hypothesis of recombination suppression, validate a mathematical framework for integrating distance-dependent effects of structural heterozygosity on crossover distribution, and highlight the need for improved modeling of crossover interference in pericentromeric regions.

A genetic locus complements resistance to Bordetella pertussis-induced histamine sensitization

Histamine plays pivotal role in normal physiology and dysregulated production of histamine or signaling through histamine receptors (HRH) can promote pathology. Previously, we showed that Bordetella pertussis or pertussis toxin can induce histamine sensitization in laboratory inbred mice and is genetically controlled by Hrh1/HRH1. HRH1 allotypes differ at three amino acid residues with P263-V313-L331 and L263-M313-S331, imparting sensitization and resistance respectively. Unexpectedly, we found several wild-derived inbred strains that carry the resistant HRH1 allotype (L263-M313-S331) but exhibit histamine sensitization. This suggests the existence of a locus modifying pertussis-dependent histamine sensitization. Congenic mapping identified the location of this modifier locus on mouse chromosome 6 within a functional linkage disequilibrium domain encoding multiple loci controlling sensitization to histamine. We utilized interval-specific single-nucleotide polymorphism (SNP) based association testing across laboratory and wild-derived inbred mouse strains and functional prioritization analyses to identify candidate genes for this modifier locus. Atg7, Plxnd1, Tmcc1, Mkrn2, Il17re, Pparg, Lhfpl4, Vgll4, Rho and Syn2 are candidate genes within this modifier locus, which we named Bphse, enhancer of Bordetella pertussis induced histamine sensitization. Taken together, these results identify, using the evolutionarily significant diversity of wild-derived inbred mice, additional genetic mechanisms controlling histamine sensitization.

Predicting recombination frequency from map distance

Map distance is one of the key measures in genetics and indicates the expected number of crossovers between two loci. Map distance is estimated from the observed recombination frequency using mapping functions, the most widely used of those, Haldane and Kosambi, being developed at the time when the number of markers was low and unobserved crossovers had a substantial effect on the recombination fractions. In contemporary high-density marker data, the probability of multiple crossovers between adjacent loci is negligible and different mapping functions yield the same result, that is, the recombination frequency between adjacent loci is equal to the map distance in Morgans. However, high-density linkage maps contain an interpretation problem: the map distance over a long interval is additive and its association with recombination frequency is not defined. Here, we demonstrate with high-density linkage maps from humans and stickleback fishes that the inverses of Haldane's and Kosambi's mapping functions systematically underpredict recombination frequencies from map distance. To remedy this, we formulate a piecewise function that yields more accurate predictions of recombination frequency from map distance. Our results demonstrate that the association between map distance and recombination frequency is context-dependent and without a universal solution.

Genetic mapping of microbial and host traits reveals production of immunomodulatory lipids by Akkermansia muciniphila in the murine gut

The molecular bases of how host genetic variation impacts the gut microbiome remain largely unknown. Here we used a genetically diverse mouse population and applied systems genetics strategies to identify interactions between host and microbe phenotypes including microbial functions, using faecal metagenomics, small intestinal transcripts and caecal lipids that influence microbe-host dynamics. Quantitative trait locus (QTL) mapping identified murine genomic regions associated with variations in bacterial taxa; bacterial functions including motility, sporulation and lipopolysaccharide production and levels of bacterial- and host-derived lipids. We found overlapping QTL for the abundance of Akkermansia muciniphila and caecal levels of ornithine lipids. Follow-up in vitro and in vivo studies revealed that A. muciniphila is a major source of these lipids in the gut, provided evidence that ornithine lipids have immunomodulatory effects and identified intestinal transcripts co-regulated with these traits including Atf3, which encodes for a transcription factor that plays vital roles in modulating metabolism and immunity. Collectively, these results suggest that ornithine lipids are potentially important for A. muciniphila-host interactions and support the role of host genetics as a determinant of responses to gut microbes.

QTL mapping for seed density per silique in Brassica napus

Seed density per silique (SDPS) and valid silique length (VSL) are two important yield-influencing traits in rapeseed. SDPS has a direct or indirect effect on rapeseed yield through its effect on seed per silique. In this study, a quantitative trait locus (QTL) for SDPS was detected on chromosome A09 using the QTL-seq approach and confirmed via linkage analysis in the mapping population obtained from 4263 × 3001 cross. Furthermore, one major QTL for SDPS (qSD.A9-1) was mapped to a 401.8 kb genomic interval between SSR markers Nys9A190 and Nys9A531. In the same genomic region, a QTL (qSL.A9) linked to VSL was also detected. The phenotypic variation of qSD.A9-1 and qSL.A9 was 53.1% and 47.6%, respectively. Results of the additive and dominant effects demonstrated that the expression of genes controlling SDPS and VSL were derived from a different parent in this population. Subsequently, we identified 56 genes that included 45 specific genes with exonic (splicing) variants. Further analysis identified specific genes containing mutations that may be related to seed density as well as silique length. These genes could be used for further studies to understand the details of these traits of rapeseed.

Compatibility between snails and schistosomes: insights from new genetic resources, comparative genomics, and genetic mapping

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.

Mycobacterium abscessus pathogenesis identified by phenogenomic analyses

The medical and scientific response to emerging and established pathogens is often severely hampered by ignorance of the genetic determinants of virulence, drug resistance and clinical outcomes that could be used to identify therapeutic drug targets and forecast patient trajectories. Taking the newly emergent multidrug-resistant bacteria Mycobacterium abscessus as an example, we show that combining high-dimensional phenotyping with whole-genome sequencing in a phenogenomic analysis can rapidly reveal actionable systems-level insights into bacterial pathobiology. Through phenotyping of 331 clinical isolates, we discovered three distinct clusters of isolates, each with different virulence traits and associated with a different clinical outcome. We combined genome-wide association studies with proteome-wide computational structural modelling to define likely causal variants, and employed direct coupling analysis to identify co-evolving, and therefore potentially epistatic, gene networks. We then used in vivo CRISPR-based silencing to validate our findings and discover clinically relevant M. abscessus virulence factors including a secretion system, thus illustrating how phenogenomics can reveal critical pathways within emerging pathogenic bacteria.

A dominant negative ADIPOQ mutation in a diabetic family with renal disease, hypoadiponectinemia, and hyperceramidemia

Adiponectin, encoded by ADIPOQ, is an insulin-sensitizing, anti-inflammatory, and renoprotective adipokine that activates receptors with intrinsic ceramidase activity. We identified a family harboring a 10-nucleotide deletion mutation in ADIPOQ that cosegregates with diabetes and end-stage renal disease. This mutation introduces a frameshift in exon 3, resulting in a premature termination codon that disrupts translation of adiponectin's globular domain. Subjects with the mutation had dramatically reduced circulating adiponectin and increased long-chain ceramides levels. Functional studies suggest that the mutated protein acts as a dominant negative through its interaction with non-mutated adiponectin, decreasing circulating adiponectin levels, and correlating with metabolic disease.

Genome-wide linkage analysis combined with genome sequencing in large families with intracranial aneurysms

Rupture of an intracranial aneurysm (IA) leads to aneurysmal subarachnoid haemorrhage (ASAH), a severe type of stroke. Some rare variants that cause IA in families have been identified, but still, the majority of genetic causes, as well as the biological mechanisms of IA development and rupture, remain unknown. We aimed to identify rare, damaging variants for IA in three large Dutch families with multiple affected members with IA (N = 9, 11, and 6). By combining linkage analysis and genome sequencing (GS), we identified six rare and damaging variants for which all cases within one of the families were heterozygous. These variants were p.Tyr87Cys in SYCP1, p.Phe1077Leu in FMNL2, p.Thr754Lys in TBC1D2, p.Arg321His in ZNF782, p.Arg979Trp in CCDC180, and p.Val125Met in NCBP1. None of the variants showed association with IA status in a large cohort of 937 patients from the general IA patient population and 1046 controls. Gene expression in IA and cerebral artery tissue further prioritized FMNL2 and TBC1D2 as potential important players in IA pathophysiology. Further studies are needed to characterize the functional consequences of the identified variants and their role in the biological mechanisms of IA.

Chromosome evolution and the genetic basis of agronomically important traits in greater yam

The nutrient-rich tubers of the greater yam, Dioscorea alata L., provide food and income security for millions of people around the world. Despite its global importance, however, greater yam remains an orphan crop. Here, we address this resource gap by presenting a highly contiguous chromosome-scale genome assembly of D. alata combined with a dense genetic map derived from African breeding populations. The genome sequence reveals an ancient allotetraploidization in the Dioscorea lineage, followed by extensive genome-wide reorganization. Using the genomic tools, we find quantitative trait loci for resistance to anthracnose, a damaging fungal pathogen of yam, and several tuber quality traits. Genomic analysis of breeding lines reveals both extensive inbreeding as well as regions of extensive heterozygosity that may represent interspecific introgression during domestication. These tools and insights will enable yam breeders to unlock the potential of this staple crop and take full advantage of its adaptability to varied environments.

The interplay of additivity, dominance, and epistasis on fitness in a diploid yeast cross

In diploid species, genetic loci can show additive, dominance, and epistatic effects. To characterize the contributions of these different types of genetic effects to heritable traits, we use a double barcoding system to generate and phenotype a panel of ~200,000 diploid yeast strains that can be partitioned into hundreds of interrelated families. This experiment enables the detection of thousands of epistatic loci, many whose effects vary across families. Here, we show traits are largely specified by a small number of hub loci with major additive and dominance effects, and pervasive epistasis. Genetic background commonly influences both the additive and dominance effects of loci, with multiple modifiers typically involved. The most prominent dominance modifier in our data is the mating locus, which has no effect on its own. Our findings show that the interplay between additivity, dominance, and epistasis underlies a complex genotype-to-phenotype map in diploids.

Identification of quantitative trait loci associated with upper temperature tolerance in turbot, Scophthalmus maximus

Temperature tolerance is an important trait from both an economic and evolutionary perspective in fish. Because of difficulties with measurements, genome-wide selection using quantitative trait loci (QTLs) affecting Upper temperature tolerance may be an alternative for genetic improvement. Turbot Scophthalmus maximus (L.) is a cold-water marine fish with high economic value in Europe and Asia. The genetic bases of upper temperature tolerance (UTTs) traits have been rarely studied. In this study, we constructed a genetic linkage map of turbot using simple sequence repeats (SSRs) and single nucleotide polymorphism (SNP) markers. A total of 190 SSR and 8,123 SNP were assigned to 22 linkage groups (LGs) of a consensus map, which spanned 3,648.29 cM of the turbot genome, with an average interval of 0.44 cM. Moreover, we re-anchored genome sequences, allowing 93.8% physical sequences to be clustered into 22 turbot pseudo-chromosomes. A high synteny was observed between two assemblies from the literature. QTL mapping and validation analysis identified thirteen QLTs which are major effect QTLs, of these, 206 linked SNP loci, and two linked SSR loci were considered to have significant QTL effects. Association analysis for UTTs with 129 QTL markers was performed for different families, results showed that eight SNP loci were significantly correlated with UTT, which markers could be helpful in selecting thermal tolerant breeds of turbot. 1,363 gene sequences were genomically annotated, and 26 QTL markers were annotated. We believe these genes could be valuable candidates affecting high temperatures, providing valuable genomic resources for the study of genetic mechanisms regulating thermal stress. Similarly, they may be used in marker-assisted selection (MAS) programs to improve turbot performance.

Major QTL confer race-nonspecific resistance in the co-evolved Cronartium quercuum f. sp. fusiforme-Pinus taeda pathosystem

Fusiform rust disease, caused by the endemic fungus Cronartium quercuum f. sp. fusiforme, is the most damaging disease affecting economically important pine species in the southeast United States. Unlike the major epidemics of agricultural crops, the co-evolved pine-rust pathosystem is characterized by steady-state dynamics and high levels of genetic diversity within environments. This poses a unique challenge and opportunity for the deployment of large-effect resistance genes. We used trait dissection to study the genetic architecture of disease resistance in two P. taeda parents that showed high resistance across multiple environments. Two mapping populations (full-sib families), each with ~1000 progeny, were challenged with a complex inoculum consisting of 150 pathogen isolates. High-density linkage mapping revealed three major-effect QTL distributed on two linkage groups. All three QTL were validated using a population of 2057 cloned pine genotypes in a 6-year-old multi-environmental field trial. As a complement to the QTL mapping approach, bulked segregant RNAseq analysis revealed a small number of candidate nucleotide binding leucine-rich repeat genes harboring SNP associated with disease resistance. The results of this study show that in P. taeda, a small number of major QTL can provide effective resistance against genetically diverse mixtures of an endemic pathogen. These QTL vary in their impact on disease liability and exhibit additivity in combination.

A 127 kb truncating deletion of PGRMC1 is a novel cause of X-linked isolated paediatric cataract

Inherited paediatric cataract is a rare Mendelian disease that results in visual impairment or blindness due to a clouding of the eye's crystalline lens. Here we report an Australian family with isolated paediatric cataract, which we had previously mapped to Xq24. Linkage at Xq24-25 (LOD = 2.53) was confirmed, and the region refined with a denser marker map. In addition, two autosomal regions with suggestive evidence of linkage were observed. A segregating 127 kb deletion (chrX:g.118373226_118500408del) in the Xq24-25 linkage region was identified from whole-genome sequencing data. This deletion completely removed a commonly deleted long non-coding RNA gene LOC101928336 and truncated the protein coding progesterone receptor membrane component 1 (PGRMC1) gene following exon 1. A literature search revealed a report of two unrelated males with non-syndromic intellectual disability, as well as congenital cataract, who had contiguous gene deletions that accounted for their intellectual disability but also disrupted the PGRMC1 gene. A morpholino-induced pgrmc1 knockdown in a zebrafish model produced significant cataract formation, supporting a role for PGRMC1 in lens development and cataract formation. We hypothesise that the loss of PGRMC1 causes cataract through disrupted PGRMC1-CYP51A1 protein-protein interactions and altered cholesterol biosynthesis. The cause of paediatric cataract in this family is the truncating deletion of PGRMC1, which we report as a novel cataract gene.

Quantitative trait locus analysis of parasitoid counteradaptation to symbiont-conferred resistance

Insect hosts and parasitoids are engaged in an intense struggle of antagonistic coevolution. Infection with heritable bacterial endosymbionts can substantially increase the resistance of aphids to parasitoid wasps, which exerts selection on parasitoids to overcome this symbiont-conferred protection (counteradaptation). Experimental evolution in the laboratory has produced counteradapted populations of the parasitoid wasp Lysiphlebus fabarum. These populations can parasitize black bean aphids (Aphis fabae) protected by the bacterial endosymbiont Hamiltonella defensa, which confers high resistance against L. fabarum. We used two experimentally evolved parasitoid populations to study the genetic architecture of the counteradaptation to symbiont-conferred resistance by QTL analysis. With simple crossing experiments, we showed that the counteradaptation is a recessive trait depending on the maternal genotype. Based on these results, we designed a customized crossing scheme to genotype a mapping population phenotyped for the ability to parasitize Hamiltonella-protected aphids. Using 1835 SNP markers obtained by ddRAD sequencing, we constructed a high-density linkage map consisting of six linkage groups (LGs) with an overall length of 828.3 cM and an average marker spacing of 0.45 cM. We identified a single QTL associated with the counteradaptation to Hamiltonella in L. fabarum on linkage group 2. Out of 120 genes located in this QTL, several genes encoding putative venoms may represent candidates for counteradaptation, as parasitoid wasps inject venoms into their hosts during oviposition.

A Decade in Review after Idiopathic Scoliosis Was First Called a Complex Trait-A Tribute to the Late Dr. Yves Cotrel for His Support in Studies of Etiology of Scoliosis

Adolescent Idiopathic Scoliosis (AIS) is a prevalent and important spine disorder in the pediatric age group. An increased family tendency was observed for a long time, but the underlying genetic mechanism was uncertain. In 1999, Dr. Yves Cotrel founded the Cotrel Foundation in the Institut de France, which supported collaboration of international researchers to work together to better understand the etiology of AIS. This new concept of AIS as a complex trait evolved in this setting among researchers who joined the annual Cotrel meetings. It is now over a decade since the first proposal of the complex trait genetic model for AIS. Here, we review in detail the vast information about the genetic and environmental factors in AIS pathogenesis gathered to date. More importantly, new insights into AIS etiology were brought to us through new research data under the perspective of a complex trait. Hopefully, future research directions may lead to better management of AIS, which has a tremendous impact on affected adolescents in terms of both physical growth and psychological development.

'Atlantic'-derived F(1) population

There are many challenges involved with the genetic analyses of autopolyploid species, such as the tetraploid potato, Solanum tuberosum (2n = 4x = 48). The development of new analytical methods has made it valuable to re-analyze an F₁ population (n = 156) derived from a cross involving 'Atlantic', a widely grown chipping variety in the USA. A fully integrated genetic map with 4285 single nucleotide polymorphisms, spanning 1630 cM, was constructed with MAPpoly software. We observed that bivalent configurations were the most abundant ones (51.0~72.4% depending on parent and linkage group), though multivalent configurations were also observed (2.2~39.2%). Seven traits were evaluated over four years (2006-8 and 2014) and quantitative trait loci (QTL) mapping was carried out using QTLpoly software. Based on a multiple-QTL model approach, we detected 21 QTL for 15 out of 27 trait-year combination phenotypes. A hotspot on linkage group 5 was identified with co-located QTL for maturity, plant yield, specific gravity, and internal heat necrosis resistance evaluated over different years. Additional QTL for specific gravity and dry matter were detected with maturity-corrected phenotypes. Among the genes around QTL peaks, we found those on chromosome 5 that have been previously implicated in maturity (StCDF1) and tuber formation (POTH1). These analyses have the potential to provide insights into the biology and breeding of tetraploid potato and other autopolyploid species.

Association between MIF gene promoter rs755622 and susceptibility to coronary artery disease and inflammatory cytokines in the Chinese Han population

Macrophage migration inhibitory factor (MIF) is an essential mediator of atherosclerotic plaque progression and instability leading to intracoronary thrombosis, therefore contributing to coronary artery disease (CAD). In this study, we investigated the relationship between MIF gene polymorphism and CAD in Chinese Han population. Three single nucleotide polymorphisms (SNP, rs755622, rs1007888 and rs2096525) of MIF gene were genotyped by TaqMan genotyping assay in 1120 control participants and 1176 CAD patients. Coronary angiography was performed in all CAD patients and Gensini score was used to assess the severity of coronary artery lesions. The plasma levels of MIF and other inflammatory mediators were measured by ELISA. The CAD patients had a higher frequency of CC genotype and C allele of rs755622 compared with that in control subjects (CC genotype: 6.5% vs. 3.9%, P = 0.008, C allele: 24.0% vs. 20.6%, P = 0.005). The rs755622 CC genotype was associated with an increased risk of CAD (OR: 1.804, 95%CI: 1.221-2.664, P = 0.003). CAD patients with a variation of rs755622 CC genotype had significantly higher Gensini score compared with patients with GG or CG genotype (all P < 0.05). In addition, the circulating MIF level was highest in CAD patients carrying rs755622 CC genotype (40.7 ± 4.2 ng/mL) and then followed by GC (37.9 ± 3.4 ng/mL) or GG genotype (36.9 ± 3.7 ng/mL, all P < 0.01). Our study showed an essential relationship between the MIF gene rs755622 variation and CAD in Chinese Han population. Individuals who carrying MIF gene rs755622 CC genotype were more susceptible to CAD and had more severe coronary artery lesion. This variation also had a potential influence in circulating MIF levels.

Genome survey and high-resolution genetic map provide valuable genetic resources for Fenneropenaeus chinensis

Fenneropenaeus chinensis is one of the most important aquaculture species in China. Research on its genomic and genetic structure not only helps us comprehend the genetic basis of complex economic traits, but also offers theoretical guidance in selective breeding. In the present study, a genome survey sequencing was performed to generate a rough reference genome utilized for groping preliminary genome characteristics and facilitate linkage and quantitative trait locus (QTL) mapping. Linkage mapping was conducted using a reduced-representation sequencing method 2b-RAD. In total, 36,762 SNPs were genotyped from 273 progenies in a mapping family, and a high-resolution linkage map was constructed. The consensus map contained 12,884 markers and spanned 5257.81 cM with an average marker interval of 0.41 cM, which was the first high-resolution genetic map in F. chinensis to our knowledge. QTL mapping and association analysis were carried out in 29 characters including body size, sex and disease resistance. 87 significant QTLs were detected in several traits and they were also evaluated by association analysis. Results of this study provide us valuable suggestions in genetic improvement and breeding of new varieties and also lay a basic foundation for further application of cloning of economic genes in selective breeding program and marker-assisted selection.

Genetic dissection of down syndrome-associated alterations in APP/amyloid-β biology using mouse models

Individuals who have Down syndrome (caused by trisomy of chromosome 21), have a greatly elevated risk of early-onset Alzheimer's disease, in which amyloid-β accumulates in the brain. Amyloid-β is a product of the chromosome 21 gene APP (amyloid precursor protein) and the extra copy or 'dose' of APP is thought to be the cause of this early-onset Alzheimer's disease. However, other chromosome 21 genes likely modulate disease when in three-copies in people with Down syndrome. Here we show that an extra copy of chromosome 21 genes, other than APP, influences APP/Aβ biology. We crossed Down syndrome mouse models with partial trisomies, to an APP transgenic model and found that extra copies of subgroups of chromosome 21 gene(s) modulate amyloid-β aggregation and APP transgene-associated mortality, independently of changing amyloid precursor protein abundance. Thus, genes on chromosome 21, other than APP, likely modulate Alzheimer's disease in people who have Down syndrome.

Detection of breeding signatures in wheat using a linkage disequilibrium-corrected mapping approach

Marker assisted breeding, facilitated by reference genome assemblies, can help to produce cultivars adapted to changing environmental conditions. However, anomalous linkage disequilibrium (LD), where single markers show high LD with markers on other chromosomes but low LD with adjacent markers, is a serious impediment for genetic studies. We used a LD-correction approach to overcome these drawbacks, correcting the physical position of markers derived from 15 and 135 K arrays in a diversity panel of bread wheat representing 50 years of breeding history. We detected putative mismapping of 11.7% markers and improved the physical alignment of 5.4% markers. Population analysis indicated reduced genetic diversity over time as a result of breeding efforts. By analysis of outlier loci and allele frequency change over time we traced back the 2NS/2AS translocation of Aegilops ventricosa to one cultivar, "Cardos" (registered in 1998) which was the first among the panel to contain this translocation. A "selective sweep" for this important translocation region on chromosome 2AS was found, putatively linked to plant response to biotic stress factors. Our approach helps in overcoming the drawbacks of incorrectly anchored markers on the wheat reference assembly and facilitates detection of selective sweeps for important agronomic traits.

Identification and mapping of quantitative trait loci for resistance to Liriomyza trifolii in romaine lettuce cultivar 'Valmaine'

Liriomyza trifolii (Diptera: Agromyzidae) is a leafminer that causes ruinous damage to many leafy vegetables including lettuce (Lactuca sativa L.) by stippling and tunneling the leaves. In this study, a population of 125 F3 families was developed from the intraspecific cross of 'Valmaine' (resistant) and 'Okeechobee' (susceptible) romaine cultivars for inheritance analysis and molecular mapping of the resistance loci controlling stippling damage. The experiments were conducted in an insectarium (controlled environment). Stippling damage proved to be heritable because the broad-sense heritability (H2) was 0.58. A segregation analysis suggested that a single dominant allele, Sd1 locus, controls resistance against L. trifolii. Furthermore, a quantitative trait loci (QTL) analysis identified one novel QTL, named Stippling on LG5 (qSTP5), flanked by two SNPs that were mapped to a 5.2 cM (8.5 Mb region) interval, explaining over 13% of the total phenotypic variance. Desirable allele for resistance to L. trifolii was derived from resistant cultivar Valmaine. Identification of SNPs closely linked to the QTL responsible for L. trifolii resistance should facilitate plant breeders to develop resistant romaine lettuce cultivars.

A quantitative trait rare variant nonparametric linkage method with application to age-at-onset of Alzheimer's disease

To analyze pedigrees with quantitative trait (QT) and sequence data, we developed a rare variant (RV) quantitative nonparametric linkage (QNPL) method, which evaluates sharing of minor alleles. RV-QNPL has greater power than the traditional QNPL that tests for excess sharing of minor and major alleles. RV-QNPL is robust to population substructure and admixture, locus heterogeneity, and inclusion of nonpathogenic variants and can be readily applied outside of coding regions. When QNPL was used to analyze common variants, it often led to loci mapping to large intervals, e.g., >40 Mb. In contrast, when RVs are analyzed, regions are well defined, e.g., a gene. Using simulation studies, we demonstrate that RV-QNPL is substantially more powerful than applying traditional QNPL methods to analyze RVs. RV-QNPL was also applied to analyze age-at-onset (AAO) data for 107 late-onset Alzheimer's disease (LOAD) pedigrees of Caribbean Hispanic and European ancestry with whole-genome sequence data. When AAO of AD was analyzed regardless of APOE ε4 status, suggestive linkage (LOD = 2.4) was observed with RVs in KNDC1 and nominally significant linkage (p < 0.05) was observed with RVs in LOAD genes ABCA7 and IQCK. When AAO of AD was analyzed for APOE ε4 positive family members, nominally significant linkage was observed with RVs in APOE, while when AAO of AD was analyzed for APOE ε4 negative family members, nominal significance was observed for IQCK and ADAMTS1. RV-QNPL provides a powerful resource to analyze QTs in families to elucidate their genetic etiology.

The bHLH transcription factor VfTT8 underlies zt2, the locus determining zero tannin content in faba bean (Vicia faba L.)

Faba bean (Vicia faba L.) is an important protein-rich fodder crop, which is widely cultivated in temperate areas. However, antinutritional compounds such as condensed tannins, limit the use of this protein source in monogastric feed formulations. Previous studies demonstrated that two recessive and complementary genes, zt1 and zt2, control absence of tannin and white flower colour in faba bean. An ortholog of the Medicago WD40 transcription factor TTG1 was reported to encode the zt1 phenotype, but the responsible gene for zt2 is still unknown. Here we used a candidate gene approach combined with linkage mapping, comparative genomics and gene expression to fine map the zt2 genomic region and to identify the regulatory gene controlling both traits. Seventy-two genes, including 23 MYB and bHLH regulatory genes predicted to be associated with anthocyanin expression together with WRKY proteins, were screened and genotyped in three mapping populations. The linkage groups constructed identified the regulatory gene, TRANSPARENT TESTA8 (TT8), encoding a basic helix-loop-helix (bHLH) transcription factor, as the candidate for zt2. This finding was supported by qPCR analysis and further validated in different genetic backgrounds. Accordingly, VfTT8 was downregulated in white flowered types while showing high levels of expression in wild genotypes. Our results provide new insights on the regulatory mechanisms of tannin biosynthesis in faba bean and will facilitate the development of an ultimate zt2 diagnostic marker for the fast generation of new value-added cultivars free of tannins and with improved nutritional value.

Effect of heterogeneity in recombination rate on variation in realised relationship

Individuals of a specified pedigree relationship vary in the proportion of the genome they share identical by descent, i.e. in their realised or actual relationship. Predictions of the variance in realised relationship have previously been based solely on the proportion of the map length shared, which requires the implicit assumption that both recombination rate and genetic information are uniformly distributed along the genome. This ignores the possible existence of recombination hotspots, and fails to distinguish between coding and non-coding sequences. In this paper, we therefore quantify the effects of heterogeneity in recombination rate at broad and fine-scale levels on the variation in realised relationship. Variance is usually greater on a chromosome with a non-uniform recombination rate than on a chromosome with the same map length and uniform recombination rate, especially if recombination rates are higher towards chromosome ends. Reductions in variance can also be obtained, however, and the overall pattern of change is quite complex. In general, local (fine-scale) variation in recombination rate, e.g. hotspots, has a small influence on the variance in realised relationship. Differences in rates across longer regions and between chromosome ends can increase or decrease the variance in a realised relationship, depending on the genomic architecture.

Effect of marker segregation distortion on high density linkage map construction and QTL mapping in Soybean (Glycine max L.)

Marker segregation distortion is a natural phenomenon. Severely distorted markers are usually excluded in the construction of linkage maps. We investigated the effect of marker segregation distortion on linkage map construction and quantitative trait locus (QTL) mapping. A total of 519 recombinant inbred lines of soybean from orthogonal and reciprocal crosses between LSZZH and NN493-1 were genotyped by specific length amplified fragment markers and seed linoleic acid content was measured in three environments. As a result, twenty linkage groups were constructed with 11,846 markers, including 1513 (12.77%) significantly distorted markers, on 20 chromosomes, and the map length was 2475.86 cM with an average marker-interval of 0.21 cM. The inclusion of distorted markers in the analysis was shown to not only improve the grouping of the markers from the same chromosomes, and the consistency of linkage maps with genome, but also increase genome coverage by markers. Combining genotypic data from both orthogonal and reciprocal crosses decreased the proportion of distorted markers and then improved the quality of linkage maps. Validation of the linkage maps was confirmed by the high collinearity between positions of markers in the soybean reference genome and in linkage maps and by the high consistency of 24 QTL regions in this study compared with the previously reported QTLs and lipid metabolism related genes. Additionally, linkage maps that include distorted markers could add more information to the outputs from QTL mapping. These results provide important information for linkage mapping, gene cloning and marker-assisted selection in soybean.

Scrapie susceptibility-associated indel polymorphism of shadow of prion protein gene (SPRN) in Korean native black goats

Prion diseases in sheep and goats are called scrapie and belong to a group of transmissible spongiform encephalopathies (TSEs) caused by the abnormal misfolding of the prion protein encoded by the prion protein gene (PRNP). The shadow of the prion protein gene (SPRN) is the only prion gene family member that shows a protein expression profile similar to that of the PRNP gene in the central nervous system. In addition, genetic susceptibility of the SPRN gene has been reported in variant Creutzfeldt-Jakob disease (CJD), bovine spongiform encephalopathy (BSE) and scrapie. However, genetic studies of the SPRN gene have not been carried out in Korean native black goats. Here, we investigated the genotype and allele frequencies of SPRN polymorphisms in 213 Korean native black goats and compared these polymorphisms with those previously reported for scrapie-affected animals. We found a total of 6 polymorphisms including 1 nonsynonymous single nucleotide polymorphism (SNP) and 1 synonymous SNP in the open reading frame (ORF) region and 3 SNPs and 1 indel polymorphism (c.495_496insCTCCC) in the 3' untranslated region (UTR) by direct DNA sequencing. A significant difference in the allele frequency of the c.495_496insCTCCC indel polymorphism was found between the Italian scrapie-affected goats and the Korean native black goats (P < 0.001). Furthermore, there was a significant difference in the allele frequencies of the c.495_496insCTCCC indel polymorphism between Italian healthy goats and Korean native black goats (P < 0.001). To evaluate the biological impact of the novel nonsynonymous SNP c.416G > A (Arg139Gln), we carried out PROVEAN analysis. PROVEAN predicted the SNP as 'Neutral' with a score of -0.297. To the best of our knowledge, this is the first genetic study of the SPRN gene in Korean native black goats.

Identification of a gene encoding polygalacturonase expressed specifically in short styles in distylous common buckwheat (Fagopyrum esculentum)

Common buckwheat (Fagopyrum esculentum) is a heteromorphic self-incompatible (SI) species with two types of floral architecture: thrum (short style) and pin (long style). The floral morphology and intra-morph incompatibility are controlled by a single genetic locus, S. However, the molecular mechanisms underlying the heteromorphic self-incompatibility of common buckwheat remain unclear. To identify these mechanisms, we performed proteomic, quantitative reverse-transcription PCR, and linkage analyses. Comparison of protein profiles between the long and short styles revealed a protein unique to the short style. Amino-acid sequencing revealed that it was a truncated form of polygalacturonase (PG); we designated the gene encoding this protein FePG1. Phylogenetic analysis classified FePG1 into the same clade as PGs that function in pollen development and floral morphology. FePG1 expression was significantly higher in short styles than in long styles. It was expressed in flowers of a short-homostyle line but not in flowers of a long-homostyle line. Linkage analysis indicated that FePG1 was not linked to the S locus; it could be a factor downstream of this locus. Our finding of a gene putatively working under the regulation of the S locus provides useful information for elucidation of the mechanism of heteromorphic self-incompatibility.

Fast genetic mapping of complex traits in C. elegans using millions of individuals in bulk

Genetic studies of complex traits in animals have been hindered by the need to generate, maintain, and phenotype large panels of recombinant lines. We developed a new method, C. elegans eXtreme Quantitative Trait Locus (ceX-QTL) mapping, that overcomes this obstacle via bulk selection on millions of unique recombinant individuals. We use ceX-QTL to map a drug resistance locus with high resolution. We also map differences in gene expression in live worms and discovered that mutations in the co-chaperone sti-1 upregulate the transcription of HSP-90. Lastly, we use ceX-QTL to map loci that influence fitness genome-wide confirming previously reported causal variants and uncovering new fitness loci. ceX-QTL is fast, powerful and cost-effective, and will accelerate the study of complex traits in animals.

Transmission ratio distortion is frequent in Arabidopsis thaliana controlled crosses

The equal probability of transmission of alleles from either parent during sexual reproduction is a central tenet of genetics and evolutionary biology. Yet, there are many cases where this rule is violated. The preferential transmission of alleles or genotypes is termed transmission ratio distortion (TRD). Examples of TRD have been identified in many species, implying that they are universal, but the resolution of species-wide studies of TRD are limited. We have performed a species-wide screen for TRD in over 500 segregating F2 populations of Arabidopsis thaliana using pooled reduced-representation genome sequencing. TRD was evident in up to a quarter of surveyed populations. Most populations exhibited distortion at only one genomic region, with some regions being repeatedly affected in multiple populations. Our results begin to elucidate the species-level architecture of biased transmission of genetic material in A. thaliana, and serve as a springboard for future studies into the biological basis of TRD in this species.

From genome-wide associations to candidate causal variants by statistical fine-mapping

Advancing from statistical associations of complex traits with genetic markers to understanding the functional genetic variants that influence traits is often a complex process. Fine-mapping can select and prioritize genetic variants for further study, yet the multitude of analytical strategies and study designs makes it challenging to choose an optimal approach. We review the strengths and weaknesses of different fine-mapping approaches, emphasizing the main factors that affect performance. Topics include interpreting results from genome-wide association studies (GWAS), the role of linkage disequilibrium, statistical fine-mapping approaches, trans-ethnic studies, genomic annotation and data integration, and other analysis and design issues.