Molecular analysis of the parallel domestication of the common bean (Phaseolus vulgaris) in Mesoamerica and the Andes

We have studied the nucleotide diversity of common bean, Phaseolus vulgaris, which is characterized by two independent domestications in two geographically distinct areas: Mesoamerica and the Andes. This provides an important model, as domestication can be studied as a replicate experiment. We used nucleotide data from five gene fragments characterized by large introns to analyse 214 accessions (102 wild and 112 domesticated). The wild accessions represent a cross-section of the entire geographical distribution of P. vulgaris. A reduction in genetic diversity in both of these gene pools was found, which was three-fold greater in Mesoamerica compared with the Andes. This appears to be a result of a bottleneck that occurred before domestication in the Andes, which strongly impoverished this wild germplasm, leading to the minor effect of the subsequent domestication bottleneck (i.e. sequential bottleneck). These findings show the importance of considering the evolutionary history of crop species as a major factor that influences their current level and structure of genetic diversity. Furthermore, these data highlight a single domestication event within each gene pool. Although the findings should be interpreted with caution, this evidence indicates the Oaxaca valley in Mesoamerica, and southern Bolivia and northern Argentina in South America, as the origins of common bean domestication.

Genetic dissection of yield-related traits in a recombinant inbred line population created using a key breeding parent in China's wheat breeding

Understanding the genetics underlying yield formation of wheat is important for increasing wheat yield potential in breeding programs. Nanda2419 was a widely used cultivar for wheat production and breeding in China. In this study, we evaluated yield components and a few yield-related traits of a recombinant inbred line (RIL) population created by crossing Nanda2419 with the indigenous cultivar Wangshuibai in three to four trials at different geographical locations. Negative and positive correlations were found among some of these evaluated traits. Five traits had over 50 % trial-wide broad sense heritability. Using a framework marker map of the genome constructed with this population, quantitative trait loci (QTL) were identified for all traits, and epistatic loci were identified for seven of them. Our results confirmed some of the previously reported QTLs in wheat and identified several new ones, including QSn.nau-6D for effective tillers, QGn.nau-4B.2 for kernel number, QGw.nau-4D for kernel weight, QPh.nau-4B.2 and QPh.nau-4A for plant height, and QFlw.nau-5A.1 for flag leaf width. In the investigated population, Nanda2419 contributed all QTLs associated with higher kernel weight, higher leaf chlorophyll content, and a major QTL associated with wider flag leaf. Seven chromosome regions were related to more than one trait. Four QTL clusters contributed positively to breeding goal-based trait improvement through the Nanda2419 alleles and were detected in trials set in different ecological regions. The findings of this study are relevant to the molecular improvement of wheat yield and to the goal of screening cultivars for better breeding parents.

Development and Genetic Characterization of an Advanced Backcross-Nested Association Mapping (AB-NAM) Population of Wild × Cultivated Barley

The ability to access alleles from unadapted germplasm collections is a long-standing problem for geneticists and breeders. Here we developed, characterized, and demonstrated the utility of a wild barley advanced backcross-nested association mapping (AB-NAM) population. We developed this population by backcrossing 25 wild barley accessions to the six-rowed malting barley cultivar Rasmusson. The 25 wild barley parents were selected from the 318 accession Wild Barley Diversity Collection (WBDC) to maximize allelic diversity. The resulting 796 BC2F4:6 lines were genotyped with 384 SNP markers, and an additional 4022 SNPs and 263,531 sequence variants were imputed onto the population using 9K iSelect SNP genotypes and exome capture sequence of the parents, respectively. On average, 96% of each wild parent was introgressed into the Rasmusson background, and the population exhibited low population structure. While linkage disequilibrium (LD) decay (r(2) = 0.2) was lowest in the WBDC (0.36 cM), the AB-NAM (9.2 cM) exhibited more rapid LD decay than comparable advanced backcross (28.6 cM) and recombinant inbred line (32.3 cM) populations. Three qualitative traits: glossy spike, glossy sheath, and black hull color were mapped with high resolution to loci corresponding to known barley mutants for these traits. Additionally, a total of 10 QTL were identified for grain protein content. The combination of low LD, negligible population structure, and high diversity in an adapted background make the AB-NAM an important tool for high-resolution gene mapping and discovery of novel allelic variation using wild barley germplasm.

Genome-wide estimates of coancestry and inbreeding in a closed herd of ancient Iberian pigs

