Association mapping of resistance to Puccinia hordei in Australian barley breeding germplasm

Key message “To find stable resistance using association mapping tools, QTL with major and minor effects on leaf rust reactions were identified in barley breeding lines by assessing seedlings and adult plants.” Abstract Three hundred and sixty (360) elite barley (Hordeum vulgare L.) breeding lines from the Northern Region Barley Breeding Program in Australia were genotyped with 3,244 polymorphic diversity arrays technology markers and the results used to map quantitative trait loci (QTL) conferring a reaction to leaf rust (Puccinia hordei Otth). The F3:5 (Stage 2) lines were derived or sourced from different geographic origins or hubs of international barley breeding ventures representing two breeding cycles (2009 and 2011 trials) and were evaluated across eight environments for infection type at both seedling and adult plant stages. Association mapping was performed using mean scores for disease reaction, accounting for family effects using the eigenvalues from a matrix of genotype correlations. In this study, 15 QTL were detected; 5 QTL co-located with catalogued leaf rust resistance genes (Rph1, Rph3/19, Rph8/14/15, Rph20, Rph21), 6 QTL aligned with previously reported genomic regions and 4 QTL (3 on chromosome 1H and 1 on 7H) were novel. The adult plant resistance gene Rph20 was identified across the majority of environments and pathotypes. The QTL detected in this study offer opportunities for breeding for more durable resistance to leaf rust through pyramiding multiple genomic regions via marker-assisted selection.

Spot form of net blotch resistance in barley is under complex genetic control

Key message: Evaluation of resistance toPyrenophora teresf.maculatain barley breeding populations via association mapping revealed a complex genetic architecture comprising a mixture of major and minor effect genes. Abstract: In the search for stable resistance to spot form of net blotch (Pyrenophora teres f. maculata, SFNB), association mapping was conducted on four independent barley (Hordeum vulgare L.) breeding populations comprising a total of 898 unique elite breeding lines from the Northern Region Barley Breeding Program in Australia for discovery of quantitative trait loci (QTL) influencing resistance at seedling and adult plant growth stages. A total of 29 significant QTL were validated across multiple breeding populations, with 22 conferring resistance at both seedling and adult plant growth stages. The remaining 7 QTL conferred resistance at either seedling (2 QTL) or adult plant (5 QTL) growth stages only. These 29 QTL represented 24 unique genomic regions, of which five were found to co-locate with previously identified QTL for SFNB. The results indicated that SFNB resistance is controlled by a large number of QTL varying in effect size with large effects QTL on chromosome 7H. A large proportion of the QTL acted in the same direction for both seedling and adult responses, suggesting that phenotypic selection for SFNB resistance performed at either growth stage could achieve adequate levels of resistance. However, the accumulation of specific resistance alleles on several chromosomes must be considered in molecular breeding selection strategies.

Identifying the Function of Sorghum's Drought Tolerance Stay-Green QTL

The goal of this research is to understand the function of allelic variation of genes underpinning the stay-green drought adaptation trait in sorghum in order to enhance yield in water-limited environments. Stay-green, a delayed leaf senescence phenotype in sorghum, is primarily an emergent consequence of the improved balance between the supply and demand of water. Positional and functional fine-mapping of candidate genes associated with stay-green in sorghum is the focus of an international research partnership between Australian (UQ/DAFFQ) and US (Texas A&M University) scientists. Stay-green was initially mapped to four chromosomal regions (Stg1, Stg2, Stg3, and Stg4) by a number of research groups in the US and Australia. Physiological dissection of near-isolines containing single introgressions of Stg QTL (Stg1-4) indicate that these QTL reduce water demand before flowering by constricting the size of the canopy, thereby increasing water availability during grain filling and, ultimately, grain yield. Stg and root angle QTL are also co-located and, together with crop water use data, suggest the role of roots in the stay-green phenomenon. Candidate genes have been identified in Stg1-4, including genes from the PIN family of auxin efflux carriers in Stg1 and Stg2, with 10 of 11 PIN genes in sorghum co-locating with Stg QTL. Modified gene expression in some of these PIN candidates in the stay-green compared with the senescent types has been found in preliminary RNA expression profiling studies. Further proof-of-function studies are underway, including comparative genomics, SNP analysis to assess diversity at candidate genes, reverse genetics and transformation.

Fine mapping of qGW1, a major QTL for grain weight in sorghum