Maintaining genetic variation and controlling the increase in inbreeding are crucial requirements in animal conservation programs. The most widely accepted strategy for achieving these objectives is to maximize the effective population size by minimizing the global coancestry obtained from a particular pedigree. However, for most natural or captive populations genealogical information is absent. In this situation, microsatellites have been traditionally the markers of choice to characterize genetic variation, and several estimators of genealogical coefficients have been developed using marker data, with unsatisfactory results. The development of high-throughput genotyping techniques states the necessity of reviewing the paradigm that genealogical coancestry is the best parameter for measuring genetic diversity. In this study, the Illumina PorcineSNP60 BeadChip was used to obtain genome-wide estimates of rates of coancestry and inbreeding and effective population size for an ancient strain of Iberian pigs that is now in serious danger of extinction and for which very accurate genealogical information is available (the Guadyerbas strain). Genome-wide estimates were compared with those obtained from microsatellite and from pedigree data. Estimates of coancestry and inbreeding computed from the SNP chip were strongly correlated with genealogical estimates and these correlations were substantially higher than those between microsatellite and genealogical coefficients. Also, molecular coancestry computed from SNP information was a better predictor of genealogical coancestry than coancestry computed from microsatellites. Rates of change in coancestry and inbreeding and effective population size estimated from molecular data were very similar to those estimated from genealogical data. However, estimates of effective population size obtained from changes in coancestry or inbreeding differed. Our results indicate that genome-wide information represents a useful alternative to genealogical information for measuring and maintaining genetic diversity.

Inbreeding and selection shape genomic diversity in captive populations: Implications for the conservation of endangered species

Captive breeding programs are often initiated to prevent species extinction until reintroduction into the wild can occur. However, the evolution of captive populations via inbreeding, drift, and selection can impair fitness, compromising reintroduction programs. To better understand the evolutionary response of species bred in captivity, we used nearly 5500 single nucleotide polymorphisms (SNPs) in populations of white-footed mice (Peromyscus leucopus) to measure the impact of breeding regimes on genomic diversity. We bred mice in captivity for 20 generations using two replicates of three protocols: random mating (RAN), selection for docile behaviors (DOC), and minimizing mean kinship (MK). The MK protocol most effectively retained genomic diversity and reduced the effects of selection. Additionally, genomic diversity was significantly related to fitness, as assessed with pedigrees and SNPs supported with genomic sequence data. Because captive-born individuals are often less fit in wild settings compared to wild-born individuals, captive-estimated fitness correlations likely underestimate the effects in wild populations. Therefore, minimizing inbreeding and selection in captive populations is critical to increasing the probability of releasing fit individuals into the wild.

Allele frequency changes due to hitch-hiking in genomic selection programs

BACKGROUND

Genomic selection makes it possible to reduce pedigree-based inbreeding over best linear unbiased prediction (BLUP) by increasing emphasis on own rather than family information. However, pedigree inbreeding might not accurately reflect loss of genetic variation and the true level of inbreeding due to changes in allele frequencies and hitch-hiking. This study aimed at understanding the impact of using long-term genomic selection on changes in allele frequencies, genetic variation and level of inbreeding.

METHODS

Selection was performed in simulated scenarios with a population of 400 animals for 25 consecutive generations. Six genetic models were considered with different heritabilities and numbers of QTL (quantitative trait loci) affecting the trait. Four selection criteria were used, including selection on own phenotype and on estimated breeding values (EBV) derived using phenotype-BLUP, genomic BLUP and Bayesian Lasso. Changes in allele frequencies at QTL, markers and linked neutral loci were investigated for the different selection criteria and different scenarios, along with the loss of favourable alleles and the rate of inbreeding measured by pedigree and runs of homozygosity.

RESULTS

For each selection criterion, hitch-hiking in the vicinity of the QTL appeared more extensive when accuracy of selection was higher and the number of QTL was lower. When inbreeding was measured by pedigree information, selection on genomic BLUP EBV resulted in lower levels of inbreeding than selection on phenotype BLUP EBV, but this did not always apply when inbreeding was measured by runs of homozygosity. Compared to genomic BLUP, selection on EBV from Bayesian Lasso led to less genetic drift, reduced loss of favourable alleles and more effectively controlled the rate of both pedigree and genomic inbreeding in all simulated scenarios. In addition, selection on EBV from Bayesian Lasso showed a higher selection differential for mendelian sampling terms than selection on genomic BLUP EBV.

CONCLUSIONS

Neutral variation can be shaped to a great extent by the hitch-hiking effects associated with selection, rather than just by genetic drift. When implementing long-term genomic selection, strategies for genomic control of inbreeding are essential, due to a considerable hitch-hiking effect, regardless of the method that is used for prediction of EBV.

The value of genomic relationship matrices to estimate levels of inbreeding

BACKGROUND

Genomic relationship matrices are used to obtain genomic inbreeding coefficients. However, there are several methodologies to compute these matrices and there is still an unresolved debate on which one provides the best estimate of inbreeding. In this study, we investigated measures of inbreeding obtained from five genomic matrices, including the Nejati-Javaremi allelic relationship matrix (FNEJ), the Li and Horvitz matrix based on excess of homozygosity (FL&H), and the VanRaden (methods 1, FVR1, and 2, FVR2) and Yang (FYAN) genomic relationship matrices. We derived expectations for each inbreeding coefficient, assuming a single locus model, and used these expectations to explain the patterns of the coefficients that were computed from thousands of single nucleotide polymorphism genotypes in a population of Iberian pigs.