Key message We detected seven QTLs for 100-grain weight in sorghum using an F 2 population, and delimited qGW1 to a 101-kb region on the short arm of chromosome 1, which contained 13 putative genes. Abstract Sorghum is one of the most important cereal crops. Breeding high-yielding sorghum varieties will have a profound impact on global food security. Grain weight is an important component of grain yield. It is a quantitative trait controlled by multiple quantitative trait loci (QTLs); however, the genetic basis of grain weight in sorghum is not well understood. In the present study, using an F2 population derived from a cross between the grain sorghum variety SA2313 (Sorghum bicolor) and the Sudan-grass variety Hiro-1 (S. bicolor), we detected seven QTLs for 100-grain weight. One of them, qGW1, was detected consistently over 2 years and contributed between 20 and 40 % of the phenotypic variation across multiple genetic backgrounds. Using extreme recombinants from a fine-mapping F3 population, we delimited qGW1 to a 101-kb region on the short arm of chromosome 1, containing 13 predicted gene models, one of which was found to be under purifying selection during domestication. However, none of the grain size candidate genes shared sequence similarity with previously cloned grain weight-related genes from rice. This study will facilitate isolation of the gene underlying qGW1 and advance our understanding of the regulatory mechanisms of grain weight. SSR markers linked to the qGW1 locus can be used for improving sorghum grain yield through marker-assisted selection.

Bayesian Inference for Retrospective Population Genetics Models Using Markov Chain Monte Carlo Methods

Genetics, the science of heredity and variation in living organisms, has a central role in medicine, in breeding crops and livestock, and in studying fundamental topics of biological sciences such as evolution and cell functioning. Currently the field of genetics is under a rapid development because of the recent advances in technologies by which molecular data can be obtained from living organisms. In order that most information from such data can be extracted, the analyses need to be carried out using statistical models that are tailored to take account of the particular genetic processes. In this thesis we formulate and analyze Bayesian models for genetic marker data of contemporary individuals. The major focus is on the modeling of the unobserved recent ancestry of the sampled individuals (say, for tens of generations or so), which is carried out by using explicit probabilistic reconstructions of the pedigree structures accompanied by the gene flows at the marker loci. For such a recent history, the recombination process is the major genetic force that shapes the genomes of the individuals, and it is included in the model by assuming that the recombination fractions between the adjacent markers are known. The posterior distribution of the unobserved history of the individuals is studied conditionally on the observed marker data by using a Markov chain Monte Carlo algorithm (MCMC). The example analyses consider estimation of the population structure, relatedness structure (both at the level of whole genomes as well as at each marker separately), and haplotype configurations. For situations where the pedigree structure is partially known, an algorithm to create an initial state for the MCMC algorithm is given. Furthermore, the thesis includes an extension of the model for the recent genetic history to situations where also a quantitative phenotype has been measured from the contemporary individuals. In that case the goal is to identify positions on the genome that affect the observed phenotypic values. This task is carried out within the Bayesian framework, where the number and the relative effects of the quantitative trait loci are treated as random variables whose posterior distribution is studied conditionally on the observed genetic and phenotypic data. In addition, the thesis contains an extension of a widely-used haplotyping method, the PHASE algorithm, to settings where genetic material from several individuals has been pooled together, and the allele frequencies of each pool are determined in a single genotyping.

Genomic regions influencing seminal root traits in barley

Water availability is a major limiting factor for crop production, making drought adaptation and its many component traits a desirable attribute of plant cultivars. Previous studies in cereal crops indicate that root traits expressed at early plant developmental stages, such as seminal root angle and root number, are associated with water extraction at different depths. Here, we conducted the first study to map seminal root traits in barley (Hordeum vulgare L.). Using a recently developed high-throughput phenotyping method, a panel of 30 barley genotypes and a doubled-haploid (DH) population (ND24260 × 'Flagship') comprising 330 lines genotyped with diversity array technology (DArT) markers were evaluated for seminal root angle (deviation from vertical) and root number under controlled environmental conditions. A high degree of phenotypic variation was observed in the panel of 30 genotypes: 13.5 to 82.2 and 3.6 to 6.9° for root angle and root number, respectively. A similar range was observed in the DH population: 16.4 to 70.5 and 3.6 to 6.5° for root angle and number, respectively. Seven quantitative trait loci (QTL) for seminal root traits (root angle, two QTL; root number, five QTL) were detected in the DH population. A major QTL influencing both root angle and root number (RAQ2/RNQ4) was positioned on chromosome 5HL. Across-species analysis identified 10 common genes underlying root trait QTL in barley, wheat (Triticum aestivum L.), and sorghum [Sorghum bicolor (L.) Moench]. Here, we provide insight into seminal root phenotypes and provide a first look at the genetics controlling these traits in barley.

Fine Mapping of qDor7, a Major QTL Affecting Seed Dormancy in Sorghum (Sorghum bicolor (L.) Moench)