RESULTS

Except for FNEJ, the evaluated measures of inbreeding do not match with the original definitions of inbreeding coefficient of Wright (correlation) or Malécot (probability). When inbreeding coefficients are interpreted as indicators of variability (heterozygosity) that was gained or lost relative to a base population, both FNEJ and FL&H led to sensible results but this was not the case for FVR1, FVR2 and FYAN. When variability has increased relative to the base, FVR1, FVR2 and FYAN can indicate that it decreased. In fact, based on FYAN, variability is not expected to increase. When variability has decreased, FVR1 and FVR2 can indicate that it has increased. Finally, these three coefficients can indicate that more variability than that present in the base population can be lost, which is also unreasonable. The patterns for these coefficients observed in the pig population were very different, following the derived expectations. As a consequence, the rate of inbreeding depression estimated based on these inbreeding coefficients differed not only in magnitude but also in sign.

CONCLUSIONS

Genomic inbreeding coefficients obtained from the diagonal elements of genomic matrices can lead to inconsistent results in terms of gain and loss of genetic variability and inbreeding depression estimates, and thus to misleading interpretations. Although these matrices have proven to be very efficient in increasing the accuracy of genomic predictions, they do not always provide a useful measure of inbreeding.

What happens after inbreeding avoidance? Inbreeding by rejected relatives and the inclusive fitness benefit of inbreeding avoidance

Avoiding inbreeding, and therefore avoiding inbreeding depression in offspring fitness, is widely assumed to be adaptive in systems with biparental reproduction. However, inbreeding can also confer an inclusive fitness benefit stemming from increased relatedness between parents and inbred offspring. Whether or not inbreeding or avoiding inbreeding is adaptive therefore depends on a balance between inbreeding depression and increased parent-offspring relatedness. Existing models of biparental inbreeding predict threshold values of inbreeding depression above which males and females should avoid inbreeding, and predict sexual conflict over inbreeding because these thresholds diverge. However, these models implicitly assume that if a focal individual avoids inbreeding, then both it and its rejected relative will subsequently outbreed. We show that relaxing this assumption of reciprocal outbreeding, and the assumption that focal individuals are themselves outbred, can substantially alter the predicted thresholds for inbreeding avoidance for focal males. Specifically, the magnitude of inbreeding depression below which inbreeding increases a focal male’s inclusive fitness increases with increasing depression in the offspring of a focal female and her alternative mate, and it decreases with increasing relatedness between a focal male and a focal female’s alternative mate, thereby altering the predicted zone of sexual conflict. Furthermore, a focal male’s inclusive fitness gain from avoiding inbreeding is reduced by indirect opportunity costs if his rejected relative breeds with another relative of his. By demonstrating that variation in relatedness and inbreeding can affect intra- and inter-sexual conflict over inbreeding, our models lead to novel predictions for family dynamics. Specifically, parent-offspring conflict over inbreeding might depend on the alternative mates of rejected relatives, and male-male competition over inbreeding might lead to mixed inbreeding strategies. Making testable quantitative predictions regarding inbreeding strategies occurring in nature will therefore require new models that explicitly capture variation in relatedness and inbreeding among interacting population members.

Introgression of null allele of Kunitz trypsin inhibitor through marker-assisted backcross breeding in soybean (Glycine max L. Merr.)

BACKGROUND

Presence of Kunitz trypsin inhibitor (KTI) in soybean seeds necessitates pre-heat treatment of the soy-flour for its inactivation before using it in food and feed products. The heat treatment not only enhances processing costs of the soy-based foods and feeds but also affects seed-protein quality and solubility. Genetic elimination of KTI is an important and effective approach. Therefore, molecular marker-assisted backcross breeding (MABB) approach was adopted for genetic elimination of KTI from two popular soybean genotypes, DS9712 and DS9814. PI542044, an exotic germplasm line was used as donor of the kti allele which inhibits functional KTI peptide production.

RESULTS

Foreground selection for the kti allele was performed with three closely linked SSR markers while background selection was done with 93 polymorphic SSR markers. Plants in the BC1F1 generation were found to recover 70.4-87.63 % and 60.26-73.78 % of the recurrent parent genome (RPG) of DS9712 and DS9814, respectively. Similarly, selected plants in the BC2F1 generation had 93.01-98.92 % and 83.3-91.67 % recovery of their respective RPGs. Recombinant selection was performed so as to identify plants with minimal linkage drag. Biochemical analysis of the seeds of the selected plants (ktikti) confirmed absence of KTI peptides in the seeds. Phenotypically, the selected plants were comparable to the respective recurrent parent in yield and other traits.

CONCLUSIONS

MABB approach helped in speedy development of 6 KTI free breeding lines of soybean. Such lines will be suitable for the farmers and the soybean industries to use in production of soy-based foods and feeds without pre-heat treatment of the soy-flour. It would contribute towards wider acceptability of soy-based foods and feeds.

Genetic diversity and population history of eight Italian beef cattle breeds using measures of autozygosity

In the present study, GeneSeek GGP-LDv4 33k single nucleotide polymorphism chip was used to detect runs of homozygosity (ROH) in eight Italian beef cattle breeds, six breeds with distribution limited to Tuscany (Calvana, Mucca Pisana, Pontremolese) or Sardinia (Sarda, Sardo Bruna and Sardo Modicana) and two cosmopolitan breeds (Charolais and Limousine). ROH detection analyses were used to estimate autozygosity and inbreeding and to identify genomic regions with high frequency of ROH, which might reflect selection signatures. Comparative analysis among breeds revealed differences in length and distribution of ROH and inbreeding levels. The Charolais, Limousine, Sarda, and Sardo Bruna breeds were found to have a high frequency of short ROH (~ 15.000); Calvana and Mucca Pisana presented also runs longer than 16 Mbp. The highest level of average genomic inbreeding was observed in Tuscan breeds, around 0.3, while Sardinian and cosmopolitan breeds showed values around 0.2. The population structure and genetic distances were analyzed through principal component and multidimensional scaling analyses, and resulted in a clear separation among the breeds, with clusters related to productive purposes. The frequency of ROH occurrence revealed eight breed-specific genomic regions where genes of potential selective and conservative interest are located (e.g. MYOG, CHI3L1, CHIT1 (BTA16), TIMELESS, APOF, OR10P1, OR6C4, OR2AP1, OR6C2, OR6C68, CACNG2 (BTA5), COL5A2 and COL3A1 (BTA2)). In all breeds, we found the largest proportion of homozygous by descent segments to be those that represent inbreeding events that occurred around 32 generations ago, with Tuscan breeds also having a significant proportion of segments relating to more recent inbreeding.

Development of sub-tropically adapted diverse provitamin-A rich maize inbreds through marker-assisted pedigree selection, their characterization and utilization in hybrid breeding

Malnutrition has emerged as one of the major health problems worldwide. Traditional yellow maize has low provitamin-A (proA) content and its genetic base in proA biofortification breeding program of subtropics is extremely narrow. To diversify the proA rich germplasm, 10 elite low proA inbreds were crossed with a proA rich donor (HP702-22) having mutant crtRB1 gene. The F2 populations derived from these crosses were genotyped using InDel marker specific to crtRB1. Severe marker segregation distortion was observed. Seventeen crtRB1 inbreds developed through marker-assisted pedigree breeding and seven inbreds generated using marker-assisted backcross breeding were characterized using 77 SSRs. Wide variation in gene diversity (0.08 to 0.79) and dissimilarity coefficient (0.28 to 0.84) was observed. The inbreds were grouped into three major clusters depicting the existing genetic diversity. The crtRB1-based inbreds possessed high β-carotene (BC: 8.72μg/g), β-cryptoxanthin (BCX: 4.58μg/g) and proA (11.01μg/g), while it was 2.35μg/g, 1.24μg/g and 2.97μg/g in checks, respectively. Based on their genetic relationships, 15 newly developed crtRB1-based inbreds were crossed with five testers (having crtRB1 gene) using line × tester mating design. 75 experimental hybrids with crtRB1 gene were evaluated over three locations. These experimental hybrids possessed higher BC (8.02μg/g), BCX (4.69μg/g), proA (10.37μg/g) compared to traditional hybrids used as check (BC: 2.36 μg/g, BCX: 1.53μg/g, proA: 3.13μg/g). Environment and genotypes × environment interaction had minor effects on proA content. Both additive and dominance gene action were significant for proA. The mean proportion of proA to total carotenoids (TC) was 44% among crtRB1-based hybrids, while 11% in traditional hybrids. BC was found to be positively correlated with BCX (r = 0.68) and proA (r = 0.98). However, no correlation was observed between proA and grain yield. Several hybrids with >10.0 t/ha grain yield with proA content >10.0 μg/g were identified. This is the first comprehensive study on development of diverse proA rich maize hybrids through marker-assisted pedigree breeding approach. The findings provides sustainable and cost-effective solution to alleviate vitamin-A deficiency.

Suitability of Pedigree Information and Genomic Methods for Analyzing Inbreeding of Polish Cold-Blooded Horses Covered by Conservation Programs