Seed dormancy is a key domestication trait for major crops, which is acquired in long-term systems development processes and enables the survival of plants in adverse natural conditions. It is a complex trait under polygenic control and is affected by endogenous and environmental factors. In the present study, a major seed dormancy QTL in sorghum (Sorghum bicolor (L.) Moench), qDor7, detected previously, was fine mapped using a large, multi-generational population. The qDor7 locus was delimited to a 96-kb region which contains 16 predicted gene models. These results lay a solid foundation for cloning qDor7. In addition, the functional markers tightly linked to the seed dormancy QTL may be used in marker-assisted selection for seed dormancy in sorghum.

Directional dominance on stature and cognition in diverse human populations

Homozygosity has long been associated with rare, often devastating, Mendelian disorders1, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness2. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power3, 4. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one second, general cognitive ability and educational attainment (P < 1 × 10−300, 2.1 × 10−6, 2.5 × 10−10 and 1.8 × 10−10, respectively). In each case, increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months’ less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing evidence that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples5, 6, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection7, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been.

A preliminary genetic map of Zhikong scallop (Chlamys farreri Jones et Preston 1904)

Zhikong scallop (Chlamys farreri Jones et Preston 1904) is one of the most important aquaculture species in China. The development of a genetic linkage map would provide a powerful tool for the genetic improvement of this species. Amplified fragment length polymorphism (AFLP) is a PCR-based technique that has proven to be powerful in genome fingerprinting and mapping, and population analysis. Genetic maps of C. farreri were constructed using AFLP markers and a full-sib family with 60 progeny. A total of 503 segregating AFLP markers were obtained, with 472 following the Mendelian segregation ratio of 1:1 and 31 markers showing significant (P< 0.05) segregation distortion. The male map contained 166 informative AFLP markers in 23 linkage groups covering 2468 cM. The average distance between markers was 14.9 cM. The female genetic map consisted of 198 markers in 25 linkage groups spanning 3130 cM with an average inter-marker spacing of 15.8 cM. DNA polymorphisms that segregated in a 3:1 ratio as well as the AFLP markers that were heterozygous in both parents were included to construct combined linkage genetic map. Five shared linkage groups, ranging from 61.1 to 162.5 cM, were identified between the male and female maps, covering 431 cM. Amplified fragment length polymorphism markers appeared to be evenly distributed within the linkage groups. Although preliminary, these maps provide a starting point for the mapping of the functional genes and quantitative trait loci in C. farreri.

A genetic linkage map of Pacific white shrimp (Litopenaeus vannamei): sex-linked microsatellite markers and high recombination rates

Pacific white shrimp (Litopenaeus vannamei) is the leading species farmed in the Western Hemisphere and an economically important aquaculture species in China. In this project, a genetic linkage map was constructed using amplified fragment length polymorphism (AFLP) and microsatellite markers. One hundred and eight select AFLP primer combinations and 30 polymorphic microsatellite markers produced 2071 markers that were polymorphic in either of the parents and segregated in the progeny. Of these segregating markers, 319 were mapped to 45 linkage groups of the female framework map, covering a total of 4134.4 cM; and 267 markers were assigned to 45 linkage groups of the male map, covering a total of 3220.9 cM. High recombination rates were found in both parental maps. A sex-linked microsatellite marker was mapped on the female map with 6.6 cM to sex and a LOD of 17.8, two other microsatellite markers were also linked with both 8.6 cM to sex and LOD score of 14.3 and 16.4. The genetic maps presented here will serve as a basis for the construction of a high-resolution genetic map, quantitative trait loci (QTLs) detection, marker-assisted selection (MAS) and comparative genome mapping.

Epistasis among Drosophila persimilis factors conferring hybrid male sterility with D. pseudoobscura bogotana.

The Bateson-Dobzhansky-Muller model posits that hybrid incompatibilities result from genetic changes that accumulate during population divergence. Indeed, much effort in recent years has been devoted to identifying genes associated with hybrid incompatibilities, often with limited success, suggesting that hybrid sterility and inviability are frequently caused by complex interactions between multiple loci and not by single or a small number of gene pairs. Our previous study showed that the nature of epistasis between sterility-conferring QTL in the Drosophila persimilis-D. pseudoobscura bogotana species pair is highly specific. Here, we further dissect one of the three QTL underlying hybrid male sterility between these species and provide evidence for multiple factors within this QTL. This result indicates that the number of loci thought to contribute to hybrid dysfunction may have been underestimated, and we discuss how linkage and complex epistasis may be characteristic of the genetics of hybrid incompatibilities. We further pinpoint the location of one locus that confers hybrid male sterility when homozygous, dubbed "mule-like", to roughly 250 kilobases.

Two genes on A/J chromosome 18 are associated with susceptibility to Staphylococcus aureus infection by combined microarray and QTL analyses.