Traditionally, pedigree-based relationship coefficients were used to manage inbreeding and control inbreeding depression that occurs within populations. The extensive incorporation of genomic data in livestock breeding creates the opportunity to develop and implement methods to manage populations at the genomic level. Consequently, the realized proportion of the genome that two individuals share can be more accurately estimated instead of using pedigree information to estimate the expected proportion of shared alleles. To make use of this improvement, in this study we evaluated the genomic inbreeding measures in the Polish conserved cold-blooded horse population and compared the data with the traditional measures of inbreeding. Additionally, an ancestry fractions/proportions from Admixture software were tested as an estimate of lineage (ancestry coefficient) used for horses qualifying for the conservation program. The highest correlation of pedigree-based (F_PED) and genomic inbreeding estimates was found for F_ROH (runs of homozygosity-based F coefficient) and F_UNI (F coefficient based on the correlation between uniting gametes). F_ROH correlation with F_PED tended to increase as the number of generations registered as pedigree increased. While lineage and gene contributions (Q) from Admixture software correlated, they showed poor direct compliance; hence, Q-value cannot be recommended as the estimate of pedigree-based lineage. All these findings suggest that the methods of genomics should be considered as an alternative or support in the analysis of population structure in conservative breeding that can help control inbreeding in rare horse populations.

Laboratory rearing of Anopheles arabiensis: impact on genetic variability and implications for Sterile Insect Technique (SIT) based mosquito control in northern Sudan

BACKGROUND

Mosquito colony populations often show significant changes in their population genetic make-up compared to the field populations that were used as founding source. Most of the changes that have been reported are indicators of depletion in the overall genetic diversity of the colony populations. The Sterile Insect Techniques programme of mosquito control that is underway in Northern Sudan uses sterilized males produced from a laboratory-maintained colony population. The genetic diversity of an advanced generation of this colony population was quantitatively assessed and compared to the field population from which the colony was derived.

METHODS

Anopheles arabiensis mosquito samples from the 13th generation of the colony, and from the locality that was the source of the first generation of the colony, were genotyped at 11 microsatellite loci distributed throughout the species' genome. Standard population genetic analyses were carried out to quantify and compare their population genetic make-up and diversities.

RESULTS

The colony samples showed significant reduction in the total number of alleles, the numbers of rare and private alleles, and the fractions of heterozygote individuals at all the loci. The pattern of change is consistent with the expected effect of the use of a small number of mosquitoes when the colony was established. Departure from Hardy-Weinberg equilibrium in the direction of homozygote excess was observed at some loci and attributed to the presence of null-alleles.

CONCLUSIONS

This study highlights the need for broad sampling when initiating colony populations and for ongoing assessment of the population genetic make-up of colony populations. Previous assessments of survivorship, dispersive behaviour and swarm formation indicate that the inbreeding and reduced genetic variability reported in this study may not have had direct fitness consequences yet. However, noting the lessons learned in other SIT programmes about the impact of colonization on male sexual behaviour and longevity, as well as other inbreeding related adverse effects, a systematic investigation of these potential effects is recommended because they have direct impact on the ultimate success of the programme.

Development of chickpea near-isogenic lines for Fusarium wilt

Four pairs of near-isogenic lines (NILs) of chickpea with resistance/susceptibility to Fusarium oxysporum f. sp. ciceris (Foc) have been developed in this study. These lines were produced by searching in advanced recombinant inbred lines (RILs) that are segregating for Foc race 5 based on a phenotypic screening. The sequence tagged microsatellite (STMS) marker TA59, closely linked to wilt resistance genes on linkage group 2 (LG2) of the chickpea map, was used to assist the selection of resistant or susceptible genotypes. The NILs were also characterized for disease reaction to Foc races 1A, 2, 3 and 4. Resistance, susceptibility and slow wilting reactions were found in these NILs. Our results suggest that more than one gene controls the resistance to race 5. Combination of the major gene foc-5 linked to TA59 with other gene/s appears to be required to complete resistance, and the absence of these unknown genes leads to slow wilting reactions. The independent differential responses to races 2 and 3 observed in three NILs could be explained as recombination events. This result suggests that foc-2 and foc-3 are delimiting points at opposite ends of a genomic region that includes the remaining foc genes and the TA59 marker. This set of NILs has great potential for studying the genetics and mechanisms of wilt resistance. In addition, the NIL RIP8-94-11 can be used as differential line for Foc race 3; it showed a clear resistance reaction to race 3 and susceptibility to the other Foc races.

Increased rate of close-kin unions in the central Andes in the half millennium before European contact

Analyzing ancient DNA of the central Andes, Ringbauer and colleagues identify a markedly elevated rate of unions of closely related parents after ca. 1000 CE. This change of mating preferences sheds new light on a unique system of social organization based on ancestry ("ayllu") whereby within-group unions were preferred to facilitate sharing of resources.

Assessment of the Genetic Diversity of a Local Pig Breed Using Pedigree and SNP Data

Herein, the genetic diversity of the local Přeštice Black-Pied pig breed was assessed by the simultaneous analysis of the pedigree and single nucleotide polymorphism (SNP) data. The information about sire line, dam, date of birth, sex, breeding line, and herd for 1971 individuals was considered in the pedigree analysis. The SNP analysis (n = 181) was performed using the Illumina PorcineSNP60 BeadChip kit. The quality of pedigree and SNPs and the inbreeding coefficients (F) and effective population size (Ne) were evaluated. The correlations between inbreeding based on the runs of homozygosity (FROH) and pedigree (FPED) were also calculated. The average FPED for all animals was 3.44%, while the FROH varied from 10.81% for a minimum size of 1 Mbp to 3.98% for a minimum size of 16 Mbp. The average minor allele frequency was 0.28 ± 0.11. The observed and expected within breed heterozygosities were 0.38 ± 0.13 and 0.37 ± 0.12, respectively. The Ne, obtained using both the data sources, reached values around 50 animals. Moderate correlation coefficients (0.49-0.54) were observed between FPED and FROH. It is necessary to make decisions that stabilize the inbreeding rate in the long-term using optimal contribution selection based on the available SNP data.

High-Speed Mouse Backcrossing Through the Female Germ Line

Transferring mouse mutations into specific mouse strain backgrounds can be critical for appropriate analysis of phenotypic effects of targeted genomic alterations and quantitative trait loci. Speed congenic breeding strategies incorporating marker-assisted selection of progeny with the highest percentage target background as breeders for the next generation can produce congenic strains within approximately 5 generations. When mating selected donor males to target strain females, this may require more than 1 year, with each generation lasting 10 to 11 weeks including 3 weeks of gestation and 7 to 8 weeks until the males reach sexual maturity. Because ovulation can be induced in female mice as early as 3 weeks of age, superovulation-aided backcrossing of marker-selected females could accelerate the production of congenic animals by approximately 4 weeks per generation, reducing time and cost. Using this approach, we transferred a transgenic strain of undefined genetic background to >99% C57BL/6J within 10 months, with most generations lasting 7 weeks. This involved less than 60 mice in total, with 9 to 18 animals per generation. Our data demonstrate that high-speed backcrossing through the female germline is feasible and practical with small mouse numbers.

Introgression of blast resistance genes into the elite rice variety MR263 through marker-assisted backcrossing

BACKGROUND

Blast caused by the fungus Magnaporthe oryzae is a significant disease threat to rice across the world and is especially prevalent in Malaysia. An elite, early-maturing, high-yielding Malaysian rice variety, MR263, is susceptible to blast and was used as the recurrent parent in this study. To improve MR263 disease resistance, the Pongsu Seribu 1 rice variety was used as donor of the blast resistance Pi-7(t), Pi-d(t)1 and Pir2-3(t) genes and qLN2 quantitative trait locus (QTL). The objective was to introgress these blast resistance genes into the background of MR263 using marker-assisted backcrossing with both foreground and background selection.

RESULTS

Improved MR263-BR-3-2, MR263-BR-4-3, MR263-BR-13-1 and MR263-BR-26-4 lines carrying the Pi-7(t), Pi-d(t)1 and Pir2-3(t) genes and qLN2 QTL were developed using the simple sequence repeat (SSR) markers RM5961 and RM263 (linked to the blast resistance genes and QTL) for foreground selection and a collection of 65 polymorphic SSR markers for background selection in backcrossed and selfed generations. A background analysis revealed that the highest rate of recurrent parent genome recovery was 96.1% in MR263-BR-4-3 and 94.3% in MR263-BR-3-2.

CONCLUSION

The addition of blast resistance genes can be used to improve several Malaysian rice varieties to combat this major disease.

Fertile introgression products generated via somatic hybridization between wheat and Thinopyrum intermedium

Fertile hybrids were produced with genetic material transferred from Th. intermedium into a wheat background and supply a source of genetic variation to wheat improvement. Both symmetric and asymmetric somatic hybrids have been obtained from the combination of wheatgrass (Thinopyrum intermedium) and bread wheat (Triticum aestivum). Two wheat protoplast populations, one derived from embryogenic calli and the other from a non-regenerable, rapidly dividing cell line, were fused with Th. intermedium protoplasts which had been (or not been) pre-irradiated with UV. Among the 124 regenerated calli, 64 could be categorized as being of hybrid origin on the basis of plant morphology, peroxidase isozyme, RAPD DNA profiling and karyological analysis. Numerous green plantlets were regenerated from 13 calli recovered from either the symmetric hybrid (no UV pre-treatment) or the asymmetric one (30 s UV irradiation). One of these hybrid plants proved to be vigorous and self-fertile. The regenerants were all closer in phenotype to wheat than to Th. intermedium. Genomic in situ hybridization analysis showed that the chromosomes in the hybrids were largely intact wheat ones, although a few Th. intermedium chromosome fragments had been incorporated within them.

Genetic diversity and the application of runs of homozygosity-based methods for inbreeding estimation in German White-headed Mutton sheep