Although it has recently been shown that A/J mice are highly susceptible to Staphylococcus aureus sepsis as compared to C57BL/6J, the specific genes responsible for this differential phenotype are unknown. Using chromosome substitution strains (CSS), we found that loci on chromosomes 8, 11, and 18 influence susceptibility to S. aureus sepsis in A/J mice. We then used two candidate gene selection strategies to identify genes on these three chromosomes associated with S. aureus susceptibility, and targeted genes identified by both gene selection strategies. First, we used whole genome transcription profiling to identify 191 (56 on chr. 8, 100 on chr. 11, and 35 on chr. 18) genes on our three chromosomes of interest that are differentially expressed between S. aureus-infected A/J and C57BL/6J. Second, we identified two significant quantitative trait loci (QTL) for survival post-infection on chr. 18 using N(2) backcross mice (F(1) [C18A]xC57BL/6J). Ten genes on chr. 18 (March3, Cep120, Chmp1b, Dcp2, Dtwd2, Isoc1, Lman1, Spire1, Tnfaip8, and Seh1l) mapped to the two significant QTL regions and were also identified by the expression array selection strategy. Using real-time PCR, 6 of these 10 genes (Chmp1b, Dtwd2, Isoc1, Lman1, Tnfaip8, and Seh1l) showed significantly different expression levels between S. aureus-infected A/J and C57BL/6J. For two (Tnfaip8 and Seh1l) of these 6 genes, siRNA-mediated knockdown of gene expression in S. aureus-challenged RAW264.7 macrophages induced significant changes in the cytokine response (IL-1 beta and GM-CSF) compared to negative controls. These cytokine response changes were consistent with those seen in S. aureus-challenged peritoneal macrophages from CSS 18 mice (which contain A/J chromosome 18 but are otherwise C57BL/6J), but not C57BL/6J mice. These findings suggest that two genes, Tnfaip8 and Seh1l, may contribute to susceptibility to S. aureus in A/J mice, and represent promising candidates for human genetic susceptibility studies.

Resolving candidate genes of mouse skeletal muscle QTL via RNA-Seq and expression network analyses

BACKGROUND:

We have recently identified a number of Quantitative Trait Loci (QTL) contributing to the 2-fold muscle weight difference between the LG/J and SM/J mouse strains and refined their confidence intervals. To facilitate nomination of the candidate genes responsible for these differences we examined the transcriptome of the tibialis anterior (TA) muscle of each strain by RNA-Seq.

13,726 genes were expressed in mouse skeletal muscle. Intersection of a set of 1061 differentially expressed transcripts with a mouse muscle Bayesian Network identified a coherent set of differentially expressed genes that we term the LG/J and SM/J Regulatory Network (LSRN). The integration of the QTL, transcriptome and the network analyses identified eight key drivers of the LSRN (Kdr, Plbd1, Mgp, Fah, Prss23, 2310014F06Rik, Grtp1, Stk10) residing within five QTL regions, which were either polymorphic or differentially expressed between the two strains and are strong candidates for quantitative trait genes (QTGs) underlying muscle mass. The insight gained from network analysis including the ability to make testable predictions is illustrated by annotating the LSRN with knowledge-based signatures and showing that the SM/J state of the network corresponds to a more oxidative state. We validated this prediction by NADH tetrazolium reductase staining in the TA muscle revealing higher oxidative potential of the SM/J compared to the LG/J strain (p<0.03).

Thus, integration of fine resolution QTL mapping, RNA-Seq transcriptome information and mouse muscle Bayesian Network analysis provides a novel and unbiased strategy for nomination of muscle QTGs.

Natural variation of rice strigolactone biosynthesis is associated with the deletion of two MAX1 orthologs

Rice (Oryza sativa) cultivar Azucena--belonging to the Japonica subspecies--exudes high strigolactone (SL) levels and induces high germination of the root parasitic plant Striga hermonthica. Consistent with the fact that SLs also inhibit shoot branching, Azucena is a low-tillering variety. In contrast, Bala, an Indica cultivar, is a low-SL producer, stimulates less Striga germination, and is highly tillered. Using a Bala × Azucena F6 population, a major quantitative trait loci--qSLB1.1--for the exudation of SL, tillering, and induction of Striga germination was detected on chromosome 1. Sequence analysis of the corresponding locus revealed a rearrangement of a 51- to 59-kbp stretch between 28.9 and 29 Mbp in the Bala genome, resulting in the deletion of two cytochrome P450 genes--SLB1 and SLB2--with high homology to the Arabidopsis SL biosynthesis gene, MAX1. Both rice genes rescue the Arabidopsis max1-1 highly branched mutant phenotype and increase the production of the SL, ent-2'-epi-5-deoxystrigol, when overexpressed in Bala. Furthermore, analysis of this region in 367 cultivars of the publicly available Rice Diversity Panel population shows that the rearrangement at this locus is a recurrent natural trait associated with the Indica/Japonica divide in rice.