The German White-headed Mutton (GWM) sheep is a monitoring population believed to have been improved through crosses with other breeds, e.g., Texel (TXL) and French Berrichone du Cher (BDC). The primary aim of the study was to analyse genetic diversity and breed composition of GWM sheep. Furthermore, different measures of computing inbreeding from the runs of homozygosity (ROH) were investigated. Data for GWM consisted of pedigree information on 19,000 animals and 40,753 quality filtered SNPs on 46 individuals. Additionally, publicly available genotype data on 209 individuals belonging to nine sheep breeds were included in the analysis. Due to evenness of SNPs spacing and proportionality of the number of SNPs in each autosome to autosome length, a high correlation (rp = 0.99) was found between genomic inbreeding coefficients computed based on the length of ROH (FROH_L) and those computed relative to the number of SNPs in ROH (FROH_N). Total inbreeding was partitioned into values for individual chromosomes revealing the highest levels of inbreeding on chromosomes 1, 2 and 3. Correlations between the ROH-based inbreeding measures and pedigree inbreeding reached 0.82. The observed heterozygosity estimate in GWM was high (0.39), however, the breed suffered low level of effective population size (~50) from a genomic viewpoint. Moreover, effective number of founders (186), and effective number of ancestors (144) implied disequilibrium of founder contribution and a genetic bottleneck in the breed. Multidimensional scaling and network visualisation analyses revealed close connectedness of GWM to BDC and German Texel (GTX). A model-based admixture analysis consistently indicated the flow of genes from other breeds, particularly BDC to GWM. Our analyses highlight the mixed genetic background of GWM sheep and furthermore, suggest a close monitoring of the breed to consolidate its genetic diversity while averting further reduction in the effective population size.

Improvement of a Yairipok Chujak Maize Landrace from North Eastern Himalayan Region for β-Carotene Content through Molecular Marker-Assisted Backcross Breeding

In the North Eastern Himalayan region (NEHR) of India, maize is an important food crop. The local people cultivate the maize landraces and consume them as food. However, these landraces are deficient in β-carotene content. Thus, we aimed to incorporate the crtRB1 gene from UMI285β⁺ into the genetic background of the NEHR maize landrace, Yairipok Chujak (CAUM66), and thereby enhance the β-carotene content through marker-assisted backcrossing (MABC). In this regard, we backcrossed and screened BC₁F₁ and BC₂F₁ plants possessing the heterozygous allele for crtRB1 and then screened with 106 polymorphic simple sequence repeat (SSR) markers. The plants having maximum recurrent parent genome recovery (RPGR) were selected in each generation and selfed to produce BC₂F₂ seeds. In the BC₂F₂ generation, four plants (CAUM66-54-9-12-2, CAUM66-54-9-12-11, CAUM66-54-9-12-13, and CAUM66-54-9-12-24) having homozygous crtRB1-favorable allele with maximum RPGR (86.74–90.16%) were selected and advanced to BC₂F₃. The four selected plants were selfed to produce BC₂F₃ and then evaluated for agronomic traits and β-carotene content. The agronomic performance of the four lines was similar (78.83–99.44%) to that of the recurrent parent, and β-carotene content (7.541–8.711 μg/g) was on par with the donor parent. Our study is the first to improve the β-carotene content in NEHR maize landrace through MABC. The newly developed lines could serve as potential resources to further develop nutrition-rich maize lines and could provide genetic stock for use in breeding programs.

Estimation of genomic inbreeding coefficients based on high-density SNP markers in Chinese Holstein cattle

In livestock, inbreeding coefficient based on pedigree information is usually used to evaluate the level of inbreeding. Recently, with cost reduction of high-density SNP genotyping, it's possible to analyze real genomic inbreeding degree using genomic information. In this study, utilizing high-density SNP chip data, we analyzed the frequency and distribution of runs of homozygosity (ROH) in 2107 Chinese Holstein cattle in Beijing area, and calculated 2 genomic inbreeding coefficients, i.e., 1) the proportion of ROH length in the total length of autosomal genome (Froh), and 2) the percentage of homozygous SNPs (Fhom). Then we analyzed the correlation between 2 genomic inbreeding coefficients and the correlation between genomic and pedigree inbreeding coefficients. We totally detected 44 676 ROHs that mainly ranged from 1 to 10 Mb. Various lengths of ROHs existed in the genome. There were more short ROHs than long ROHs. ROHs aren't evenly distributed in chromosomes. The area with most ROHs is in the middle part of chromosome 10. Strong correlation (r > 0.90) existed between 2 kinds of genomic inbreeding coefficients, but the correlation between pedigree and genomic inbreeding coefficients were much lower (r < 0.50). Our finding suggests that pedigree completeness influences the correlation between genomic and pedigree inbreeding. Genomic inbreeding measures may reflect individuals' real inbreeding, which could be a useful tool to evaluate population inbreeding.

Genotyping both live and dead animals to improve post-weaning survival of pigs in breeding programs

BACKGROUND

In this study, we tested whether genotyping both live and dead animals (GSD) realises more genetic gain for post-weaning survival (PWS) in pigs compared to genotyping only live animals (GOS).

METHODS

Stochastic simulation was used to estimate the rate of genetic gain realised by GSD and GOS at a 0.01 rate of pedigree-based inbreeding in three breeding schemes, which differed in PWS (95%, 90% and 50%) and litter size (6 and 10). Pedigree-based selection was conducted as a point of reference. Variance components were estimated and then estimated breeding values (EBV) were obtained in each breeding scheme using a linear or a threshold model. Selection was for a single trait, i.e. PWS with a heritability of 0.02 on the observed scale. The trait was simulated on the underlying scale and was recorded as binary (0/1). Selection candidates were genotyped and phenotyped before selection, with only live candidates eligible for selection. Genotyping strategies differed in the proportion of live and dead animals genotyped, but the phenotypes of all animals were used for predicting EBV of the selection candidates.

RESULTS

Based on a 0.01 rate of pedigree-based inbreeding, GSD realised 14 to 33% more genetic gain than GOS for all breeding schemes depending on PWS and litter size. GSD increased the prediction accuracy of EBV for PWS by at least 14% compared to GOS. The use of a linear versus a threshold model did not have an impact on genetic gain for PWS regardless of the genotyping strategy and the bias of the EBV did not differ significantly among genotyping strategies.

CONCLUSIONS

Genotyping both dead and live animals was more informative than genotyping only live animals to predict the EBV for PWS of selection candidates, but with marginal increases in genetic gain when the proportion of dead animals genotyped was 60% or greater. Therefore, it would be worthwhile to use genomic information on both live and more than 20% dead animals to compute EBV for the genetic improvement of PWS under the assumption that dead animals reflect increased liability on the underlying scale.

Genetic mapping for grain quality and yield-attributed traits in Basmati rice using SSR-based genetic map

Rice grain shape and nutritional quality traits have high economic value for commercial production of rice and largely determine the market price, besides influencing the global food demand for high-quality rice. Detection, mapping and exploitation of quantitative trait loci (QTL) associated with kernel elongation and grain quality in Basmati rice is considered as an efficient strategy for improving the kernel elongation and grain quality trait in rice varieties. Genetic information in rice for most of these traits is scanty and needed interventions through the use of molecular markers. A recombinant inbred lines (RIL) population consisting of 130 lines generated from the cross involving Basmati 370, a superior quality Basmati variety and Pusa Basmati 1121, a Basmati derived variety were used to map the QTLs for 9 important grain quality and yield related traits. Correlation studies showed that various components of yield show a significant positive relationship with grain yield. A genetic map was constructed using 70 polymorphic simple sequence repeat (SSR) markers spanning a genetic distance of 689.3 cM distributed over 12 rice chromosomes. Significant variation was observed and showed transgressive segregation for grain quality traits in RIL population. A total of 20 QTLs were identified associated with nine yield and quality traits. Epistatic interactions were also identified for grain quality related traits indicating complex genetic nature inheritance. Therefore, the identified QTLs and flanking marker information could be utilized in the marker-assisted selection to improve kernel elongation and nutritional grain quality traits in rice varieties.

Genomic Prediction of Individual Inbreeding Levels for the Management of Genetic Diversity in Populations With Small Effective Size

In populations of small effective size (Ne), such as those in conservation programmes, companion animals or livestock species, inbreeding control is essential. Homozygosity-by-descent (HBD) segments provide relevant information in that context, as they allow accurate estimation of the inbreeding coefficient, provide locus-specific information and their length is informative about the "age" of inbreeding. Our objective was to evaluate tools for predicting HBD in future offspring based on parental genotypes, a problem equivalent to identifying segments identical-by-descent (IBD) among the four parental chromosomes. In total, we reviewed and evaluated 16 approaches using simulated and real data from populations with small Ne. The methods included model-based approaches as well as more computationally efficient rule-based approaches. The accuracy of the methods was then evaluated, including with low-density marker panels, genotyping-by-sequencing data and small groups of individuals, typical features of such populations. Two model-based approaches performed consistently well, while some rule-based approaches proved accurate for genome-wide predictions. The model-based approaches were particularly efficient when genomic information was sparse or degraded. Methods using phased data proved to be more accurate, while some approaches relying on unphased genotype data were sensitive to the assumed allele frequencies. In some settings, pedigree-based predictions ranked high for recent inbreeding levels. Finally, we showed that our evaluation is also informative about the accuracy of the methods for estimating relatedness and identifying IBD segments between pairs of present-day individuals. This study shows that future inbreeding can be accurately predicted, including at specific loci, but not all methods perform equally well.